Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,076,604 G→A 5.7% L153L (CTG→CTA) ‡ yihY → BrkB family putative transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,076,6040GA 5.7% 390.8 / 6.7 122L153L (CTG→CTA) ‡yihYBrkB family putative transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (61/53);  new base A (2/5);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 134 total aligned reads displayed.

GTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAT  >  NC_000913/4076455‑4076753
                                                                                                                                                     |                                                                                                                                                     
gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGGCGctgctg                                                                                                                                                       <  1:526812/150‑1 (MQ=255)
gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGctgctg                                                                                                                                                       <  2:1210080/150‑1 (MQ=255)
gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGctgctg                                                                                                                                                       <  1:1566172/150‑1 (MQ=255)
 ttGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTgg                                                                                                                                                      <  1:325559/150‑1 (MQ=255)
  tGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGc                                                                                                                                                     <  2:627344/150‑1 (MQ=255)
    gggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAg                                                                                                                                                   <  2:748111/150‑1 (MQ=255)
    gggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAg                                                                                                                                                   >  2:486515/1‑150 (MQ=255)
     ggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAgg                                                                                                                                                  <  2:241502/150‑1 (MQ=255)
        gTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGTCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGc                                                                                                                                               <  2:842408/150‑1 (MQ=255)
         tGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGcc                                                                                                                                              <  1:197228/150‑1 (MQ=255)
          gCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccc                                                                                                                                             <  2:1036309/150‑2 (MQ=255)
          gCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACAATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccc                                                                                                                                             <  2:1101628/150‑2 (MQ=255)
            gggCTGAGCGTCACCGCGGTCTTGTTGATCTACTCCATCGCTAGCGCGTGGAATACCATCTGGCGCAGTACACGAGCGCGACCAAAAATTTCCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGt                                                                                                                                           <  1:131104/150‑1 (MQ=255)
             ggCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTc                                                                                                                                          <  2:966076/150‑1 (MQ=255)
              gCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCt                                                                                                                                         <  2:1581513/150‑1 (MQ=255)
              gCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccccagc                                                                                                                                         <  1:1149957/150‑6 (MQ=255)
               cTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGACGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTg                                                                                                                                        >  1:237495/1‑150 (MQ=255)
                tGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTgg                                                                                                                                       >  1:1028594/1‑150 (MQ=255)
                tGATCGTCACCGCGGTATTGTTGATGTACTCCATCGATATCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTgg                                                                                                                                       <  1:340399/150‑1 (MQ=255)
                 gATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGAGTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGc                                                                                                                                      >  2:463017/1‑150 (MQ=255)
                     gTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGagc                                                                                                                                  >  1:163665/1‑148 (MQ=255)
                     gTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATc                                                                                                                                  >  1:1357573/1‑150 (MQ=255)
                        aCCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGt                                                                                                                               <  1:1661798/150‑1 (MQ=255)
                         ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                              >  2:1263684/1‑150 (MQ=255)
                         ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                              >  1:493932/1‑150 (MQ=255)
                         ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                              >  2:1696047/1‑150 (MQ=255)
                           gcgTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTcc                                                                                                                            <  1:1263684/150‑1 (MQ=255)
                            cgTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTAAACGCTGGGGCCGCTGCTGGCAGGGACCAGTCTGGCGATCAGTTCCt                                                                                                                           >  1:886822/1‑150 (MQ=255)
                             gTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTa                                                                                                                          <  1:1169007/150‑1 (MQ=255)
                              ttATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGAGTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAt                                                                                                                         >  2:1676341/1‑150 (MQ=255)
                                aTTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAttt                                                                                                                       <  1:1198091/150‑1 (MQ=255)
                                 ttgttgATGTACTCCATCGATAGCGCGTTGACTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTg                                                                                                                      <  1:1106235/150‑1 (MQ=255)
                                 ttgttgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTg                                                                                                                      >  1:1234627/1‑150 (MQ=255)
                                  tgttgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGc                                                                                                                     >  1:92689/1‑150 (MQ=255)
                                    ttgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGctc                                                                                                                   <  2:943916/150‑1 (MQ=255)
                                              ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAATTTTACTCGTTCGCCGTGTAATGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                         >  1:273407/1‑150 (MQ=255)
                                              ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                         >  2:495685/1‑150 (MQ=255)
                                              ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                         >  2:438423/1‑150 (MQ=255)
                                                aTCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTgg                                                                                                       >  1:1169032/1‑150 (MQ=255)
                                                 tCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGCGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                      >  1:944311/1‑150 (MQ=255)
                                                 tCGATAGCGCGTTGAATACCAGCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                      <  1:163341/150‑1 (MQ=255)
                                                  cGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgc                                                                                                     <  1:1195381/150‑1 (MQ=255)
                                                  cGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgc                                                                                                     <  2:492804/150‑1 (MQ=255)
                                                   gATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcg                                                                                                    <  1:1683517/150‑1 (MQ=255)
                                                     tAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                  <  2:132565/150‑1 (MQ=255)
                                                        cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgagcga                                                                                               >  1:793392/1‑150 (MQ=255)
                                                         gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                              <  2:1708313/150‑1 (MQ=255)
                                                         gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                              >  1:1038250/1‑150 (MQ=255)
                                                           gTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtct                                                                                            >  2:927973/1‑150 (MQ=255)
                                                             tGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCa                                                                                          >  2:1430361/1‑150 (MQ=255)
                                                               aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGCGGCCAGTCTGGCGATCAGTTCCTAATTGCTCTCCCTGCGCTGGGCGAGCGATCTCaaa                                                                                        >  2:700900/1‑149 (MQ=255)
                                                                 tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATAc                                                                                      >  1:376846/1‑150 (MQ=255)
                                                                 tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATAc                                                                                      >  2:1591963/1‑150 (MQ=255)
                                                                   ccATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTg                                                                                    >  1:1196574/1‑150 (MQ=255)
                                                                   ccATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCAGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTg                                                                                    >  1:1346041/1‑150 (MQ=255)
                                                                   ccATATGGCGAAGTAAACGAGCGCTACCAAAAATTTACTCGTTCCCCGTGTACTGGATGATTATAACGCTGGCGCCGCTGCTGGAAGGGGCCAGTCTGTCGATCAGTTCCTATTTGCTCTCACATCGCTGGGCGAGCGATCTCaaatctg                                                                                    >  1:351664/1‑145 (MQ=255)
                                                                        tGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATc                                                                               <  1:1495088/150‑1 (MQ=255)
                                                                        tGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATc                                                                               >  2:314110/1‑150 (MQ=255)
                                                                          gcgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGa                                                                             >  1:1513140/1‑150 (MQ=255)
                                                                             cAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACaa                                                                          <  1:1620076/150‑1 (MQ=255)
                                                                              aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                         >  1:1016798/1‑150 (MQ=255)
                                                                                 aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtg                                                                      >  1:122336/1‑150 (MQ=255)
                                                                                   aCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgct                                                                    <  1:314110/150‑1 (MQ=255)
                                                                                     gAGCGCGACCAAAAATTTACGCGTTCGCCGTCTACTGGATGATTTAAACGCTGGGGCCGCTGCTGGCTGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                                                  <  2:714299/150‑1 (MQ=255)
                                                                                      aGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCg                                                                 <  1:47885/150‑1 (MQ=255)
                                                                                        cgcgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTCCGTa                                                               >  2:1020869/1‑150 (MQ=255)
                                                                                          cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAtt                                                             >  1:798475/1‑150 (MQ=255)
                                                                                            aCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAtttt                                                           >  2:70753/1‑150 (MQ=255)
                                                                                                    tttACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGc                                                   <  1:438423/150‑1 (MQ=255)
                                                                                                       aCTCGTTCGCCGTGTACTGGATGATTTAAACGCTGGGGCCGCTGCTGGCAGGGACCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgt                                                <  2:886822/150‑1 (MQ=255)
                                                                                                          cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                             <  2:1028594/150‑1 (MQ=255)
                                                                                                          cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                             >  1:486605/1‑150 (MQ=255)
                                                                                                          cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTCGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                             <  2:686190/150‑1 (MQ=255)
                                                                                                           gTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTc                                            >  2:949130/1‑150 (MQ=255)
                                                                                                            ttCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCg                                           <  2:122336/150‑1 (MQ=255)
                                                                                                                ccGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                       <  2:8201/150‑1 (MQ=255)
                                                                                                                ccGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                       <  1:431773/150‑1 (MQ=255)
                                                                                                                ccGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGAGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                       >  2:958367/1‑150 (MQ=255)
                                                                                                                ccGGGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                       <  2:333269/150‑1 (MQ=255)
                                                                                                                  gtgtACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtc                                     <  2:776119/150‑1 (MQ=255)
                                                                                                                     tACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTcc                                  >  1:680720/1‑150 (MQ=255)
                                                                                                                         ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                                   >  2:117803/1‑113 (MQ=255)
                                                                                                                         ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCt                              >  1:812471/1‑150 (MQ=255)
                                                                                                                                  tAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCGGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGATTCTGTTGTCGTGGATCTCCTTCTGGTTGCTGt                     <  1:90275/150‑1 (MQ=255)
                                                                                                                                   aaCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTa                    >  1:1120548/1‑150 (MQ=255)
                                                                                                                                    aCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTAc                   >  1:249304/1‑150 (MQ=255)
                                                                                                                                    aCGCTGGGGCCGCAGCTGGAAGGGGCCAGTCTGGCGATCAGTAACTATTAGCTCTCCCTGCGCTCGGCGAGCGATCTCAACACTGACATCGACAACTTGCTGCGTATATTTACGTTGCTGTTGTCGAGGAACTCCTTCTGGTTGATGTAc                   >  2:1282278/1‑150 (MQ=255)
                                                                                                                                     cagactggCCCCTGCCAGCAGCGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACa                  >  1:257750/7‑150 (MQ=255)
                                                                                                                                      gCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAg                 >  2:1705457/1‑150 (MQ=255)
                                                                                                                                      gCTGGGGCCGCTGCTGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGc                >  2:1394311/1‑150 (MQ=255)
                                                                                                                                        tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTACTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTGGTGGATCTCCTTCTGGTTGCTGTACAGct               >  1:721985/1‑149 (MQ=255)
                                                                                                                                         ctggCCCCTGCCAGCAGCGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCAt              >  1:202954/3‑150 (MQ=255)
                                                                                                                                           ggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTg            >  2:792742/1‑150 (MQ=255)
                                                                                                                                            gCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTAAATACTGTCATCTACAACGTGCTGCGTATTTTTCCGTTTCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGt           <  1:203594/150‑1 (MQ=255)
                                                                                                                                                  gctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAc     <  2:376846/150‑1 (MQ=255)
                                                                                                                                                  gctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAc     <  2:1272889/150‑1 (MQ=255)
                                                                                                                                                  gctgGCAGGGGCCAGTCTGGCGATAAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCCATACTGTCATCGACAACGTGCGGCGGAGTTTTCCTTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTCTTCCTAc     >  1:1113269/1‑150 (MQ=255)
                                                                                                                                                    tgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCa   <  1:1696047/150‑1 (MQ=255)
                                                                                                                                                     ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAt  >  1:70828/1‑150 (MQ=255)
                                                                                                                                                     ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCCGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAt  >  2:535124/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAT  >  NC_000913/4076455‑4076753

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