Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,076,609 G→C 6.1% G155A (GGG→GCG)  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,6090GC 6.1% 422.0 / 12.9 127G155A (GGG→GCG) yihYBrkB family putative transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (65/54);  new base C (3/5);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 131 total aligned reads displayed.

GGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCCGCG  >  NC_000913/4076460‑4076758
                                                                                                                                                     |                                                                                                                                                     
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)
        ggCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTc                                                                                                                                               <  2:966076/150‑1 (MQ=255)
         gCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCt                                                                                                                                              <  2:1581513/150‑1 (MQ=255)
          cTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGACGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTg                                                                                                                                             >  1:237495/1‑150 (MQ=255)
           tGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTgg                                                                                                                                            >  1:1028594/1‑150 (MQ=255)
           tGATCGTCACCGCGGTATTGTTGATGTACTCCATCGATATCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTgg                                                                                                                                            <  1:340399/150‑1 (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)
                  cACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAg                                                                                                                                     >  2:1133409/1‑150 (MQ=255)
                    ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                                   >  2:1696047/1‑150 (MQ=255)
                    ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                                   >  2:1263684/1‑150 (MQ=255)
                    ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                                   >  1:493932/1‑150 (MQ=255)
                      gcgTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTcc                                                                                                                                 <  1:1263684/150‑1 (MQ=255)
                        gTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTa                                                                                                                               <  1:1169007/150‑1 (MQ=255)
                         ttATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGAGTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAt                                                                                                                              >  2:1676341/1‑150 (MQ=255)
                           aTTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAttt                                                                                                                            <  1:1198091/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:438423/1‑150 (MQ=255)
                                         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                              >  2:495685/1‑150 (MQ=255)
                                           aTCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTgg                                                                                                            >  1:1169032/1‑150 (MQ=255)
                                            tCGATAGCGCGTTGAATACCAGCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                           <  1:163341/150‑1 (MQ=255)
                                             cGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgc                                                                                                          <  2:492804/150‑1 (MQ=255)
                                             cGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgc                                                                                                          <  1:1195381/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)
                                                    gcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                                   >  1:1038250/1‑150 (MQ=255)
                                                      gTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtct                                                                                                 >  2:927973/1‑150 (MQ=255)
                                                          aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGCGGCCAGTCTGGCGATCAGTTCCTAATTGCTCTCCCTGCGCTGGGCGAGCGATCTCaaa                                                                                             >  2:700900/1‑149 (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)
                                                                  cTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCAt                                                                                     <  1:654392/150‑1 (MQ=255)
                                                                   tGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATc                                                                                    >  2:314110/1‑150 (MQ=255)
                                                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCg                                                                                   >  2:1026714/1‑150 (MQ=255)
                                                                     gcgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGa                                                                                  >  1:1513140/1‑150 (MQ=255)
                                                                         aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                              >  1:1016798/1‑150 (MQ=255)
                                                                         aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                              >  1:810393/1‑150 (MQ=255)
                                                                         aGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAAc                                                                              <  2:949447/150‑1 (MQ=255)
                                                                          gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                             >  2:695097/1‑150 (MQ=255)
                                                                          gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                             <  1:1656218/150‑1 (MQ=255)
                                                                          gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                             >  1:1434666/1‑150 (MQ=255)
                                                                            aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtg                                                                           <  2:736800/150‑1 (MQ=255)
                                                                            aaaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtg                                                                           >  1:122336/1‑150 (MQ=255)
                                                                               cGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                                        <  1:1133409/150‑1 (MQ=255)
                                                                               cGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                                        >  2:1692602/1‑150 (MQ=255)
                                                                                gAGCGCGACCAAAAATTTACGCGTTCGCCGTCTACTGGATGATTTAAACGCTGGGGCCGCTGCTGGCTGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                                                       <  2:714299/150‑1 (MQ=255)
                                                                                 aGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCg                                                                      <  2:92689/150‑1 (MQ=255)
                                                                                   cgcgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTCCGTa                                                                    >  2:1020869/1‑150 (MQ=255)
                                                                                     cgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAtt                                                                  >  1:798475/1‑150 (MQ=255)
                                                                                               tttACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGc                                                        <  1:438423/150‑1 (MQ=255)
                                                                                                  aCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgt                                                     <  1:1591963/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)
                                                                                                         cGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTg                                              <  2:1234627/150‑1 (MQ=255)
                                                                                                         cGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTg                                              <  2:793392/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)
                                                                                                             gtgtACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtc                                          <  2:776119/150‑1 (MQ=255)
                                                                                                                 aCTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCt                                      >  2:1322218/1‑150 (MQ=255)
                                                                                                                   tGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTc                                    <  1:486515/150‑1 (MQ=255)
                                                                                                                   tGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCAAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTc                                    >  1:337348/1‑150 (MQ=255)
                                                                                                                    ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                                        >  2:117803/1‑113 (MQ=255)
                                                                                                                    ggATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCt                                   >  1:812471/1‑150 (MQ=255)
                                                                                                                     gaggatTTTAACGCTGGGGCCGATGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTg                                  <  2:237495/147‑1 (MQ=255)
                                                                                                                             tAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCGGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGATTCTGTTGTCGTGGATCTCCTTCTGGTTGCTGt                          <  1:90275/150‑1 (MQ=255)
                                                                                                                              aaCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTa                         >  1:1120548/1‑150 (MQ=255)
                                                                                                                                cagactggCCCCTGCCAGCAGCGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACa                       >  1:257750/7‑150 (MQ=255)
                                                                                                                                cGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCGCTGCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGATGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACa                       <  1:263903/150‑1 (MQ=255)
                                                                                                                                 gCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAg                      >  2:1705457/1‑150 (MQ=255)
                                                                                                                                 gCTGGGGCCGCTGCTGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGc                     >  2:1394311/1‑150 (MQ=255)
                                                                                                                                   tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTACTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTGGTGGATCTCCTTCTGGTTGCTGTACAGct                    >  1:721985/1‑149 (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)
                                                                                                                                                ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCaa       >  2:465950/1‑149 (MQ=255)
                                                                                                                                                ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAt       >  1:70828/1‑150 (MQ=255)
                                                                                                                                                ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCCGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAt       >  2:535124/1‑150 (MQ=255)
                                                                                                                                                  cAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTTTCGTGGATCTCCTTCTGGTTGATGTACAGCATTGTTACTACCATcc     >  1:116756/1‑150 (MQ=255)
                                                                                                                                                   aGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCCTCTGGTTGCTGTACAGCATTGTTCCTACCATCcg    <  1:379807/150‑1 (MQ=255)
                                                                                                                                                     gggCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCcgcg  <  2:1539845/150‑1 (MQ=255)
                                                                                                                                                     gggCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGGTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCcgcg  <  2:1203369/150‑1 (MQ=255)
                                                                                                                                                     gggCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGATGTTGTCGTGGATCTACTTCTGGTTGCTGTACAGCATTGTTCCTACCATCcgcg  >  1:159910/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCATCCGCG  >  NC_000913/4076460‑4076758

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: