Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,172,941 C→A E342* (GAA→TAA)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,172,9410CA96.5% 174.8 / NA 57E342* (GAA→TAA) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (29/26);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCCACGGTCGCTGGCTTCATCAAGACGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTCAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGTGGTTCGCTGAAGAACAAT  >  NC_000913/1172695‑1173179
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
aaCGCCCACGGTCGCTGGCGTCATCAAGACGCATAATGCGGTGCGGAGCAATTTTAATTCGTTCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACGTACGCAGCGCATCCAGCGGCGCCGTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGtt                                                                                                                                                                                                                                            <  1:584745/251‑1 (MQ=255)
                ggCTTCATCAAGACGAATTATGCGTTGCGGAGCAATTTTAAGTCTTGCGAGCAGTGCACCCAGCGCTTCACGGCGACCTTCTCTTCCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTATGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCATTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCCCCCggg                                                                                                                                                                                                                            <  2:376719/251‑1 (MQ=255)
                    tcatcaAGACGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGcc                                                                                                                                                                                                                         >  2:362872/1‑250 (MQ=255)
                     catcaAGACGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCa                                                                                                                                                                                                                        <  1:193597/250‑1 (MQ=255)
                      atcaAGACGCATAATCCGTTGCGGAGCAATTTTAATTCGTACGAGCAGTTCACCCAGCGCTTCACGGCGACATTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGCGGTACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAACTTGGACTTCGCGGCTTTTGTCGGGAATTGGTAAGGTTTTAGCTGCACCCGGGGccggt                                                                                                                                                                                                                      >  2:93285/1‑248 (MQ=255)
                        caAGACGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGGCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGttt                                                                                                                                                                                                                    <  2:206114/251‑1 (MQ=255)
                             cgcctaacgCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACATTCCCTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCCGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCAGTCGCGGCTTTTGTCGGTAAATGTTCAGTTTTTAGCTGCCCCCGGGGCCAGTGTTTc                                                                                                                                                                                                                 >  1:523812/9‑249 (MQ=255)
                             cGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTc                                                                                                                                                                                                                 >  1:523829/1‑249 (MQ=255)
                                aTAATGCGTTGCGGAGCCATTAAAATTCGTGCAAGCCGATCTCCCAGAGCTTCAAGGCGACTTACACTTTCTACCGCGAACACCACCGGACCGTCGAAAGTCTCGAGGACCTTCCGCAGCGCATCCAGCGGCGCTTATTGTTGTGCCTGCACGGCCAGGACTGGCAGTTTCTGGAAACCTCAACTGGCGTTCGCGGCTTTGGTCGGTGAATGTTAAAGTTTGAACTGCACCCGGGGCCAGGTTTCCAGccc                                                                                                                                                                                                            >  2:62525/1‑249 (MQ=255)
                                    tGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTg                                                                                                                                                                                                          >  1:556052/1‑249 (MQ=255)
                                            ggAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTc                                                                                                                                                                                                >  1:670477/1‑251 (MQ=255)
                                                cAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGCAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTACCTGTTACGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGGGCTCGTcc                                                                                                                                                                                            >  2:607719/1‑251 (MQ=255)
                                                 aaTTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGACAGTCTCGAGGAACTTACGCAGCGCATCCAGAGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTcc                                                                                                                                                                                            >  1:697752/1‑250 (MQ=255)
                                                 aaTTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTccc                                                                                                                                                                                           >  1:389163/1‑250 (MQ=255)
                                                       aaTTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCCCCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCa                                                                                                                                                                                     <  2:136566/251‑1 (MQ=255)
                                                        aTTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTGCAGCTCTGAGAAGAGCTCGGCCACCCGCAg                                                                                                                                                                                    <  2:664185/251‑1 (MQ=255)
                                                                 aGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCCGGAACCCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGAACCAGGGGCCAGTTTTTCAGCTCTGAGAAGCGCTCGTCCACCCGCACCCAGcgcgg                                                                                                                                                                           >  2:216701/1‑250 (MQ=255)
                                                                   cAGTTCACACAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGGCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCGGAGAAGAGCTCGGCCACCCGCAGCCAGAGCGAtt                                                                                                                                                                         <  2:192290/251‑1 (MQ=255)
                                                                      ttCACCCAGCGCTTCACGGCGACCTTCCCTTTCTACCGAGAACACCACCGGACCGTCGAACGTCTCGAGGAACTTACGCAGCGCATCCCGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGGCGGTCAATGTTAAGCTTTTAGCTGCCCCCGGGGCCAGTTTTTCAGCCCTGAGCAGAGCTCGTCCACCCGCAGCCAGAGCGATtg                                                                                                                                                                        >  2:366785/1‑249 (MQ=255)
                                                                          cccAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGtggtgg                                                                                                                                                                  >  1:452883/1‑251 (MQ=255)
                                                                                      cGGCGACTTTCACTTTCTCCCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCACCCAGAGCGATTGTGGTGGCAACAGCGGgcg                                                                                                                                                      >  2:240816/1‑251 (MQ=255)
                                                                                       ggCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTTTGGTGGCAACAGCGGgcgc                                                                                                                                                     >  2:182104/1‑251 (MQ=255)
                                                                                         cGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCAt                                                                                                                                                   >  1:84078/1‑251 (MQ=255)
                                                                                          gACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCAACCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATc                                                                                                                                                  >  2:503390/1‑251 (MQ=255)
                                                                                               tCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGa                                                                                                                                               <  1:658394/249‑1 (MQ=255)
                                                                                                          cGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGAt                                                                                                                                  <  1:368253/251‑1 (MQ=255)
                                                                                                           gagaACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGAtt                                                                                                                                 >  2:151085/1‑251 (MQ=255)
                                                                                                                      cGGACCGTCCAAACTCCAGAGGAAAGTACGCAGCGCCTCCAGCGACGCTTTTTGTTGTGCCTGAACGGCCAGCTCTGGCAGCTTCTGGCAACCTAACTTCGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAACAGCTCGTCCACCCCCCGCCAGAGCGATTGTGGTGGCAACCGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAAcgcg                                                                                                                      >  2:420037/1‑251 (MQ=255)
                                                                                                                       ggACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGCAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTCAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGCTTCTCaaaccgagc                                                                                                                     >  2:368295/1‑245 (MQ=255)
                                                                                                                          ccGTCGAAAGGCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGGAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCGCAAAACGCGCCAg                                                                                                                  <  2:576686/251‑1 (MQ=255)
                                                                                                                             tCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCagag                                                                                                                                                                              <  1:381251/188‑1 (MQ=255)
                                                                                                                             tCGAAAGTCTAGAGGCACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGATTCTGGAACCCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTTAGAAGAGCTCGTCCACCCGCAGCCagag                                                                                                                                                                              >  2:381251/1‑188 (MQ=255)
                                                                                                                              cGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCTATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGGCAGCgtg                                                                                                              <  2:623480/251‑1 (MQ=255)
                                                                                                                                  aGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTc                                                                                                            >  1:215877/1‑249 (MQ=255)
                                                                                                                                                tACGCAGGGAATCCAGCGGCGATTTCTGTAGAGCCTGAAAGGCCAGGTCTGGCAGCTTCTGGAAAACTAAATTGGCATTAGCGGCTTTTGTCGGTCAATGTTAAGGTTTTAGATGCATCCAGGGCCAGTTTTTCAGCTCTGAGAGGAGCTCGTCCACACGCAGCCCGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCACGATTCTCAAAACGCGTCAGCTTGTCAGACTGGaaa                                                                                                  >  2:690683/1‑245 (MQ=255)
                                                                                                                                                      gcgcATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGGTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGTTTGTGGTGGCAACAGCGGGCGCATCTGATCGACGCCGCGATGCTCAAAACGCGCCAGCGTTTCAGCCTGGAAACGTTCGGCACTg                                                                                      <  2:371985/251‑1 (MQ=255)
                                                                                                                                                       cgcATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCACAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTgg                                                                                     >  1:461803/1‑251 (MQ=255)
                                                                                                                                                                       ttttGTTGTTCCTGAACGGCCAGGTCTGGTAGTTTCTGGAAACCTAAAGTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCGGCACCCGGGGCCAGTGTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCa                                                                       <  2:270261/249‑1 (MQ=255)
                                                                                                                                                                                    gAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCa                                                          >  1:197921/1‑249 (MQ=255)
                                                                                                                                                                                         gCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGGATTCACGAGCAAGGTAtt                                                   <  1:182104/251‑1 (MQ=255)
                                                                                                                                                                                             ggTCTGGCAGTTTCTGGAACCCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGc                                                >  1:230199/1‑250 (MQ=255)
                                                                                                                                                                                              gTCTGGCAGTTTCTGGTAACCTAAATTGGCATTCGCGGCTTTTGGCGGTAAATGTTAAGTTTTTAGCTGCGCCCGGGGACAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAg                                              <  2:523812/251‑1 (MQ=255)
                                                                                                                                                                                              gTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAg                                              <  2:523829/251‑1 (MQ=255)
                                                                                                                                                                                                  ggCAGTTTCTTGAAACCTAAATTGTCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTGTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCTATTGTGGTGGCAACAGCGGGCGCATCGGATCGATGCCGCGATTCTCAAAACGCGCTAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAgg                                             <  1:294693/248‑1 (MQ=255)
                                                                                                                                                                                                     aGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGGGCACGGAATTGGCCGGGAAAt                                        >  2:604071/1‑250 (MQ=255)
                                                                                                                                                                                                     aGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATa                                       <  1:93285/251‑1 (MQ=255)
                                                                                                                                                                                                            gtaaaCCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGTTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTTGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGTGAAATAACGGAAc                                <  2:129242/249‑1 (MQ=255)
                                                                                                                                                                                                            ggAAACCTAAAGTGGCAGTCGCGGCTTTTGGCGGTAAATGTGAAGGTGTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGaa                                 <  2:556052/250‑1 (MQ=255)
                                                                                                                                                                                                                aCCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAgcg                            <  1:362872/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccTAAATTGGCATTCGCAGCTTTCGTCGGCAAATGTTAAGTTTTTAGCTGCACACGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACACGCAGACAGAGCGCTTGTGGTGGCAAAAGCGGGCGCATCGGATCGACGACGCGCTTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCCAGAGCAAGGAATTGGCAGGGAAATAACTGAACAgcgg                           >  2:466476/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttGGCATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCggcggc                       >  1:4238/1‑249 (MQ=255)
                                                                                                                                                                                                                           cATTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTAAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAAAGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGGGGt                 >  1:584859/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCCGTGGtt                >  2:580827/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTTCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGTGGtt                <  1:505917/251‑1 (MQ=255)
                                                                                                                                                                                                                              tCGCGGCTTTTGTCGGTAAATGTTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGTGTTTc               <  2:670477/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tggTCGGTAAATGTTAAGTTTTTAGCGGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGGCCACCCGCAGCCAGAGCGATTGTGGGGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAAGAACGGAATAGCGGCGGCAGTGGTTCGCTgaagaa      <  2:461803/248‑1 (MQ=255)
                                                                                                                                                                                                                                           cGGTAAATGGTAAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCGCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGGATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGTGGTTCGCTGAAGAACAAt  <  1:607719/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
AACGCCCACGGTCGCTGGCTTCATCAAGACGCATAATGCGTTGCGGAGCAATTTTAATTCGTGCGAGCAGTTCACCCAGCGCTTCACGGCGACCTTCACTTTCTACCGAGAACACCACCGGACCGTCGAAAGTCTCGAGGAACTTACGCAGCGCATCCAGCGGCGCTTTTTGTTGTGCCTGAACGGCCAGGTCTGGCAGTTTCTGGAAACCTAAATTGGCATTCGCGGCTTTTGTCGGTAAATGTTCAGTTTTTAGCTGCACCCGGGGCCAGTTTTTCAGCTCTGAGAAGAGCTCGTCCACCCGCAGCCAGAGCGATTGTGGTGGCAACAGCGGGCGCATCGGATCGACGCCGCGATTCTCAAAACGCGCCAGCGTGTCAGCCTGGAAACGTTCGGCACTGGTTTCCAGATCGCCAGTATTCACGAGCAAGGTATTGGCAGGGAAATAACTGAACAGCGGCGGCAGTGGTTCGCTGAAGAACAAT  >  NC_000913/1172695‑1173179

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