Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,411,903 T→C 5.2% N16D (AAC→GAC)  ttcA ← tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,411,9030TC 5.2% 240.3 / 2.6 77N16D (AAC→GAC) ttcAtRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCAT  >  NC_000913/1411761‑1412049
                                                                                                                                              |                                                                                                                                                  
tgcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGtt                                                                                                                                              >  3:185189/1‑149 (MQ=255)
tgcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGtt                                                                                                                                              >  1:353121/1‑149 (MQ=255)
 gcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTc                                                                                                                                                   <  8:145711/143‑1 (MQ=255)
 gcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTc                                                                                                                                                   >  7:145711/1‑143 (MQ=255)
 gcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                                                                                                                                             >  8:336352/1‑149 (MQ=255)
 gcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTGTTTGTAATTTGTTCAGGTTg                                                                                                                                             >  8:381160/1‑149 (MQ=255)
  cAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGtt                                                                                                                                              <  8:194871/147‑1 (MQ=255)
  cAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                                                                                                                                             >  2:24790/1‑148 (MQ=255)
   aaaTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTa                                                                                                                                           <  6:317580/149‑1 (MQ=255)
    aaTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTAt                                                                                                                                          <  6:223057/149‑1 (MQ=255)
      ttGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTg                                                                                                                                        <  4:316055/149‑1 (MQ=255)
         cgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttc                                                                                                                                     <  5:290268/149‑1 (MQ=255)
         cgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGtt                                                                                                                                      <  1:24790/148‑1 (MQ=255)
         cgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGtt                                                                                                                                      <  6:6641/148‑1 (MQ=255)
          gcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttct                                                                                                                                    >  2:127127/1‑149 (MQ=255)
              aaTCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttctttct                                                                                                                                <  2:11343/149‑1 (MQ=255)
                   ccAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGt                                                                                                                            >  1:250405/1‑148 (MQ=255)
                     aGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAAt                                                                                                                         <  2:250405/149‑1 (MQ=255)
                         tGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttgt                                                                                                                     <  8:69124/149‑1 (MQ=255)
                           gTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttgttg                                                                                                                   <  1:127127/149‑1 (MQ=255)
                             ataGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGAt                                                                                                                 >  5:319744/1‑149 (MQ=255)
                             ataGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGAt                                                                                                                 >  7:52910/1‑149 (MQ=255)
                              taGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACATTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGAtt                                                                                                                >  3:252475/1‑149 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                                                                                                                                                >  1:377699/1‑114 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                                                                                                                                                <  2:377699/114‑1 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTAtt                                                                                                                                         <  6:99927/121‑1 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTAtt                                                                                                                                         >  5:99927/1‑121 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTc                                                                                                             >  4:182638/1‑149 (MQ=255)
                                 cTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTc                                                                                                             <  3:201114/149‑1 (MQ=255)
                                   gTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCtt                                                                                                           >  3:356157/1‑149 (MQ=255)
                                     cTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTg                                                                                                          >  2:317614/1‑148 (MQ=255)
                                         aCCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCatta                                                                                                     >  5:302898/1‑149 (MQ=255)
                                           ccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCattatt                                                                                                   <  7:336352/149‑1 (MQ=255)
                                             ccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATCTTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTc                                                                                                 <  2:422314/149‑1 (MQ=255)
                                              cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTTTTTGTAATTTGTTGATTTTCTTGCATTATTTCa                                                                                                <  5:26544/149‑1 (MQ=255)
                                              cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTcc                                                                                                <  7:297460/149‑2 (MQ=255)
                                                gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGGTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCAGTATTTCAGt                                                                                              >  8:328277/1‑149 (MQ=255)
                                                     gCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTctct                                                                                         >  4:387882/1‑149 (MQ=255)
                                                       aaaCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                                       <  7:359855/149‑1 (MQ=255)
                                                       aaaCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTgg                                                                                       <  3:387882/149‑1 (MQ=255)
                                                        aaCCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGt                                                                                      <  7:381160/149‑1 (MQ=255)
                                                                    aTCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATgggg                                                                           >  1:129971/1‑148 (MQ=255)
                                                                         cttcttCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAAtt                                                                     >  4:326389/1‑149 (MQ=255)
                                                                             gtcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGCCGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTTAATGGGGCAAATTgggg                                                                 <  5:334548/148‑1 (MQ=255)
                                                                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTggggg                                                                >  2:9982/1‑148 (MQ=255)
                                                                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTggg                                                                  <  4:356157/146‑1 (MQ=255)
                                                                                aaTCATATTGAAGTCAGCAAAGGCTTCGCCCACGTTACGACGCAGACGGTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGCTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGc                                                               >  4:314893/1‑148 (MQ=255)
                                                                                     tatTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACtt                                                         >  3:365445/1‑149 (MQ=255)
                                                                                      atTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCatta                                                                                                     <  3:113832/104‑1 (MQ=255)
                                                                                      atTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCatta                                                                                                     >  4:113832/1‑104 (MQ=255)
                                                                                      atTGAAGTCAGCAATGGCTTCGCCAACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACttt                                                        >  5:370977/1‑149 (MQ=255)
                                                                                       ttGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACttt                                                        <  3:182638/148‑1 (MQ=255)
                                                                                         gAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGc                                                      <  8:52910/148‑1 (MQ=255)
                                                                                         gAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGc                                                      <  6:319744/148‑1 (MQ=255)
                                                                                          aaGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGTTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCaa                                                    >  6:333426/1‑149 (MQ=255)
                                                                                             gcAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTa                                                 <  4:356946/148‑1 (MQ=255)
                                                                                                 cAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTggg                                                                  >  1:153443/1‑128 (MQ=255)
                                                                                                 cAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTggg                                                                  <  2:153443/128‑1 (MQ=255)
                                                                                                  aaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGAt                                             >  5:181292/1‑148 (MQ=255)
                                                                                                    tGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAAc                                          >  1:264823/1‑149 (MQ=255)
                                                                                                       cTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACcgc                                       <  2:353121/149‑1 (MQ=255)
                                                                                                        ttCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACcgcg                                      >  1:132353/1‑149 (MQ=255)
                                                                                                            cccACGTTACGACGCAGACGTTTTTGTAATTTGTCCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCAtt                                  >  1:268845/1‑149 (MQ=255)
                                                                                                            cccACGTTACGACGCAGACGTTTTTGTAATTTGTCCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCAt                                   >  1:172886/1‑148 (MQ=255)
                                                                                                               aCGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCaa                               >  6:354923/1‑149 (MQ=255)
                                                                                                               aCGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCaa                               >  5:91036/1‑149 (MQ=255)
                                                                                                                cGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAAc                              <  3:326389/149‑1 (MQ=255)
                                                                                                                 gTTACTACGCAGACGTTTTTGTAATTTGTCCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACa                             <  2:268845/149‑1 (MQ=255)
                                                                                                                 gTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACa                             >  8:229430/1‑149 (MQ=255)
                                                                                                                 gTTACGACGCAGACGTTTTTGTAATTTGTCCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACa                             <  2:172886/149‑1 (MQ=255)
                                                                                                                   tACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATg                           >  3:351140/1‑149 (MQ=255)
                                                                                                                     cgacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATgg                          <  4:365445/148‑1 (MQ=255)
                                                                                                                      gacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGc                         >  1:228985/1‑148 (MQ=255)
                                                                                                                      gacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCt                        <  7:170614/149‑1 (MQ=255)
                                                                                                                         gCAGTCGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATc                     <  1:162657/149‑1 (MQ=255)
                                                                                                                                         tttGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGTTGCAAATTGGGGGCATACTTTGCAACTACGATAACCGCGCATTCACCATGGCTATCTGTTCGTCGTTCATGt     >  3:257682/1‑149 (MQ=255)
                                                                                                                                            gTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca   <  2:132353/148‑1 (MQ=255)
                                                                                                                                             ttCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcat  <  2:264823/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCAT  >  NC_000913/1411761‑1412049

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