Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,647,672 G→C 6.2% intergenic (+21/+33) gor → / ← dinQ glutathione reductase (NADPH)/UV inducible membrane toxin DinQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,6720GC 6.2% 382.5 / 12.7 112intergenic (+21/+33)gor/dinQglutathione reductase (NADPH)/UV inducible membrane toxin DinQ
Reads supporting (aligned to +/- strand):  ref base G (51/54);  new base C (4/3);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 121 total aligned reads displayed.

GGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCA  >  NC_000913/3647523‑3647820
                                                                                                                                                     |                                                                                                                                                    
ggCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagtag                                                                                                                                                      <  2:1092748/150‑1 (MQ=255)
  cTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCTAGAAAAAAATACTTCGACAATAGCGTCGTCATTCAGCCTTGGGCGGCAGAAGAGTTGGTGACAATGCGTTTAATGTTAAAGGGCTAATAGTAgtg                                                                                                                                                    <  2:326054/150‑1 (MQ=255)
  cTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAgtg                                                                                                                                                    <  2:781427/150‑1 (MQ=255)
   tttGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAgtgt                                                                                                                                                   <  2:439687/150‑1 (MQ=255)
           ggACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGTCTAAGAGTAGTGTGCTCTTAg                                                                                                                                           >  2:397354/1‑150 (MQ=255)
            gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                          >  1:53483/1‑150 (MQ=255)
            gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                          <  1:834477/150‑1 (MQ=255)
                aaaTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCtt                                                                                                                                      <  2:36875/150‑1 (MQ=255)
                    gTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGCCCtttaac                                                                                                                                  <  1:27101/150‑5 (MQ=255)
                     ttGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGTAACTAAAATTTACTTGGAAAATACCTTCGCAATTCACCATACGTCGGAAGAAGATTTCGTGTCAATTCGTTAAATTTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTa                                                                                                                                 <  2:1339256/150‑1 (MQ=255)
                     ttGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCACCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTa                                                                                                                                 <  1:789146/150‑1 (MQ=255)
                       gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGtgtg                                                                                                                                                  <  2:53229/131‑1 (MQ=255)
                       gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGtgtg                                                                                                                                                  >  1:53229/1‑131 (MQ=255)
                       gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                               >  1:1333260/1‑150 (MQ=255)
                        cAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGtt                                                                                                                             >  1:1035554/1‑150 (MQ=255)
                         aGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGtt                                                                                                                             >  2:1077559/1‑150 (MQ=255)
                            gCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTcc                                                                                                                          >  1:597932/1‑150 (MQ=255)
                             cTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACTATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGGTCTTAGCCCTTAATTACGTTTCCg                                                                                                                         >  1:489750/1‑150 (MQ=255)
                             cTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCg                                                                                                                         >  1:1390143/1‑150 (MQ=255)
                              ttCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGc                                                                                                                        <  2:1196853/150‑1 (MQ=255)
                              ttCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGc                                                                                                                        <  1:1037570/150‑1 (MQ=255)
                                      ggCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGtt                                                                                                                >  1:1241774/1‑150 (MQ=255)
                                       gcgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTc                                                                                                               <  2:781733/150‑1 (MQ=255)
                                       gcgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTc                                                                                                               <  1:1403284/150‑1 (MQ=255)
                                        cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                              >  1:687624/1‑150 (MQ=255)
                                        cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGcc                                                                                                                                         >  1:1161510/1‑123 (MQ=255)
                                        cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGcc                                                                                                                                         <  2:1161510/123‑1 (MQ=255)
                                        cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGGTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                              >  2:71050/1‑150 (MQ=255)
                                         gcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAg                                                                                                             <  2:1557332/150‑1 (MQ=255)
                                         gcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAg                                                                                                             <  2:53483/150‑1 (MQ=255)
                                         gcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAg                                                                                                             <  2:1145938/150‑1 (MQ=255)
                                           tGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGttt                                                                                                                            >  1:1420469/1‑133 (MQ=255)
                                           tGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGttt                                                                                                                            >  1:1414622/1‑133 (MQ=255)
                                           tGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGttt                                                                                                                            <  2:1420469/133‑1 (MQ=255)
                                           tGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGttt                                                                                                                            <  2:1414622/133‑1 (MQ=255)
                                            gAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAg                                                                                                          >  1:470072/1‑150 (MQ=255)
                                             aaGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGc                                                                                                         <  1:90218/150‑1 (MQ=255)
                                             aaGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCt                                                                                                        <  2:1035554/150‑1 (MQ=255)
                                               gATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTg                                                                                                       <  2:121716/150‑1 (MQ=255)
                                               gATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTg                                                                                                       >  2:257470/1‑150 (MQ=255)
                                                aTGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGa                                                                                                      <  1:700292/150‑1 (MQ=255)
                                                 tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCATAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGaa                                                                                                     <  2:633552/150‑1 (MQ=255)
                                                 tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                               <  2:882974/124‑1 (MQ=255)
                                                 tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                               >  1:882974/1‑124 (MQ=255)
                                                   ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  2:1090254/150‑1 (MQ=255)
                                                         aCCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAg                                                                                             >  1:1272401/1‑150 (MQ=255)
                                                         aCCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAg                                                                                             >  1:468389/1‑150 (MQ=255)
                                                          ccAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGc                                                                                            >  1:438696/1‑150 (MQ=255)
                                                            aaaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCgg                                                                                          <  2:357247/150‑1 (MQ=255)
                                                            aaaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCgg                                                                                          <  1:245203/150‑1 (MQ=255)
                                                            aaaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCgg                                                                                          >  1:933445/1‑150 (MQ=255)
                                                               aaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATc                                                                                       <  2:753874/150‑1 (MQ=255)
                                                                aaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCa                                                                                      <  1:311475/150‑1 (MQ=255)
                                                                  gACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGt                                                                                    >  2:472798/1‑150 (MQ=255)
                                                                    cTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTc                                                                                  >  1:1400427/1‑150 (MQ=255)
                                                                      tCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                                <  2:519550/150‑1 (MQ=255)
                                                                      tCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                                <  2:1025641/150‑1 (MQ=255)
                                                                      tCGACAATACCGTCGCCATTCACCAAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                                <  2:392944/150‑1 (MQ=255)
                                                                       cGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcag                                                                               >  1:126703/1‑150 (MQ=255)
                                                                             tACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAgcgc                                                                         <  2:1545191/150‑1 (MQ=255)
                                                                                 gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAAtt                                                                     <  1:965766/150‑1 (MQ=255)
                                                                                    gCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTGAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAAc                                                                  >  1:1445769/1‑150 (MQ=255)
                                                                                    gCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAAc                                                                  >  1:984273/1‑150 (MQ=255)
                                                                                     ccATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACg                                                                 >  1:938/1‑150 (MQ=255)
                                                                                     ccATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACg                                                                 <  2:1272401/150‑1 (MQ=255)
                                                                                      cATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGc                                                                >  2:1608237/1‑150 (MQ=255)
                                                                                       aTTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGcc                                                               <  2:1300081/150‑1 (MQ=255)
                                                                                         tCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  1:190107/112‑1 (MQ=255)
                                                                                         tCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   >  2:190107/1‑112 (MQ=255)
                                                                                         tCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGcccc                                                             >  2:455678/1‑150 (MQ=255)
                                                                                            cccAACGGCGGCAGAAGAGTTCGTGACAATACGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAg                                                          <  2:535141/150‑1 (MQ=255)
                                                                                             ccAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGa                                                         >  1:1657149/1‑150 (MQ=255)
                                                                                                aCGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGa                                                                                                <  2:1467348/108‑1 (MQ=255)
                                                                                                 cggcggCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAACCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACga                                                     <  2:952772/150‑1 (MQ=255)
                                                                                                    cggcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatga                                                  >  2:609491/1‑150 (MQ=255)
                                                                                                     ggcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatgat                                                 >  1:529249/1‑150 (MQ=255)
                                                                                                       cAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTg                                               >  2:1392490/1‑150 (MQ=255)
                                                                                                       cAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTg                                               <  1:1317501/150‑1 (MQ=255)
                                                                                                         gaagaGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCt                                             <  1:1077559/150‑1 (MQ=255)
                                                                                                            gagTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTa                                          >  1:1493276/1‑150 (MQ=255)
                                                                                                              gTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatc                                        <  1:257470/150‑1 (MQ=255)
                                                                                                                tCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaa                                      <  2:470072/150‑1 (MQ=255)
                                                                                                                 cGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaat                                     <  2:933445/150‑1 (MQ=255)
                                                                                                                   tGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaatca                                   >  2:357001/1‑150 (MQ=255)
                                                                                                                   tGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaatca                                   >  2:757264/1‑150 (MQ=255)
                                                                                                                      cAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGATATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCAccc                                >  1:1563442/1‑150 (MQ=255)
                                                                                                                       aaTGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCg                               <  2:597932/150‑1 (MQ=255)
                                                                                                                          gCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGttt                            >  1:1563859/1‑150 (MQ=255)
                                                                                                                                aaTGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCa                      >  2:1355246/1‑150 (MQ=255)
                                                                                                                                   gTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGc                   <  2:612324/150‑1 (MQ=255)
                                                                                                                                   gTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGc                   <  1:543296/150‑1 (MQ=255)
                                                                                                                                        aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              <  2:51815/150‑1 (MQ=255)
                                                                                                                                        aGGGCTAAGAGCACACTACTCTTAGCCCTTAATTAAGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCACTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              >  1:965112/1‑150 (MQ=255)
                                                                                                                                        aGGGCTAAGAGCACACTACTCTTAGCCCTTAAATACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCGGTTTTCTCCATGCGATgg              >  1:1531961/1‑150 (MQ=255)
                                                                                                                                           gCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGt           >  2:840820/1‑150 (MQ=255)
                                                                                                                                               agagTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTgaga       <  2:1057140/150‑1 (MQ=255)
                                                                                                                                                gagTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATCGCATGGAgaaaacg      <  2:583822/150‑7 (MQ=255)
                                                                                                                                                gagTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATCGCATGGAgaaaacg      <  1:609491/150‑7 (MQ=255)
                                                                                                                                                  gtagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATg    <  1:560708/150‑1 (MQ=255)
                                                                                                                                                    agtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCa  >  1:673249/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCA  >  NC_000913/3647523‑3647820

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