Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,193,287 T→C 6.9% D97G (GAC→GGC) ‡ thiF ← adenylyltransferase, modifies ThiS C‑terminus

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,2870TC 6.9% 290.6 / 9.9 102D97G (GAC→GGC) ‡thiFadenylyltransferase, modifies ThiS C‑terminus
Reads supporting (aligned to +/- strand):  ref base T (48/45);  new base C (5/2);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

Only 100 of 102 total aligned reads displayed.

AAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCC  >  NC_000913/4193139‑4193433
                                                                                                                                                    |                                                                                                                                                  
aaGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTc                                                                                                                                                   >  2:1350845/1‑150 (MQ=255)
aaGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTc                                                                                                                                                   >  2:1545067/1‑150 (MQ=255)
 aGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCg                                                                                                                                                  >  2:626959/1‑150 (MQ=255)
  gcggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTATTGCTGTCAGTTGAATGTCgg                                                                                                                                                 <  1:238789/150‑1 (MQ=255)
  cctgcgTGTTGAGAGCCACGCAGGCGGCATTAATCTCCTGCCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCgg                                                                                                                                                 <  2:327222/147‑1 (MQ=255)
      cgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAATGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGAtt                                                                                                                                             >  2:1247637/1‑150 (MQ=255)
      cgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGagt                                                                                                                                             >  1:747038/1‑148 (MQ=255)
      cgTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGAtt                                                                                                                                             >  2:514180/1‑150 (MQ=255)
       gtgtTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTc                                                                                                                                            <  2:1034438/150‑1 (MQ=255)
        tgtTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCa                                                                                                                                           >  2:939005/1‑150 (MQ=255)
          ttGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAc                                                                                                                                         <  1:247315/150‑1 (MQ=255)
            gAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACtg                                                                                                                                       <  1:1354292/150‑1 (MQ=255)
             aGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCATCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                                      <  1:1664649/150‑1 (MQ=255)
              gTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgtg                                                                                                                                     >  2:1268414/1‑150 (MQ=255)
               tGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCTCATCGGCCCGGGCAACCGCAGCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTgtgt                                                                                                                                    <  1:802624/150‑1 (MQ=255)
                gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGCGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTGGTTGTAATGCTGTCAGTTGAATGTCGGGATTCACCTGTGTc                                                                                                                                   >  1:118319/1‑150 (MQ=255)
                gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTc                                                                                                                                   >  2:503311/1‑150 (MQ=255)
                gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTAGGGATTAAACTttgta                                                                                                                                   >  1:365506/1‑149 (MQ=255)
                gCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATCCCGACATTCAACtgacag                                                                                                                                   <  1:1244713/150‑5 (MQ=255)
                      cAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAGCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGTATTCAACTGTGTCAGTCGc                                                                                                                             >  2:54126/1‑150 (MQ=255)
                         gcggcATTAATCTCCTGGCGAGTCGCCATATTCTCGGTACAGTCGAGCACCACAACGGCCCGTGCAACCGCATCTTTTAAAGCCTCACACGTTAACCGTTGTTGTAATGCAGTCAGTAGAATGTCGGGATTCAACTGTGTCAGTGGCTGt                                                                                                                          >  2:284586/1‑150 (MQ=255)
                            gcATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTgg                                                                                                                       >  1:1454739/1‑150 (MQ=255)
                               ttAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCAGCTTTTAACGCCTCCCCCCTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTg                                                                                                                    <  2:1081797/150‑1 (MQ=255)
                                 aaTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGATTGTCGGGATTCAACTGTGTc                                                                                                                                   >  1:825957/1‑133 (MQ=255)
                                 aaTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTc                                                                                                                                   <  2:825957/133‑1 (MQ=255)
                                 aaTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAc                                                                                                                  <  2:698774/150‑1 (MQ=255)
                                  aTCTCCTGGCGAGTCTCCATATTGTCGGTACAGTCTAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCAGCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAct                                                                                                                 <  1:284586/150‑2 (MQ=255)
                                  aTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                                 >  1:1320321/1‑150 (MQ=255)
                                  aTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                                 <  2:34805/150‑1 (MQ=255)
                                  aTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                                 <  2:370969/150‑1 (MQ=255)
                                  aTCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                                 <  1:626959/150‑1 (MQ=255)
                                    ctcCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTg                                                                                                               >  2:880457/1‑150 (MQ=255)
                                      ccTTGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGGCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCg                                                                                                             <  1:1082358/150‑1 (MQ=255)
                                      ccTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGGTGTAATGCTGTCAGTGGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCg                                                                                                             >  1:1424175/1‑150 (MQ=255)
                                               gcGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                                    <  2:224626/149‑1 (MQ=255)
                                                   cATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                                >  1:537549/1‑150 (MQ=255)
                                                   cATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                                <  1:122795/150‑1 (MQ=255)
                                                        tGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGatat                                                                                           >  1:826984/1‑150 (MQ=255)
                                                         gTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACACGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATc                                                                                          >  1:784473/1‑150 (MQ=255)
                                                          tCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCt                                                                                         <  2:1336411/150‑1 (MQ=255)
                                                            ggTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGGCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTc                                                                                       <  1:1101581/150‑1 (MQ=255)
                                                            ggTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCATCACCCGGTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTTTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTc                                                                                       >  1:856220/1‑150 (MQ=255)
                                                             gTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCa                                                                                      <  2:491484/150‑1 (MQ=255)
                                                                  gTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGt                                                                                 >  2:714143/1‑150 (MQ=255)
                                                                  gTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGt                                                                                 >  1:427461/1‑150 (MQ=255)
                                                                        caccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAagag                                                                           >  2:339067/1‑150 (MQ=255)
                                                                         accacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTATTGCTGTCAGTTGAATGTCGGGCTTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                                          >  1:1581004/1‑150 (MQ=255)
                                                                          ccacATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGa                                                                         >  2:799728/1‑150 (MQ=255)
                                                                            acaTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGTCGATCGATATCTTCAGTGGTAAAGAGGAtt                                                                       <  2:736714/150‑1 (MQ=255)
                                                                              aTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTg                                                                     <  1:1725835/150‑1 (MQ=255)
                                                                              aTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTg                                                                     <  2:1320321/150‑1 (MQ=255)
                                                                                cGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTc                                                                   <  2:1538441/150‑1 (MQ=255)
                                                                                 ggCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCgg                                                                                                      <  2:912844/114‑1 (MQ=255)
                                                                                 ggCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCgg                                                                                                      >  1:912844/1‑114 (MQ=255)
                                                                                 ggCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCg                                                                  <  2:397255/150‑1 (MQ=255)
                                                                                  gCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGACTGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGt                                                                 <  2:926133/150‑1 (MQ=255)
                                                                                   cccGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATTTCTTCAGTGGTAAAGAGGATTTGTCGtt                                                                >  1:1549801/1‑150 (MQ=255)
                                                                                     cGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGc                                                              >  2:255737/1‑150 (MQ=255)
                                                                                      gTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCa                                                             <  2:1139597/150‑1 (MQ=255)
                                                                                        gCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGa                                                           <  1:1158369/150‑1 (MQ=255)
                                                                                              gCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCtt                                                     <  1:1214086/150‑1 (MQ=255)
                                                                                               cATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                                    <  1:1438237/150‑1 (MQ=255)
                                                                                               cATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGATGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                                    >  1:1553391/1‑150 (MQ=255)
                                                                                                   ttttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                <  1:1440258/150‑1 (MQ=255)
                                                                                                   ttttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                                >  1:150914/1‑150 (MQ=255)
                                                                                                     ttAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGca                                              <  1:604625/150‑1 (MQ=255)
                                                                                                      tAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                                             <  2:1442055/150‑1 (MQ=255)
                                                                                                       aaCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcaca                                            >  2:716259/1‑150 (MQ=255)
                                                                                                              cACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATcgtcgt                                     <  1:541406/150‑1 (MQ=255)
                                                                                                                  cGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                                 >  2:1494125/1‑150 (MQ=255)
                                                                                                                    ttAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                               <  1:1162873/150‑1 (MQ=255)
                                                                                                                    ttAACCGCTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                               <  2:645739/150‑1 (MQ=255)
                                                                                                                      aaCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGcc                             <  2:1100636/150‑1 (MQ=255)
                                                                                                                        ccGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAg                           >  1:1172335/1‑150 (MQ=255)
                                                                                                                        ccGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAg                           <  1:832547/150‑1 (MQ=255)
                                                                                                                          gttgttGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTa                         >  2:1660833/1‑150 (MQ=255)
                                                                                                                             gttgtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCa                      >  2:413263/1‑150 (MQ=255)
                                                                                                                              ttgtACTGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAg                     <  1:1247637/150‑1 (MQ=255)
                                                                                                                                gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                   >  2:523389/1‑150 (MQ=255)
                                                                                                                                gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                   >  1:1085479/1‑150 (MQ=255)
                                                                                                                                gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                   <  1:1350845/150‑1 (MQ=255)
                                                                                                                                gtAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg                   <  1:1545067/150‑1 (MQ=255)
                                                                                                                                  aaTGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTc                 >  2:891873/1‑150 (MQ=255)
                                                                                                                                  aaTGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTc                 <  2:1561399/150‑1 (MQ=255)
                                                                                                                                     gCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCg              <  2:1454739/150‑1 (MQ=255)
                                                                                                                                     gCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCg              <  1:880457/150‑1 (MQ=255)
                                                                                                                                       acacAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGAc            >  2:782335/4‑150 (MQ=255)
                                                                                                                                        cacAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg           >  1:1305784/3‑150 (MQ=255)
                                                                                                                                        cacAGTTGAATCCCGACATCCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACg           >  1:424972/3‑150 (MQ=255)
                                                                                                                                          cAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGcc         <  2:1227328/150‑1 (MQ=255)
                                                                                                                                          cAGTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGcc         >  2:1661428/1‑150 (MQ=255)
                                                                                                                                           aGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGccc        >  2:1463179/1‑149 (MQ=255)
                                                                                                                                           aGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCa        >  2:813581/1‑150 (MQ=255)
                                                                                                                                           aGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGAACATCGTCGTCATCTGCAAGTACCAGCGTCCCGACGCCa        >  2:231164/1‑150 (MQ=255)
                                                                                                                                            gTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCTGACGCCAg       >  2:1223663/1‑150 (MQ=255)
                                                                                                                                            gTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAg       >  1:815303/1‑150 (MQ=255)
                                                                                                                                            gTTGAATCCCGACATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAg       >  2:1671730/1‑150 (MQ=255)
                                                                                                                                                aaTGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGcgcc   <  2:537549/150‑1 (MQ=255)
                                                                                                                                                aaTGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGcgcc   >  2:448187/1‑150 (MQ=255)
                                                                                                                                                 aTGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGcgccc  >  2:188733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCC  >  NC_000913/4193139‑4193433

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