Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,033,627 T→A 5.4% intergenic (+14/+30) idi → / ← lysS isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,033,6270TA 5.4% 431.9 / 10.0 147intergenic (+14/+30)idi/lysSisopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (69/70);  new base A (5/3);  total (74/73)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 151 total aligned reads displayed.

TTATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATG  >  NC_000913/3033478‑3033776
                                                                                                                                                      |                                                                                                                                                     
ttATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACAt                                                                                                                                                       <  2:1132050/150‑1 (MQ=255)
 tATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACAtt                                                                                                                                                      >  1:1239574/1‑150 (MQ=255)
  aTCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACAttt                                                                                                                                                     >  1:701498/1‑150 (MQ=255)
    cAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGATCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGc                                                                                                                                                   >  2:1512804/1‑150 (MQ=255)
     taTGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGcc                                                                                                                                                  <  2:388793/149‑1 (MQ=255)
     aaTGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGcc                                                                                                                                                  >  1:819738/1‑150 (MQ=255)
      aTGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAAACCCCGACATTTGcc                                                                                                                                                  <  2:1264311/150‑1 (MQ=255)
        ggTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCggg                                                                                                                                               >  2:28454/1‑150 (MQ=255)
          tgtgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGAAAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGCCTTTTGCCTGGGt                                                                                                                                             >  2:1084166/1‑150 (MQ=255)
            tgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTtg                                                                                                                                           <  1:1461607/150‑1 (MQ=255)
            tgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTtg                                                                                                                                           >  1:1188511/1‑150 (MQ=255)
             taTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTtgt                                                                                                                                          <  1:1595533/149‑1 (MQ=255)
             gATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTtgt                                                                                                                                          >  2:411028/1‑150 (MQ=255)
              aTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGGGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGGCAAATGTCGGGGttttt                                                                                                                                         <  2:1203305/150‑4 (MQ=255)
               tttAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGa                                                                                                                                        >  1:413131/1‑150 (MQ=255)
               tttAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACACCGACATTTGCCGGGGTTGTGa                                                                                                                                        >  2:111199/1‑150 (MQ=255)
                  aGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGGCAAATGTCGGGGtttttttat                                                                                                                                     <  2:1696472/150‑8 (MQ=255)
                    cAGATGTACTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATa                                                                                                                                   <  1:734560/150‑1 (MQ=255)
                     aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATaa                                                                                                                                  >  2:568107/1‑150 (MQ=255)
                     aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATaa                                                                                                                                  <  1:455040/150‑1 (MQ=255)
                     aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATaa                                                                                                                                  >  2:832333/1‑150 (MQ=255)
                      gATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGCCATTTGCCGGGGTTGTGAGCATAAc                                                                                                                                 >  2:940643/1‑150 (MQ=255)
                           aTTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAAt                                                                                                                            <  1:271422/150‑1 (MQ=255)
                           aTTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAAt                                                                                                                            <  1:597355/150‑1 (MQ=255)
                           aTTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTAACCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGAATAACGTAAt                                                                                                                            >  1:411306/1‑150 (MQ=255)
                             tACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGc                                                                                                                          >  2:79419/1‑150 (MQ=255)
                                acGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTa                                                                                                                       >  2:1289327/1‑150 (MQ=255)
                                 cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                      >  1:1581284/1‑150 (MQ=255)
                                 cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                      >  1:877880/1‑150 (MQ=255)
                                 cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                      >  1:681273/1‑150 (MQ=255)
                                   gTATTGATGCCACGCCGTTGGCGTTCAGTCCGTTGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAttt                                                                                                                    <  1:1162545/150‑1 (MQ=255)
                                   gTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAttt                                                                                                                    <  1:731219/150‑1 (MQ=255)
                                       tGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTAcc                                                                                                                <  1:172230/150‑1 (MQ=255)
                                           gCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGAc                                                                                                            >  1:155588/1‑150 (MQ=255)
                                            ccACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGCAGCCAGAAAACGAGTATCTGCATTTACCCAGCTGAAAT‑AAAAAAACCCCCACATTTGCCCGGGTTGTCCGCATAACGGAATGCTAATTTTACCGGCCg                                                                                                           >  1:1044355/1‑150 (MQ=255)
                                            ccACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACg                                                                                                           >  2:670286/1‑150 (MQ=255)
                                            ccACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACg                                                                                                           >  2:895544/1‑150 (MQ=255)
                                             cACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGc                                                                                                          >  1:765589/1‑150 (MQ=255)
                                             cACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGc                                                                                                          >  1:560335/1‑150 (MQ=255)
                                                    tGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCg                                                                                                   >  1:32640/1‑150 (MQ=255)
                                                      ggCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCggg                                                                                                 <  2:188932/150‑1 (MQ=255)
                                                         gTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAAc                                                                                              <  1:329130/150‑1 (MQ=255)
                                                           tCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAg                                                                                            <  2:1297279/150‑1 (MQ=255)
                                                                 cGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAAc                                                                                      <  1:388931/150‑1 (MQ=255)
                                                                 cGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATATTCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGTGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAAc                                                                                      <  2:125053/150‑1 (MQ=255)
                                                                  gTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACg                                                                                     <  1:22021/150‑1 (MQ=255)
                                                                     gATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcg                                                                                  <  1:712399/150‑1 (MQ=255)
                                                                       tGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcgcg                                                                                <  2:219922/150‑1 (MQ=255)
                                                                            aTGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATgg                                                                           <  2:847340/150‑1 (MQ=255)
                                                                                 ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTGTgg                                                                      >  2:223491/1‑150 (MQ=255)
                                                                                 ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCC‑GGAACAGAATAACGTCGCGGATGGTATGGc                                                                     <  1:1050219/150‑1 (MQ=255)
                                                                                  gacaCAAATCGCTCAGCAAGAAATCGATTATCTGCATTTTCCCAGCTTAAAT‑AACAAACCCCCGACATTTGCCGGGGGTGTGAGCATAACTTAATGCTTATTTTACCGGACGCCTCGCCGGGAACAGACTACCGTCGCGGATGGTATGGc                                                                     <  2:289051/147‑1 (MQ=255)
                                                                                      cAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGCCGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                                 >  2:1604779/1‑150 (MQ=255)
                                                                                      cAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                                        <  1:1532038/143‑1 (MQ=255)
                                                                                      cAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                                 <  1:930682/150‑1 (MQ=255)
                                                                                       aaaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                                <  2:1640885/150‑1 (MQ=255)
                                                                                        aaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                               <  2:385932/150‑1 (MQ=255)
                                                                                        aaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                               <  2:155789/150‑1 (MQ=255)
                                                                                         aTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGt                                                              <  2:500371/150‑1 (MQ=255)
                                                                                         aTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGt                                                              <  1:593070/150‑1 (MQ=255)
                                                                                          tCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                             <  2:488609/150‑1 (MQ=255)
                                                                                          tCGCGAAGCCAGAAAACGATTATCAGCATTTACCCAGCTTAAAT‑AAACAACCCCCGACATTTGCCGGGGTTGTGAGCCTAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                             <  2:339374/150‑1 (MQ=255)
                                                                                           cgcgAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGa                                                            >  2:799330/1‑150 (MQ=255)
                                                                                           cgcgAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTAACGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGATGTTGGTGa                                                            >  2:1038326/1‑150 (MQ=255)
                                                                                            gcgAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGAATAACGTAATGCTTATTATACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                                           >  2:201625/1‑150 (MQ=255)
                                                                                             ggAAGGCAGTAAACGTTTTTCTTGATTTTACCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTT‑CCGGACGCATCGCCGGTAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                         <  1:1020466/149‑1 (MQ=255)
                                                                                                aGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGc                                                       >  1:985594/1‑150 (MQ=255)
                                                                                                aGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGc                                                       <  1:832333/150‑1 (MQ=255)
                                                                                                   cAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                    <  2:1691426/150‑1 (MQ=255)
                                                                                                      aaaaCGATTATCTTCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAcc                                                 <  2:701498/150‑1 (MQ=255)
                                                                                                       aaaCGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                >  2:1642879/1‑150 (MQ=255)
                                                                                                          cGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                             >  2:730351/1‑150 (MQ=255)
                                                                                                          cGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                             >  1:1109479/1‑150 (MQ=255)
                                                                                                             ttATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGt                                          <  2:777691/150‑1 (MQ=255)
                                                                                                              tATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                         <  1:565575/150‑1 (MQ=255)
                                                                                                                   gCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                    <  2:32640/150‑1 (MQ=255)
                                                                                                                     aTTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTAGATACCa                                  >  1:1299044/1‑150 (MQ=255)
                                                                                                                       ttACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                                >  2:1230088/1‑150 (MQ=255)
                                                                                                                        tACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                               >  1:963278/1‑150 (MQ=255)
                                                                                                                         aCCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTc                              >  1:1540270/1‑150 (MQ=255)
                                                                                                                          cccAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                             <  2:978019/150‑1 (MQ=255)
                                                                                                                          cccAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACACAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                             <  2:1347746/150‑1 (MQ=255)
                                                                                                                             aGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAg                          >  1:1344943/1‑150 (MQ=255)
                                                                                                                               cTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                        <  1:1508728/150‑1 (MQ=255)
                                                                                                                               cTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                        <  2:1488069/150‑1 (MQ=255)
                                                                                                                                ttAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAcc                       <  1:988429/150‑1 (MQ=255)
                                                                                                                                 taaat‑aaaAAAACCCCGACATTTGCCGGGGTTGAGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGAAGGTATGGCTGATGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAcca                      >  1:1573560/1‑149 (MQ=255)
                                                                                                                                   aat‑aaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAActgc                    >  1:1410609/1‑150 (MQ=255)
                                                                                                                                   aat‑aaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGGGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgc                    <  1:79419/150‑1 (MQ=255)
                                                                                                                                      gctcacaaCCCCGGCAAATGTCGGGGTTTTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgctg                  >  2:1296912/7‑150 (MQ=255)
                                                                                                                                       aaaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt                 >  1:718637/1‑150 (MQ=255)
                                                                                                                                       aaaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt                 >  2:1150535/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTc                <  2:411306/150‑1 (MQ=255)
                                                                                                                                           caaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAAAAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTTTcggc             >  2:1159214/2‑150 (MQ=255)
                                                                                                                                           aaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggc             <  2:59239/150‑1 (MQ=255)
                                                                                                                                           aaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggc             <  1:460017/150‑1 (MQ=255)
                                                                                                                                           aaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggc             <  2:1581284/150‑1 (MQ=255)
                                                                                                                                             aCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGCCGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTACCAGACCTGCTGTcggcgg           >  1:891556/1‑150 (MQ=255)
                                                                                                                                              ccccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt          <  2:155588/150‑1 (MQ=255)
                                                                                                                                                      tttGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg  >  2:1341034/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAAT‑AAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATG  >  NC_000913/3033478‑3033776

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