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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,3910AT 6.9% 296.5 / 10.4 87K196I (AAA→ATA) ‡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 = 8.43e-01

CACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTT  >  NC_000913/1895243‑1895537
                                                                                                                                                    |                                                                                                                                                  
caccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGaaaa                                                                                                                                                    >  3:448942/1‑149 (MQ=255)
 accaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGaaaa                                                                                                                                                    <  2:915179/148‑1 (MQ=255)
    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                                                                        >  4:156076/1‑136 (MQ=255)
                                                                                         aGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGAcc                                                                        <  3:156076/136‑1 (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                                                            <  3:658616/133‑1 (MQ=255)
                                                                                                        ttCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                            >  4:658616/1‑133 (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                                                         >  1:544676/1‑123 (MQ=255)
                                                                                                                     tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGaa                                                         <  2:544676/123‑1 (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                                                                    <  1:399541/91‑1 (MQ=255)
                                                                                                                                          tACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTAtt                                                                    >  2:399541/1‑91 (MQ=255)
                                                                                                                                           aCGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgc         <  4:504710/149‑1 (MQ=255)
                                                                                                                                            tggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGc                                                                                           >  2:716714/2‑66 (MQ=255)
                                                                                                                                            tggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGc                                                                                           <  1:716714/65‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
CACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTT  >  NC_000913/1895243‑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: