Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,092,774:1 +A 5.1% coding (897/951 nt) gshB → glutathione synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,092,7741.A 5.1% 504.5 / 22.2 118coding (897/951 nt)gshBglutathione synthetase
Reads supporting (aligned to +/- strand):  ref base . (54/58);  new base A (4/2);  total (58/60)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 122 total aligned reads displayed.

GTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092626‑3092920
                                                                                                                                                     |                                                                                                                                                  
gTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtg                                                                                                                                                     >  3:678264/1‑149 (MQ=255)
  gACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTc                                                                                                                                                  <  3:412063/149‑1 (MQ=255)
           aTCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCgg                                                                                                                                         >  2:685424/1‑149 (MQ=255)
             cGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGaa                                                                                                                                       <  1:537278/149‑1 (MQ=255)
              gCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAAt                                                                                                                                      <  1:840674/149‑1 (MQ=255)
                ccGTCAGATCGCGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGt                                                                                                                                    <  4:634139/149‑1 (MQ=255)
                 cGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGtt                                                                                                                                   >  1:11685/1‑149 (MQ=255)
                 cGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGtt                                                                                                                                   <  4:193041/149‑1 (MQ=255)
                   tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTaa                                                                                                                                 >  4:940184/1‑149 (MQ=255)
                   tCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTaa                                                                                                                                 >  2:94509/1‑149 (MQ=255)
                    cAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                                                                                                                                <  4:64741/149‑1 (MQ=255)
                     aGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATg                                                                                                                               >  3:131382/1‑149 (MQ=255)
                       aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGa                                                                                                                             >  1:567272/1‑149 (MQ=255)
                       aTCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGa                                                                                                                             <  2:905707/149‑1 (MQ=255)
                        tCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGAt                                                                                                                            <  4:942852/149‑1 (MQ=255)
                            gCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCa                                                                                                                        >  2:390498/1‑149 (MQ=255)
                             ccGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACTAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCAt                                                                                                                       >  4:244775/1‑149 (MQ=255)
                             ccGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCAt                                                                                                                       >  1:586566/1‑149 (MQ=255)
                               gACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCg                                                                                                                     >  3:682041/1‑149 (MQ=255)
                                    tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGcg                                                                                                                >  1:641832/1‑148 (MQ=255)
                                    tGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCa                                                                                                                >  1:734808/1‑149 (MQ=255)
                                     gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                               >  3:270744/1‑149 (MQ=255)
                                     gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                               >  4:719527/1‑149 (MQ=255)
                                     gaaagaaaAAGGGCTGATTTTTGTTGGTCTGGATAGCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGCGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCAc                                                                                                               <  4:841811/149‑1 (MQ=255)
                                      aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg                                                                                                              >  4:49962/1‑149 (MQ=255)
                                      aaagaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACg                                                                                                              >  4:461215/1‑149 (MQ=255)
                                        agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                                                                                                                                <  4:541775/129‑1 (MQ=255)
                                        agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAAt                                                                                                                                >  3:541775/1‑129 (MQ=255)
                                        agaaaAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGtt                                                                                                            >  4:851366/1‑149 (MQ=255)
                                            aaaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAca                                                                                                        <  2:428626/149‑1 (MQ=255)
                                             aaGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcag                                                                                                       <  2:586566/149‑1 (MQ=255)
                                                ggCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTCAATGGATGCCATCGAAGCACGTTTAcagcag                                                                                                    <  3:220804/149‑1 (MQ=255)
                                                 gCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagc                                                                                                   >  3:581414/1‑149 (MQ=255)
                                                  cTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagca                                                                                                  <  2:691646/149‑1 (MQ=255)
                                                   tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGACACCGGaaa                                                                                                                                      >  4:13779/1‑111 (MQ=255)
                                                   tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGACACCGGaaa                                                                                                                                      <  3:13779/112‑2 (MQ=255)
                                                   tGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTAcagcagcag                                                                                                 >  3:827903/1‑149 (MQ=255)
                                                     aTTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTa                                                                                               >  2:50637/1‑149 (MQ=255)
                                                     aTTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTa                                                                                               >  1:883579/1‑149 (MQ=255)
                                                      tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGAt                                                                                                                            <  3:664158/119‑1 (MQ=255)
                                                      tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCGCGTTTACAGCAGCAGTaa                                                                                              <  2:641832/149‑1 (MQ=255)
                                                      tttttGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGATTCCGGTGTCGATCACCGGAATGTTAATGGAt                                                                                                                            >  4:664158/1‑119 (MQ=255)
                                                            ttgGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCAcc                                                                                        >  1:876614/1‑149 (MQ=255)
                                                               gTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTa                                                                                     >  4:653162/1‑149 (MQ=255)
                                                                 cTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGc                                                                                   <  3:117056/149‑1 (MQ=255)
                                                                      tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                              <  1:768628/149‑1 (MQ=255)
                                                                      tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                              <  4:682041/149‑1 (MQ=255)
                                                                      tatCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGaa                                                                              <  3:1033297/149‑1 (MQ=255)
                                                                        tcatcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAgg                                                                            <  2:11685/149‑1 (MQ=255)
                                                                          cttggCCACCGTCTGACTGTAACTAACGTCACCAGCGCAACCTGTATTCGTGAGATTGCAGCATAGTTTCCGGTGACGATCACCGGTATGTAAGGGGAGGCGATCTAAGAAGGTTTACAGCAGCAGTAAACTAACGTAGCCAGAAGGAt                                                                          <  2:470053/146‑1 (MQ=255)
                                                                          atcGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAt                                                                          <  4:803707/149‑1 (MQ=255)
                                                                               cGACCGTCTGACTGAAATTAACGTCCCCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGt                                                                     <  3:1059706/149‑1 (MQ=255)
                                                                                gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  4:413193/149‑1 (MQ=255)
                                                                                gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  2:77123/149‑1 (MQ=255)
                                                                                gACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  1:590434/149‑1 (MQ=255)
                                                                                  ccGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTga                                                                  <  3:867026/149‑1 (MQ=255)
                                                                                    gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                <  4:613709/149‑1 (MQ=255)
                                                                                    gTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTgaga                                                                >  1:525728/1‑149 (MQ=255)
                                                                                     tCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtgt                                                                                                                                                   >  1:272279/1‑65 (MQ=255)
                                                                                     tCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGgtgt                                                                                                                                                   <  2:272279/65‑1 (MQ=255)
                                                                                       tgactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGActc                                                             <  2:540385/149‑1 (MQ=255)
                                                                                       tgactgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGActc                                                             <  4:679805/149‑1 (MQ=255)
                                                                                          ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGa                                                          <  4:830348/149‑1 (MQ=255)
                                                                                          ctgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTAGAGACTCTGa                                                          <  4:681912/149‑1 (MQ=255)
                                                                                           tgaAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCa                                                                                          >  3:831117/1‑116 (MQ=255)
                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTaa                                                                                                                                 <  3:132097/73‑1 (MQ=255)
                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTaa                                                                                                                                 >  4:132097/1‑73 (MQ=255)
                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtc                                                                         <  3:466176/129‑1 (MQ=255)
                                                                                               aTTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGAtc                                                                         >  4:466176/1‑129 (MQ=255)
                                                                                                    cGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAgcgc                                                >  3:509423/1‑149 (MQ=255)
                                                                                                     gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGACACCGGaaa                                                                                                                                      >  1:1049349/1‑61 (MQ=37)
                                                                                                     gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGATCGACACCGGaaa                                                                                                                                      <  2:1049349/62‑2 (MQ=37)
                                                                                                     gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGcc                                               >  2:979471/1‑149 (MQ=255)
                                                                                                     gTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGATGAGACTCTGAGTGACAGCGcc                                               >  3:662687/1‑149 (MQ=255)
                                                                                                           aGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGTATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCCGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCcttctt                                         >  1:494922/1‑149 (MQ=255)
                                                                                                            gCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCttt                                        >  1:531356/1‑149 (MQ=255)
                                                                                                              ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                      >  2:1055748/1‑149 (MQ=255)
                                                                                                              ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                      >  1:141379/1‑149 (MQ=255)
                                                                                                              ccAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTcc                                      <  4:761461/149‑1 (MQ=255)
                                                                                                                aaCCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCAc                                    <  1:794274/149‑1 (MQ=255)
                                                                                                                   cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                 >  2:540109/1‑149 (MQ=255)
                                                                                                                   cTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCa                                 >  1:216275/1‑149 (MQ=255)
                                                                                                                    tGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCa                                     >  3:69160/1‑144 (MQ=255)
                                                                                                                     gTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATa                               >  4:639430/1‑149 (MQ=255)
                                                                                                                      tATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATAc                              >  1:843315/1‑149 (MQ=255)
                                                                                                                        ttCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTg                            >  3:285467/1‑149 (MQ=255)
                                                                                                                        ttCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTg                            >  1:792334/1‑149 (MQ=255)
                                                                                                                                 ttGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    >  1:743617/1‑100 (MQ=255)
                                                                                                                                 ttGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGtt                                                                    <  2:743617/100‑1 (MQ=255)
                                                                                                                                 ttGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGtt                   >  1:468035/1‑149 (MQ=255)
                                                                                                                                  tGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTg                  >  3:434603/1‑149 (MQ=255)
                                                                                                                                    aaGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCt                <  3:461215/149‑1 (MQ=255)
                                                                                                                                         gagTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATGTCGTTGAGACTCAGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCtttttt           >  2:903143/1‑149 (MQ=255)
                                                                                                                                            tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        <  4:827903/149‑1 (MQ=255)
                                                                                                                                            tttCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  3:988820/1‑149 (MQ=255)
                                                                                                                                            attCCGGTGATCGACACCGGAATGTTAATGGATGTCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  3:537254/2‑149 (MQ=255)
                                                                                                                                            attCCGGTGATCGACACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGaa        >  2:634177/2‑149 (MQ=255)
                                                                                                                                               ccGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            >  2:958850/1‑142 (MQ=255)
                                                                                                                                               ccGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCttttt            <  1:958850/142‑1 (MQ=255)
                                                                                                                                                  gtgtCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGa  >  3:353592/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTGTTGGTCTGGATATCATCGGCGACCGTCTGACTGAAATTAACGTCACCAGCCCAACCTGTATTCGTGAGATTGAAGCAGAGTTTCCGGTGTCGATCACCGGAATGTTAATGGATGCCATCGAAGCACGTTTACAGCAGCAGTAACCCACCTTAGCGAGAAGGATCTCGTTGAGACTCTGAGTGACAGCGCCCTTCTTTCCACGCATACTGGGCGCTGTTGCTTTTTTGAACCAGGA  >  NC_000913/3092626‑3092920

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: