Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,314,877 C→T 6.1% D380N (GAC→AAC) ‡ infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,314,8770CT 6.1% 495.4 / 15.8 147D380N (GAC→AAC) ‡infBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base C (66/72);  new base T (5/4);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 148 total aligned reads displayed.

CAATGTGCTGGGTAATGCCGCCCGCTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTTCATCATCGCTTTGATGACCTGA  >  NC_000913/3314730‑3315025
                                                                                                                                                   |                                                                                                                                                    
cAATGTGCTGGGTAATGCCGCCCGCTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCg                                                                                                                                                     >  2:625333/1‑149 (MQ=255)
              aTGCCGCCCGCTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCctc                                                                                                                                       <  3:1145815/149‑1 (MQ=255)
               tGCCGCCCGCTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCctct                                                                                                                                      <  1:625333/149‑1 (MQ=255)
                 ccgcccgcTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTc                                                                                                                                    <  4:253757/149‑1 (MQ=255)
                       gcTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTc                                                                                                                              <  4:519678/149‑1 (MQ=255)
                         ttCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGc                                                                                                                                                    >  3:1430725/1‑125 (MQ=255)
                             ccAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCAtt                                                                                                                                               >  2:1128218/1‑126 (MQ=255)
                             ccAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCAtt                                                                                                                                               <  1:1128218/126‑1 (MQ=255)
                             ccAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTc                                                                                                                        >  2:527433/1‑149 (MQ=255)
                               agagGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCac                                                                                                                      >  4:1343012/1‑149 (MQ=255)
                                        tttCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGa                                                                                                             <  4:932613/149‑1 (MQ=255)
                                        tttCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGa                                                                                                             >  2:1443946/1‑149 (MQ=255)
                                        tttCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGa                                                                                                             >  1:529787/1‑149 (MQ=255)
                                        tttCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGa                                                                                                             >  3:1151768/1‑149 (MQ=255)
                                         ttCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGc                                                                                                                                                    >  1:365131/1‑109 (MQ=255)
                                            gTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAAc                                                                                                         <  3:1064290/149‑1 (MQ=255)
                                              tGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACtt                                                                                                       <  4:462849/149‑1 (MQ=255)
                                                aaCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTc                                                                                                                                                  >  1:464031/1‑104 (MQ=255)
                                                aaCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTc                                                                                                                                                  <  2:464031/104‑1 (MQ=255)
                                                aaCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCAt                                                                                                                                                <  3:880578/106‑1 (MQ=255)
                                                aaCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCAt                                                                                                                                                >  4:880578/1‑106 (MQ=255)
                                                 aCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCAtt                                                                                                                                               >  3:983378/1‑106 (MQ=255)
                                                 aCGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTAt                                                                                                    >  1:112597/1‑149 (MQ=255)
                                                  cGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATg                                                                                                   <  4:1368963/149‑1 (MQ=255)
                                                  cGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATg                                                                                                   <  3:1343012/149‑1 (MQ=255)
                                                   gAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACtt                                                                                                       >  3:612252/1‑144 (MQ=255)
                                                   gAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACtt                                                                                                       <  4:612252/144‑1 (MQ=255)
                                                   gAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATgg                                                                                                  >  2:526141/1‑149 (MQ=255)
                                                     aTGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTAt                                                                                                    <  4:477030/145‑1 (MQ=255)
                                                     aTGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTAt                                                                                                    >  3:477030/1‑145 (MQ=255)
                                                      tGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGccc                                                                                               <  1:1424459/149‑1 (MQ=255)
                                                          gTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCAtct                                                                                           <  1:122266/149‑1 (MQ=255)
                                                          gTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCAtct                                                                                           >  1:906494/1‑149 (MQ=255)
                                                           tCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAg                                                                                                               <  4:743270/128‑1 (MQ=255)
                                                           tCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAg                                                                                                               >  3:743270/1‑128 (MQ=255)
                                                              agcagAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCctct                                                                                                                                      >  4:485518/1‑102 (MQ=255)
                                                                   agagGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGttt                                                                                                                          <  3:488734/109‑1 (MQ=255)
                                                                    gagGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAAc                                                                                 <  1:1443946/149‑1 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAg                                                                                                                                 <  1:784567/97‑1 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAg                                                                                                                                 >  2:784567/1‑97 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGtt                                                                                                                           <  3:422889/103‑1 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGtt                                                                                                                           >  4:422889/1‑103 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCAtctc                                                                                          <  3:1272094/136‑1 (MQ=255)
                                                                        ttttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCAtctc                                                                                          >  4:1272094/1‑136 (MQ=255)
                                                                         tttACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGATCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCt                                                                            <  4:129768/149‑1 (MQ=255)
                                                                             ccGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAATTTTATGGCCCa                                                                                              <  3:836241/127‑1 (MQ=255)
                                                                             ccGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCa                                                                                              >  4:836241/1‑127 (MQ=255)
                                                                             ccGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGc                                                                        <  4:157484/149‑1 (MQ=255)
                                                                                   tCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCg                                                                       >  1:1365559/1‑144 (MQ=255)
                                                                                     aaCGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCacg                                                                                                                     >  4:1224394/1‑96 (MQ=255)
                                                                                     aaCGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCt                                                                <  2:112597/149‑1 (MQ=255)
                                                                                      aCGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTc                                                                                                                                    >  3:15925/1‑80 (MQ=255)
                                                                                      aCGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTg                                                               >  3:156968/1‑149 (MQ=255)
                                                                                         tGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATc                                                            <  1:527433/149‑1 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCtt                                                                                                                                     <  3:39932/72‑1 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTcc                                                                                                                                   <  1:482003/74‑1 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGttt                                                                                                                          >  4:391026/1‑83 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCa                                                                                              >  4:532678/1‑111 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCa                                                                                              <  3:532678/111‑1 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATa                                                        >  3:916573/1‑149 (MQ=255)
                                                                                             ccATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATa                                                        >  3:876690/1‑149 (MQ=255)
                                                                                              cATGATTGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTa                                                                                                     >  3:1388778/1‑103 (MQ=255)
                                                                                              cATGATTGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCGTCCAGCTCGTTTTCACGACGCAGGATAACTTTa                                                                                                     <  4:1388778/103‑1 (MQ=255)
                                                                                               atgatgGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATaa                                                                                                          <  2:205097/97‑1 (MQ=255)
                                                                                               atgatgGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGAt                                                         >  2:1421920/1‑146 (MQ=255)
                                                                                                tgatgGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATg                                                                                                   >  4:1442281/1‑103 (MQ=255)
                                                                                                    ggTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGt                                                 <  4:163236/149‑1 (MQ=255)
                                                                                                     gTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGtt                                                <  3:901222/149‑1 (MQ=255)
                                                                                                      tCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTg                                               >  3:7887/1‑149 (MQ=255)
                                                                                                      tCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTg                                               >  3:1243032/1‑149 (MQ=255)
                                                                                                        acaACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGttt                                                                                                                          >  3:169202/1‑72 (MQ=255)
                                                                                                        acaACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGttt                                                                                                                          <  4:169202/72‑1 (MQ=255)
                                                                                                        acaACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGAt                                             >  1:1372709/1‑149 (MQ=255)
                                                                                                            ccGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGtt                                         <  1:333007/149‑1 (MQ=255)
                                                                                                              ggCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGc                                       <  3:1329685/149‑1 (MQ=255)
                                                                                                               gcgcgcgcGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGcc                                      <  4:7887/149‑1 (MQ=255)
                                                                                                                         ttCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGccc                             <  2:599429/148‑1 (MQ=255)
                                                                                                                         ttCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGccc                             >  1:599429/1‑148 (MQ=255)
                                                                                                                         ttCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCa                            >  3:245100/1‑149 (MQ=255)
                                                                                                                         ttCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCa                            >  3:325269/1‑149 (MQ=255)
                                                                                                                          tCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTcc                                                                 <  4:689411/111‑1 (MQ=255)
                                                                                                                          tCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTcc                                                                 >  3:689411/1‑111 (MQ=255)
                                                                                                                             gCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATaa                                                       <  1:104214/118‑1 (MQ=255)
                                                                                                                              ccGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTtc                       <  4:1151768/149‑1 (MQ=255)
                                                                                                                                 gagcgaccGTGACACGGGTGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTgcgc                               <  3:830046/132‑1 (MQ=255)
                                                                                                                                 gagcgaccGTGACACGGGTGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTgcgc                               >  4:830046/7‑138 (MQ=255)
                                                                                                                                   gcgaccGTGACACGGGTGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGtt                                         >  2:37971/5‑126 (MQ=255)
                                                                                                                                   gcgaccGTGACACGGGTGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTtcatca                   <  2:789704/144‑1 (MQ=255)
                                                                                                                                     gaccGTGACACGGGTGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTc                                                                  <  4:1121160/97‑1 (MQ=255)
                                                                                                                                     aCCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTc                                                                                       <  4:887268/78‑1 (MQ=255)
                                                                                                                                      cccGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTAt                                                                                                    >  3:468212/1‑64 (MQ=255)
                                                                                                                                      cccGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTg                                                   >  2:944518/1‑113 (MQ=255)
                                                                                                                                      cccGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTg                                                   <  1:944518/113‑1 (MQ=255)
                                                                                                                                      cccGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTg                                                   >  2:1357162/1‑113 (MQ=255)
                                                                                                                                         gtgtCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTTCATCATCGCttt            <  4:916573/149‑1 (MQ=255)
                                                                                                                                           gtCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTTCATCATCGCTTtga          >  1:77060/1‑149 (MQ=255)
                                                                                                                                             cACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTgc                                 <  2:560262/124‑1 (MQ=255)
                                                                                                                                              aCGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTc                                                                  >  2:1297174/1‑90 (MQ=255)
                                                                                                                                              aCGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTc                                                                  <  1:1297174/90‑1 (MQ=255)
                                                                                                                                              aCGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTTCATCATCGCTTtgatga       >  2:579821/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAATGTGCTGGGTAATGCCGCCCGCTTCGCCAGAGGCCACTTTCGTTGAACGAATGTAGTCCAGCAGAGAGGTTTTACCGTGGTCAACGTGACCCATGATGGTCACAACCGGCGCGCGCGGTTCAGCCGCAGCACCCGTGTCACGGTCGCTCATTACCGCCTCTTCCAGCTCGTTTTCACGACGCAGGATAACTTTATGGCCCATCTCTTCAGCAACCAGCTGTGCGGTTTCCTGATCGATAACCTGGTTGATGGTTGCCATTGCGCCCAGTTTCATCATCGCTTTGATGACCTGA  >  NC_000913/3314730‑3315025

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