Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,392 A→T 6.9% K196N (AAA→AAT) ‡ pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,3920AT 6.9% 295.2 / 10.5 87K196N (AAA→AAT) ‡pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base T (3/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

AGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTT  >  NC_000913/1895247‑1895537
                                                                                                                                                 |                                                                                                                                                 
aGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAAtt                                                                                                                                                 >  1:892102/1‑148 (MQ=255)
       tACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAgg                                                                                                                                         >  3:204206/1‑149 (MQ=255)
       tACAGTTTCTTTGCTCAGTCATAATGATGTCAATGCCCGTCGTGCCTGGGTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAACCGAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAg                                                                                                                                          >  2:288840/1‑148 (MQ=255)
        aCAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAgg                                                                                                                                         >  1:435497/1‑148 (MQ=255)
            tttCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACa                                                                                                                                     <  1:607917/148‑1 (MQ=255)
             ttCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                    >  1:698706/1‑148 (MQ=255)
                tttGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGaa                                                                                                                                 <  1:240772/148‑1 (MQ=255)
                    cTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCt                                                                                                                            <  3:118336/149‑1 (MQ=255)
                     tGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCt                                                                                                                            <  2:183582/148‑1 (MQ=255)
                     tGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTg                                                                                                                           >  3:829959/1‑149 (MQ=255)
                      gAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGc                                                                                                                          >  4:9423/1‑149 (MQ=255)
                           aTAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGc                                                                                                                     >  3:43110/1‑149 (MQ=255)
                            tAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGc                                                                                                                     >  2:340285/1‑148 (MQ=255)
                             aaTGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCcaggaacaggaa                                                                                                                                 >  1:239175/1‑135 (MQ=255)
                             aaTGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGaa                                                                                                                                 <  2:239175/135‑1 (MQ=255)
                              atgatgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCgg                                                                                                                   <  2:265335/148‑1 (MQ=255)
                              atgatgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGt                                                                                                                  <  3:790555/149‑1 (MQ=255)
                                 atgTCAATGCCCGTCGGGCCTGTCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                               <  3:586408/149‑1 (MQ=255)
                                 atgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATAACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                               <  3:293095/149‑1 (MQ=255)
                                   gTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAtt                                                                                                              >  2:955096/1‑148 (MQ=255)
                                   gTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAtt                                                                                                              <  2:320284/148‑1 (MQ=255)
                                   gTCAATGCCCGTCGAGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAtt                                                                                                              <  1:500585/148‑1 (MQ=255)
                                    tCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGc                                                                                                            <  3:271807/149‑1 (MQ=255)
                                     cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccgt                                                                                                                                          <  1:178116/118‑3 (MQ=255)
                                     cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccgt                                                                                                                                          >  2:178116/1‑116 (MQ=255)
                                       aTGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                          >  1:257194/1‑148 (MQ=255)
                                         gCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCa                                                                                                       >  3:248379/1‑149 (MQ=255)
                                          cccGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCa                                                                                                       <  2:435497/148‑1 (MQ=255)
                                           ccGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAgg                                                                                                     <  3:144311/149‑1 (MQ=255)
                                            cGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAgg                                                                                                     >  2:238943/1‑148 (MQ=255)
                                                gggCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGAGTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAAt                                                                                                >  4:262139/1‑149 (MQ=255)
                                                     tggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCg                                                                                            <  2:862854/148‑1 (MQ=255)
                                                      ggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGc                                                                                           >  1:342141/1‑148 (MQ=255)
                                                       gctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGcc                                                                                          >  2:657918/1‑148 (MQ=255)
                                                         tggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                           >  4:651256/1‑96 (MQ=255)
                                                         tggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                           <  3:651256/97‑2 (MQ=255)
                                                                ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGtttt                                                                                 <  1:327007/148‑1 (MQ=255)
                                                                 cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGtt                                                                                   <  4:597295/145‑1 (MQ=255)
                                                                 cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGtt                                                                                   >  3:597295/1‑145 (MQ=255)
                                                                   gCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCAt                                                                              >  1:889099/1‑148 (MQ=255)
                                                                       ttgtcGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGGGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATCTCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAc                                                                         <  4:889125/146‑1 (MQ=255)
                                                                       ttCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGa                                                                          <  2:679933/148‑1 (MQ=255)
                                                                          tcGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTa                                                                      <  3:9423/149‑1 (MQ=255)
                                                                                   gAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGa                                                             <  4:521402/149‑1 (MQ=255)
                                                                                    aaGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCGCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                            >  3:778597/1‑149 (MQ=255)
                                                                                     aGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAcc                                                                        <  3:156076/136‑1 (MQ=255)
                                                                                     aGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAcc                                                                        >  4:156076/1‑136 (MQ=255)
                                                                                       aTTTCACGCTCACTTCCGACTGGCAATCCAATATGACCTGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGAGACCTATTCTGGCGATGa                                                          <  2:123267/148‑1 (MQ=255)
                                                                                       aTTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGaa                                                         <  3:776101/149‑1 (MQ=255)
                                                                                       aTTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGaa                                                         <  4:459911/149‑1 (MQ=255)
                                                                                        tttCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATCCTGGCGATGAAt                                                        >  4:516657/1‑149 (MQ=255)
                                                                                            aCGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggca                                                    <  3:113285/149‑1 (MQ=255)
                                                                                             cGCTCACTTCCGACTGGCTATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCAAGGTACAGGAATATCTGCACAGCGGTGATTGCTAACAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggca                                                    >  2:674176/1‑148 (MQ=255)
                                                                                             cGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggca                                                    <  2:400843/148‑1 (MQ=255)
                                                                                              gCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcagg                                                  >  3:899267/1‑149 (MQ=255)
                                                                                                 cACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAc                                                                         >  4:253788/1‑123 (MQ=255)
                                                                                                 cACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAc                                                                         <  3:253788/123‑1 (MQ=255)
                                                                                                    ttCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                            >  4:658616/1‑133 (MQ=255)
                                                                                                    ttCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                            <  3:658616/133‑1 (MQ=255)
                                                                                                      ccGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCt                                           >  1:932998/1‑148 (MQ=255)
                                                                                                        gACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTc                                         <  1:936339/148‑1 (MQ=255)
                                                                                                          cTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGc                                      <  4:43110/149‑1 (MQ=255)
                                                                                                            ggCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                     >  2:775061/1‑148 (MQ=255)
                                                                                                             gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                    >  2:43820/1‑148 (MQ=255)
                                                                                                                aTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                 <  1:790888/148‑1 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAGTCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATc                                >  1:282186/1‑148 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGaa                                                         <  2:544676/123‑1 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGaa                                                         >  1:544676/1‑123 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATc                                <  1:504059/148‑1 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATc                                <  2:428236/148‑1 (MQ=255)
                                                                                                                 tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGACGATGAATGGCAGGCATTCCTTCAGCTTAATCa                               >  4:631115/1‑149 (MQ=255)
                                                                                                                  ccAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCa                               >  1:287158/1‑148 (MQ=255)
                                                                                                                    aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAgg                             <  1:272743/148‑1 (MQ=255)
                                                                                                                      tatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCt                                                                       >  1:209797/1‑104 (MQ=255)
                                                                                                                      tatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCt                                                                       <  2:209797/104‑1 (MQ=255)
                                                                                                                        tGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCaa                         >  2:888220/1‑148 (MQ=255)
                                                                                                                            ccGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcgcg                    >  3:557705/1‑149 (MQ=255)
                                                                                                                               cgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCgcgcgc                 >  3:937383/1‑149 (MQ=255)
                                                                                                                                      tACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTAtt                                                                    >  2:399541/1‑91 (MQ=255)
                                                                                                                                      tACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTAtt                                                                    <  1:399541/91‑1 (MQ=255)
                                                                                                                                       aCGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgc         <  4:504710/149‑1 (MQ=255)
                                                                                                                                        tggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGc                                                                                           <  1:716714/65‑1 (MQ=255)
                                                                                                                                        tggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGc                                                                                           >  2:716714/2‑66 (MQ=255)
                                                                                                                                          gCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcgc       >  2:504533/1‑148 (MQ=255)
                                                                                                                                            gAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGGttt    <  3:363467/149‑1 (MQ=255)
                                                                                                                                            gAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCttt    >  3:92124/1‑149 (MQ=255)
                                                                                                                                            gAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtt     <  1:105952/148‑1 (MQ=255)
                                                                                                                                             aaaaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtttt   >  4:4413/1‑149 (MQ=255)
                                                                                                                                             aaaaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtttt   >  3:150436/1‑149 (MQ=255)
                                                                                                                                             aaaaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCttt    >  1:738198/1‑148 (MQ=255)
                                                                                                                                              aaaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCttttt  >  3:877069/1‑149 (MQ=255)
                                                                                                                                               aaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCCGGCCAACCGCGCGCCATTTAGCGCttttt  <  1:775061/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTT  >  NC_000913/1895247‑1895537

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: