Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,183,102 A→C N620H (AAC→CAC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,183,1020AC98.6% 249.0 / ‑4.3 70N620H (AAC→CAC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (36/33);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGA  >  NC_000913_3_pae_pgi/4182958‑4183246
                                                                                                                                                |                                                                                                                                                
aTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCGGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTcc                                                                                                                                             >  3:352970/1‑150 (MQ=255)
 tCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCa                                                                                                                                            >  1:120139/1‑150 (MQ=255)
     ctctctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACtt                                                                                                                                        <  2:383892/150‑1 (MQ=255)
      tctctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAATTACGTTATCGCCCAGGCGCACTCCAACTTg                                                                                                                                       <  4:231820/150‑1 (MQ=255)
       ctctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTgg                                                                                                                                      >  1:585609/1‑150 (MQ=255)
        tctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGa                                                                                                                                     >  4:273045/1‑150 (MQ=255)
        tctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGa                                                                                                                                     >  1:625293/1‑150 (MQ=255)
               gtgtACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATgaagaa                                                                                                                              >  3:273245/1‑150 (MQ=255)
                    cGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATgaag                                                                                                                                <  1:215923/143‑1 (MQ=255)
                    cGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATgaag                                                                                                                                >  2:215923/1‑143 (MQ=255)
                      cacaGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACt                                                                                                                       >  4:28380/1‑150 (MQ=255)
                       acaGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACtt                                                                                                                      >  3:27295/1‑150 (MQ=255)
                        caGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTc                                                                                                                     <  1:39947/150‑1 (MQ=255)
                           aCTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTa                                                                                                                  <  1:17882/150‑1 (MQ=255)
                              aaCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTagaa                                                                                                               <  1:115218/150‑1 (MQ=255)
                                  aaTACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTag                                                                                                                 >  2:499862/1‑144 (MQ=255)
                                  aaTACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTag                                                                                                                 <  1:499862/144‑1 (MQ=255)
                                      cGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGt                                                                                                       <  2:523021/150‑1 (MQ=255)
                                      cGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGt                                                                                                       >  1:2612/1‑150 (MQ=255)
                                       ggCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTa                                                                                                      >  4:370507/1‑150 (MQ=255)
                                           tCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTTGAAGACCTGGTAACtt                                                                                                  <  2:344428/150‑1 (MQ=255)
                                              ttGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGcc                                                                                               <  3:132533/150‑1 (MQ=255)
                                              ttGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGcc                                                                                               >  3:277964/1‑150 (MQ=255)
                                                    cTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCa                                                                                         >  4:405953/1‑150 (MQ=255)
                                                    cTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCa                                                                                         <  3:28380/150‑1 (MQ=255)
                                                        gTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAgg                                                                                     <  2:27364/150‑1 (MQ=255)
                                                         tATCGTTTAGTCAACGACGGTGTTGCCACTCACGAAATTCCTTACCTGCCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCGACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGc                                                                                    <  2:27718/150‑1 (MQ=255)
                                                         tATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAGTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCGCAGGCGCACTCCAACTTGGATGAAGAGGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGc                                                                                    <  4:352970/150‑1 (MQ=255)
                                                             gTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAAt                                                                                >  4:415048/1‑150 (MQ=255)
                                                             gTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAAt                                                                                >  2:299378/1‑150 (MQ=255)
                                                              tAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATc                                                                               <  3:423039/150‑1 (MQ=255)
                                                              tAAAGTGACCGACGGTGGTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATc                                                                               >  3:422499/1‑150 (MQ=255)
                                                                 aGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAg                                                                            <  1:424296/150‑1 (MQ=255)
                                                                  gTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGc                                                                           <  1:234329/150‑1 (MQ=255)
                                                                   tGACCTACGGTGTTGTAACTGACGAAATTACCTACCTGTCTGGTAGCGAAGAAGGCAACTACGTTATCTCCCACGCGCACTCCAACTTGGATGAAGGAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTGGCAAAGGCGAATCCAGCt                                                                          >  2:477331/1‑150 (MQ=255)
                                                                    gACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCtt                                                                         <  4:428105/150‑1 (MQ=255)
                                                                     aCCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTg                                                                        >  3:421079/1‑150 (MQ=255)
                                                                        gACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTc                                                                     <  3:309338/150‑1 (MQ=255)
                                                                               ttgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCcgcg                                                              >  2:281070/1‑150 (MQ=255)
                                                                                tgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGa                                                             <  1:550378/150‑1 (MQ=255)
                                                                                tgtAACTGACGAAATTCACAACCTGTCTGCTATCGAAGACGTCAAGTACGTTATCGCCCCGGCGCACTCCAACATGGATGAAGAAGGCCACTTCGTAGAAAACCTGGTTACTTGCCTTAGAAAAGACTAAGCGCGCTGGTTCATCCGCGa                                                             >  4:260624/1‑150 (MQ=255)
                                                                                 gtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGAc                                                                                                            <  1:440965/102‑1 (MQ=255)
                                                                                 gtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGAc                                                                                                            >  2:440965/1‑102 (MQ=255)
                                                                                 gtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCaaa                                                                                       <  4:26709/123‑1 (MQ=255)
                                                                                 gtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCaaa                                                                                       >  3:26709/1‑123 (MQ=255)
                                                                                   aaCTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCATGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCa                                                          >  3:138688/1‑150 (MQ=255)
                                                                                      tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGa                                                                                                                                     <  3:144684/72‑1 (MQ=255)
                                                                                      tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGa                                                                                                                                     >  4:144684/1‑72 (MQ=255)
                                                                                      tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGt                                                       <  4:27295/150‑1 (MQ=255)
                                                                                      tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGt                                                       >  1:443419/1‑150 (MQ=255)
                                                                                      tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGt                                                       <  4:273245/150‑1 (MQ=255)
                                                                                          gAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGAc                                                   <  4:58584/150‑1 (MQ=255)
                                                                                           aaaTTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACt                                                  <  1:6099/150‑1 (MQ=255)
                                                                                           aaaTTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACt                                                  <  1:160841/150‑1 (MQ=255)
                                                                                                   tACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGAc                                          <  1:349836/150‑1 (MQ=255)
                                                                                                           tGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCAc                                  >  2:609124/1‑150 (MQ=255)
                                                                                                             cTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCAccc                                >  2:94560/1‑150 (MQ=255)
                                                                                                              tATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAc                                                            >  2:68252/1‑121 (MQ=255)
                                                                                                              tATCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAc                                                            <  1:68252/121‑1 (MQ=255)
                                                                                                              tATCGAAGAAGGCAACTACGCTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCca                               >  2:10122/1‑150 (MQ=255)
                                                                                                                tCGAAGAAGGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCcagc                             <  1:281070/150‑1 (MQ=255)
                                                                                                                   aagaagGCAACTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAgg                          >  4:24319/1‑150 (MQ=255)
                                                                                                                           aaCTGCGTTATCGCCCAGGCGCACTCCAACGTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATcc                  <  3:400709/150‑1 (MQ=255)
                                                                                                                           aaCTACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATcc                  <  3:370507/150‑1 (MQ=255)
                                                                                                                              tACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGAc                                          >  1:299586/1‑123 (MQ=255)
                                                                                                                              tACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGAc                                          <  2:299586/123‑1 (MQ=255)
                                                                                                                              tACGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTc               >  2:102218/1‑150 (MQ=255)
                                                                                                                                cGTTATCGCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCgg             >  1:513412/1‑150 (MQ=255)
                                                                                                                                       gCCCAGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTcc      >  1:8762/1‑150 (MQ=255)
                                                                                                                                           aGGCGCACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGa  >  2:7946/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGA  >  NC_000913_3_pae_pgi/4182958‑4183246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: