Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,883 +C coding (667/687 nt) rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8791.C98.4% 241.4 / ‑3.3 62coding (671/687 nt)rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (26/35);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTC  >  NC_000913/3815741‑3816012
                                                                                                                                             |                                                                                                                                     
cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga‑t‑t‑c                                                                                                                                      <  1:563867/140‑1 (MQ=255)
   tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑cccc                                                                                                                                   >  1:112425/1‑140 (MQ=255)
   tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑cccc                                                                                                                                   >  1:1109533/1‑140 (MQ=255)
     aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑ccccct                                                                                                                                 <  2:1209326/140‑3 (MQ=255)
     aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑ccccct                                                                                                                                 <  2:430548/140‑3 (MQ=255)
          gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGc                                                                                                                             <  1:1344939/140‑1 (MQ=255)
               gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCagag                                                                                                                        <  1:10451/140‑1 (MQ=255)
               gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCagag                                                                                                                        <  1:1209544/140‑1 (MQ=255)
                atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGc                                                                                                                       <  1:706728/140‑1 (MQ=255)
                 tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGcc                                                                                                                      <  1:1336856/140‑1 (MQ=255)
                        tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGa                                                                                                               >  1:375221/1‑140 (MQ=255)
                          aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAGAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTCCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:893438/1‑140 (MQ=255)
                          aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:513846/1‑140 (MQ=255)
                          aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  1:603282/1‑140 (MQ=255)
                          aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  1:409670/1‑140 (MQ=255)
                          aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:408576/1‑140 (MQ=255)
                           tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGa                                                                                                            >  1:1009722/1‑140 (MQ=255)
                             ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                                          <  2:1109533/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCca                                                                                                                     <  1:231809/127‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCca                                                                                                                     >  2:231809/1‑127 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:104215/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  2:938543/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  2:162420/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:787960/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:628935/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:617442/140‑1 (MQ=255)
                               ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:384755/140‑1 (MQ=255)
                                 cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                                                      >  2:980759/1‑140 (MQ=255)
                                  gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctc                                                                                                     >  1:560375/1‑140 (MQ=255)
                                    tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                                   <  1:697289/140‑1 (MQ=255)
                                          tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                                                                                             >  1:492519/1‑140 (MQ=255)
                                            ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                           <  2:1326817/140‑1 (MQ=255)
                                             tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                                          <  1:1168641/140‑1 (MQ=255)
                                              aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                         >  1:776265/1‑140 (MQ=255)
                                              aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCCGCGTCGCTACAATTGGTTCGATT‑T‑CCCCTCGGGCCCGAGCCAACAAGGTGAGGAGGTCTTCATGGGGGaa                                                                                         >  2:195469/1‑140 (MQ=255)
                                                         gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                              <  2:888253/140‑1 (MQ=255)
                                                          aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                             >  2:1137346/1‑140 (MQ=255)
                                                            aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           >  1:383933/1‑140 (MQ=255)
                                                            aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           <  2:771349/140‑1 (MQ=255)
                                                              aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                         >  1:622917/1‑140 (MQ=255)
                                                              aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                         <  2:776265/140‑1 (MQ=255)
                                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 <  2:218106/140‑1 (MQ=255)
                                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTCCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 <  1:893438/140‑1 (MQ=255)
                                                                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                                                               <  1:480578/140‑1 (MQ=255)
                                                                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                    >  1:829731/1‑140 (MQ=255)
                                                                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                    >  2:704796/1‑140 (MQ=255)
                                                                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                    <  1:1121835/140‑1 (MQ=255)
                                                                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                    >  2:221192/1‑140 (MQ=255)
                                                                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                                                 >  1:1191618/1‑140 (MQ=255)
                                                                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                                                 >  2:183796/1‑140 (MQ=255)
                                                                                       tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGcc                                                <  1:1137346/140‑1 (MQ=255)
                                                                                        aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                               <  2:560375/140‑1 (MQ=255)
                                                                                               aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                        >  2:1372835/1‑140 (MQ=255)
                                                                                               aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                        >  1:981639/1‑140 (MQ=255)
                                                                                               aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                        >  2:1106006/1‑140 (MQ=255)
                                                                                                gTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TGC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                       <  1:221192/140‑1 (MQ=255)
                                                                                                 tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                      <  1:488626/140‑1 (MQ=255)
                                                                                                 tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                      <  2:1009722/140‑1 (MQ=255)
                                                                                                  ttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                     <  1:980759/140‑1 (MQ=255)
                                                                                                  ttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                     <  2:671030/140‑1 (MQ=255)
                                                                                                       aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                <  2:603282/140‑1 (MQ=255)
                                                                                                       aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                <  2:829731/140‑1 (MQ=255)
                                                                                                              ccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                         >  2:652520/1‑140 (MQ=255)
                                                                                                                  cTGCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                     <  1:1347925/140‑1 (MQ=255)
                                                                                                                    gCGTCGCTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                   <  2:1302458/140‑1 (MQ=255)
                                                                                                                          cTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct             >  2:452156/1‑140 (MQ=255)
                                                                                                                          cTACAATGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct             >  1:662373/1‑140 (MQ=255)
                                                                                                                                tGGATTCGA‑T‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa       <  2:1191618/140‑1 (MQ=255)
                                                                                                                                  gattcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga     >  2:997981/1‑140 (MQ=255)
                                                                                                                                  gattcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga     <  1:513846/140‑1 (MQ=255)
                                                                                                                                  gattcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga     <  1:1106006/140‑1 (MQ=255)
                                                                                                                                    ttcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                   >  1:989665/1‑124 (MQ=255)
                                                                                                                                    ttcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                   <  2:989665/124‑1 (MQ=255)
                                                                                                                                     tcga‑t‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc  >  1:896626/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                     
CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTC  >  NC_000913/3815741‑3816012

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