Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,537,322 T→A 5.3% intergenic (+34/‑63) greB → / → yhgF transcription elongation factor GreB/putative RNA‑binding protein YhgF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,3220TA 5.3% 300.2 / 4.9 94intergenic (+34/‑63)greB/yhgFtranscription elongation factor GreB/putative RNA‑binding protein YhgF
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base A (2/3);  total (45/49)
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.42e-01

CCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGA  >  NC_000913/3537176‑3537470
                                                                                                                                                  |                                                                                                                                                    
ccccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAtttt                                                                                                                                                   >  1:500070/1‑150 (MQ=255)
 cccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAttttt                                                                                                                                                  <  2:896302/150‑1 (MQ=255)
 cccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAttttt                                                                                                                                                  >  1:591377/1‑150 (MQ=255)
  ccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTg                                                                                                                                                 <  1:278820/150‑1 (MQ=255)
  ccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTg                                                                                                                                                 >  1:279934/1‑150 (MQ=255)
    gATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgcc                                                                                                                                               <  2:13528/150‑1 (MQ=255)
     aTGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgcca                                                                                                                                              <  1:1127571/150‑1 (MQ=255)
        gCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccagcc                                                                                                                                           <  2:1164309/150‑1 (MQ=255)
           cgcgCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatc                                                                                                                                        <  2:673618/150‑3 (MQ=255)
            gcgcATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTc                                                                                                                                       <  2:1190779/150‑1 (MQ=255)
             cgcATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTcc                                                                                                                                      <  2:1367093/150‑1 (MQ=255)
             cgcATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatcgg                                                                                                                                      <  1:961803/150‑5 (MQ=255)
                 ttGCGGAAAAAAGAAGGCGGCGATCTGGCGGTGGTGCATACCCCGGCCGGGGAAGCGAGCTGGTAGGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATGccagccagcggcagg                                                                                                                                  <  1:896963/150‑8 (MQ=255)
                  tGCTTAAAAAATAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGtata                                                                                                                                 <  1:1491362/150‑1 (MQ=255)
                  tGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGtata                                                                                                                                 <  2:123620/150‑1 (MQ=255)
                   gCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATaa                                                                                                                                >  2:1332397/1‑150 (MQ=255)
                   gCTGAAAAAAGAAGTCGGAGATCTGGCGGTGGTGAATACCCCTTCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCTTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATaa                                                                                                                                >  2:1337997/1‑150 (MQ=255)
                    cTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGGTGGCATTTTTGGCAGCCAGTCCGTATAAc                                                                                                                               >  2:160376/1‑150 (MQ=255)
                    cTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAAc                                                                                                                               >  1:1109820/1‑150 (MQ=255)
                           aagaagTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATcccc                                                                                                                        >  2:765480/1‑150 (MQ=255)
                                  cGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGAtttt                                                                                                                 >  2:280149/1‑150 (MQ=255)
                                  cGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGAtttt                                                                                                                 <  2:318553/150‑1 (MQ=255)
                                    gCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTg                                                                                                               <  1:765480/150‑1 (MQ=255)
                                    gCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTg                                                                                                               <  1:339353/150‑1 (MQ=255)
                                    gCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTg                                                                                                               <  1:345408/150‑1 (MQ=255)
                                     cGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGa                                                                                                              <  2:52765/150‑1 (MQ=255)
                                       aTCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATc                                                                                                            <  1:584484/150‑1 (MQ=255)
                                           ggCGGTGGTGAATACCCCGGCCGGGGAAGCGTGCTGGTATGTTAATGCTATCGATTAGGTTAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTTATTTTTGATCCGaa                                                                                                        <  2:1177637/150‑1 (MQ=255)
                                             cGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGTCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGaaaa                                                                                                      >  1:66837/1‑150 (MQ=255)
                                             cGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGaaaa                                                                                                      <  1:147801/150‑1 (MQ=255)
                                              ggtggtGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAg                                                                                                     <  1:1374131/150‑1 (MQ=255)
                                                 ggtgAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATg                                                                                                  <  2:1195935/150‑1 (MQ=255)
                                                 ggtgAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATg                                                                                                  >  1:654748/1‑150 (MQ=255)
                                                     gcttccccGGCGGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACt                                                                                              >  1:268746/5‑150 (MQ=255)
                                                       tACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCa                                                                                            >  1:1480992/1‑150 (MQ=255)
                                                         ccccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCaaa                                                                                          >  1:1403250/1‑150 (MQ=255)
                                                          cccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAc                                                                                         <  2:600692/150‑1 (MQ=255)
                                                          cccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAc                                                                                         >  2:1032839/1‑150 (MQ=255)
                                                           ccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAcc                                                                                        >  1:432304/1‑150 (MQ=255)
                                                           ccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAcc                                                                                        <  2:204785/150‑1 (MQ=255)
                                                             tGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCat                                                                                      <  2:500070/150‑1 (MQ=255)
                                                                 ggggAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatgatg                                                                                  >  2:1066387/1‑150 (MQ=255)
                                                                 ggggAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatgatg                                                                                  <  1:238494/150‑1 (MQ=255)
                                                                         gAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGGATGATTc                                                                          >  1:891311/1‑150 (MQ=255)
                                                                          aGCTGGTATGTTAATGCTCTCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCg                                                                         <  2:1026424/150‑1 (MQ=255)
                                                                             tGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGa                                                                                                    >  2:454090/1‑120 (MQ=255)
                                                                             tGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGa                                                                                                    <  1:454090/120‑1 (MQ=255)
                                                                             tGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTc                                                                      >  1:1389537/1‑150 (MQ=255)
                                                                              ggTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGatgaat                                                                               <  1:1309418/140‑1 (MQ=255)
                                                                              ggTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGatgaat                                                                               >  2:1309418/1‑140 (MQ=255)
                                                                               gTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATTATGAATGATTCGTTATg                                                                    >  1:1371264/1‑150 (MQ=255)
                                                                                tATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGc                                                                   >  1:1177803/1‑150 (MQ=255)
                                                                                 aTGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGcc                                                                  <  2:591377/150‑1 (MQ=255)
                                                                                 aTGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGcc                                                                  <  1:324944/150‑1 (MQ=255)
                                                                                  tGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTTCCAGCCAGTCCGTATAACTATACCCGGATTTTTTATCCGAAAAGATGAACTCAAACCATGATGAATGATTAGTTCTGCCg                                                                 >  1:400531/1‑150 (MQ=255)
                                                                                      aaTGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCatt                                                             >  2:1409309/1‑150 (MQ=255)
                                                                                       aTGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCatta                                                            >  2:644198/1‑150 (MQ=255)
                                                                                          cTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTTCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTg                                                         >  1:202629/1‑150 (MQ=255)
                                                                                           tATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGc                                                        >  1:427801/1‑150 (MQ=255)
                                                                                           tATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGc                                                        <  1:36041/150‑1 (MQ=255)
                                                                                             tCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCgg                                                      <  2:958896/150‑1 (MQ=255)
                                                                                              cGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATCAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCggg                                                     <  1:763689/150‑1 (MQ=255)
                                                                                               gAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGt                                                    >  1:489305/1‑150 (MQ=255)
                                                                                                aGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTg                                                   >  1:767238/1‑150 (MQ=255)
                                                                                                aGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTg                                                   <  2:1177803/150‑1 (MQ=255)
                                                                                                 gTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGa                                                  >  1:525482/1‑150 (MQ=255)
                                                                                                  tACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGaa                                                 <  2:637304/150‑1 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                               <  2:432304/150‑1 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                               >  2:105742/1‑150 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                               >  2:1358443/1‑150 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATACCCTCATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                               >  2:1181556/1‑150 (MQ=255)
                                                                                                     gTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTTATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAtt                                              >  1:397360/1‑150 (MQ=255)
                                                                                                      tGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTc                                             >  1:601943/1‑150 (MQ=255)
                                                                                                      tGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTc                                             <  2:489305/150‑1 (MQ=255)
                                                                                                       gAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCa                                            >  2:871894/1‑150 (MQ=255)
                                                                                                        aaaCCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAg                                           <  2:1480992/150‑1 (MQ=255)
                                                                                                         aaCCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCGTGCATTTTTGCCATCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAgg                                          <  2:982193/150‑1 (MQ=255)
                                                                                                           ccGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcg                                        >  2:994687/1‑150 (MQ=255)
                                                                                                             cgaaGGACCAGTCGTAACCTCCTGCCGATGGCGGGCATTTTTGCCAGCCAGTCCCTATAACTATCCCCTGATTTATGATCCGCAACGAGGAACTCAAACCATCCTGCATGATTCGGTCTGCCGCAGTATTGCGGGTGAAATTCAGgcgcg                                      <  2:1448774/148‑1 (MQ=255)
                                                                                                                 ggAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCg                                  >  2:287702/1‑150 (MQ=255)
                                                                                                                   aaGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGa                                <  2:409079/150‑1 (MQ=255)
                                                                                                                       gTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAg                            >  2:503010/1‑150 (MQ=255)
                                                                                                                        tCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTAAGGCGCGCCCGGAACAgg                           >  1:188486/1‑150 (MQ=255)
                                                                                                                         cTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGt                          >  2:581724/1‑150 (MQ=255)
                                                                                                                         cTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGt                          >  2:737186/1‑150 (MQ=255)
                                                                                                                          ttAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGtt                         <  1:798044/150‑1 (MQ=255)
                                                                                                                                  ctGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGc                 <  1:1066387/150‑1 (MQ=255)
                                                                                                                                  ctGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGc                 <  1:644198/150‑1 (MQ=255)
                                                                                                                                       gATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTc            <  2:654748/150‑1 (MQ=255)
                                                                                                                                       actggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTc            >  2:408921/3‑150 (MQ=255)
                                                                                                                                        ctggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCCGAACAGGTTGACGCTGCCGTTCg           >  2:436965/2‑150 (MQ=255)
                                                                                                                                          ggctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGcc         >  2:1535816/1‑150 (MQ=255)
                                                                                                                                          ggctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGcc         <  2:1290499/150‑1 (MQ=255)
                                                                                                                                            ctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTg       <  1:1605768/150‑1 (MQ=255)
                                                                                                                                            ctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTg       >  2:114137/1‑150 (MQ=255)
                                                                                                                                                 aTTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGa  <  1:1337997/150‑1 (MQ=255)
                                                                                                                                                 aTTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGa  <  1:1332397/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGA  >  NC_000913/3537176‑3537470

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