Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,076,598 G→C 5.5% P151P (CCG→CCC 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,5980GC 5.5% 438.1 / 11.9 127P151P (CCG→CCCyihYBrkB family putative transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (62/58);  new base C (2/5);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 129 total aligned reads displayed.

GCCGTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGT  >  NC_000913/4076452‑4076744
                                                                                                                                                  |                                                                                                                                                  
gccgTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAACCGAGCGCGCCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTACCGCTGGGGCCgctg                                                                                                                                                 >  2:82785/1‑150 (MQ=255)
 ccgTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCgctgc                                                                                                                                                <  2:1554392/150‑1 (MQ=255)
   gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGGCGctgctg                                                                                                                                              <  1:526812/150‑1 (MQ=255)
   gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGctgctg                                                                                                                                              <  2:1210080/150‑1 (MQ=255)
   gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGctgctg                                                                                                                                              <  1:1566172/150‑1 (MQ=255)
   gTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGctgctg                                                                                                                                              <  1:382794/150‑1 (MQ=255)
       gggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAg                                                                                                                                          <  2:748111/150‑1 (MQ=255)
        ggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAgg                                                                                                                                         <  2:241502/150‑1 (MQ=255)
           gTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGTCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGc                                                                                                                                      <  2:842408/150‑1 (MQ=255)
            tGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGcc                                                                                                                                     <  2:517596/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)
                 gCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccccagc                                                                                                                                <  1:1149957/150‑6 (MQ=255)
                  cTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGGCTg                                                                                                                               >  1:1203369/1‑150 (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)
                     aTCGTCACCGCGTTATTGTCGATGAACTCCATCGATCGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCg                                                                                                                            <  1:920720/150‑1 (MQ=255)
                       cGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGAt                                                                                                                          >  1:1708313/1‑150 (MQ=255)
                        gTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGagc                                                                                                                         >  1:163665/1‑148 (MQ=255)
                        gTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATc                                                                                                                         >  1:1357573/1‑150 (MQ=255)
                         tCACCGCGTTATTGTTGCTGTACTCCATCGATAGCGCGTTGAATAGCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTGCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCa                                                                                                                        <  1:1091437/150‑1 (MQ=255)
                            ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                     >  2:1696047/1‑150 (MQ=255)
                            ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                     >  1:493932/1‑150 (MQ=255)
                               cgTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTAAACGCTGGGGCCGCTGCTGGCAGGGACCAGTCTGGCGATCAGTTCCt                                                                                                                  >  1:886822/1‑150 (MQ=255)
                                 ttATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGAGTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAt                                                                                                                >  2:1676341/1‑150 (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:438423/1‑150 (MQ=255)
                                                 ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                >  2:495685/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)
                                                          gcgcGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgagcg                                                                                       <  2:33999/150‑1 (MQ=255)
                                                           cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgagcga                                                                                      >  1:793392/1‑150 (MQ=255)
                                                            gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                     >  1:1038250/1‑150 (MQ=255)
                                                            gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                     <  2:1708313/150‑1 (MQ=255)
                                                              gTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtct                                                                                   >  2:927973/1‑150 (MQ=255)
                                                                tGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCa                                                                                 >  2:1430361/1‑150 (MQ=255)
                                                                  aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGCGGCCAGTCTGGCGATCAGTTCCTAATTGCTCTCCCTGCGCTGGGCGAGCGATCTCaaa                                                                               >  2:700900/1‑149 (MQ=255)
                                                                    tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATAc                                                                             >  2:1591963/1‑150 (MQ=255)
                                                                      ccATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTg                                                                           >  1:1196574/1‑150 (MQ=255)
                                                                      ccATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGAACTGGATGATTTTAACGCTGCTGCCGCTGCTGGCAGGGGCCAGACAGGCGATCAGTTCCTATTTTCTCTCCCTGCGCAGGGCGAGCGATCTCAATACtc                                                                           >  1:214434/1‑149 (MQ=255)
                                                                      ccATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCAGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTg                                                                           >  1:1346041/1‑150 (MQ=255)
                                                                      ccATATGGCGAAGTAAACGAGCGCTACCAAAAATTTACTCGTTCCCCGTGTACTGGATGATTATAACGCTGGCGCCGCTGCTGGAAGGGGCCAGTCTGTCGATCAGTTCCTATTTGCTCTCACATCGCTGGGCGAGCGATCTCaaatctg                                                                           >  1:351664/1‑145 (MQ=255)
                                                                           tGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATc                                                                      <  1:1495088/150‑1 (MQ=255)
                                                                             gcgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGa                                                                    >  1:1513140/1‑150 (MQ=255)
                                                                                cAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACaa                                                                 <  1:1620076/150‑1 (MQ=255)
                                                                                 aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                >  1:1016798/1‑150 (MQ=255)
                                                                                 aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                <  2:949447/150‑1 (MQ=255)
                                                                                 aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                >  1:810393/1‑150 (MQ=255)
                                                                                  gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                               <  1:1656218/150‑1 (MQ=255)
                                                                                  gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                               >  2:695097/1‑150 (MQ=255)
                                                                                    aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtg                                                             <  2:736800/150‑1 (MQ=255)
                                                                                    aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtg                                                             >  1:122336/1‑150 (MQ=255)
                                                                                      aCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgct                                                           <  1:314110/150‑1 (MQ=255)
                                                                                       cGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                          <  1:1133409/150‑1 (MQ=255)
                                                                                       cGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                          >  2:1692602/1‑150 (MQ=255)
                                                                                         aGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCg                                                        <  2:92689/150‑1 (MQ=255)
                                                                                         aGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCg                                                        <  1:47885/150‑1 (MQ=255)
                                                                                           cgcgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTCCGTa                                                      >  2:1020869/1‑150 (MQ=255)
                                                                                               aCCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAtttt                                                  >  2:70753/1‑150 (MQ=255)
                                                                                                       tttACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGc                                          <  1:438423/150‑1 (MQ=255)
                                                                                                          aCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgt                                       >  1:791863/1‑150 (MQ=255)
                                                                                                           cTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgtt                                      >  2:1577959/1‑150 (MQ=255)
                                                                                                             cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                    <  2:1028594/150‑1 (MQ=255)
                                                                                                             cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                    >  1:486605/1‑150 (MQ=255)
                                                                                                             cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTCGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                    <  2:686190/150‑1 (MQ=255)
                                                                                                              gTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTc                                   >  1:1272889/1‑150 (MQ=255)
                                                                                                               ttCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCg                                  <  2:122336/150‑1 (MQ=255)
                                                                                                                 cGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTg                                <  2:1234627/150‑1 (MQ=255)
                                                                                                                 cGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTg                                <  2:793392/150‑1 (MQ=255)
                                                                                                                   ccGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                              <  1:431773/150‑1 (MQ=255)
                                                                                                                   ccGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                              <  2:8201/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)
                                                                                                                         aCTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCt                        >  2:1322218/1‑150 (MQ=255)
                                                                                                                           tGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCAAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTc                      >  1:337348/1‑150 (MQ=255)
                                                                                                                            ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                          <  1:117803/113‑1 (MQ=255)
                                                                                                                            ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                          >  2:117803/1‑113 (MQ=255)
                                                                                                                            ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCt                     >  1:812471/1‑150 (MQ=255)
                                                                                                                             gaggatTTTAACGCTGGGGCCGATGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTg                    <  2:237495/147‑1 (MQ=255)
                                                                                                                                      aaCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTa           >  1:1120548/1‑150 (MQ=255)
                                                                                                                                       aCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTAc          >  1:249304/1‑150 (MQ=255)
                                                                                                                                       aCGCTGGGGCCGCAGCTGGAAGGGGCCAGTCTGGCGATCAGTAACTATTAGCTCTCCCTGCGCTCGGCGAGCGATCTCAACACTGACATCGACAACTTGCTGCGTATATTTACGTTGCTGTTGTCGAGGAACTCCTTCTGGTTGATGTAc          >  2:1282278/1‑150 (MQ=255)
                                                                                                                                        cGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCGCTGCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGATGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACa         <  1:263903/150‑1 (MQ=255)
                                                                                                                                           tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTACTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTGGTGGATCTCCTTCTGGTTGCTGTACAGct      >  1:721985/1‑149 (MQ=255)
                                                                                                                                               gCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTAAATACTGTCATCTACAACGTGCTGCGTATTTTTCCGTTTCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGt  <  1:203594/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGT  >  NC_000913/4076452‑4076744

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: