Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,257,783 C→T 5.7% L10L (TTG→TTA ychF ← catalase inhibitor protein; ATPase, K+‑dependent, ribosome‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,257,7830CT 5.7% 233.4 / 3.8 70L10L (TTG→TTAychFcatalase inhibitor protein; ATPase, K+‑dependent, ribosome‑associated
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base T (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TTACGATTTCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGACAACGACTT  >  NC_000913/1257638‑1257930
                                                                                                                                                 |                                                                                                                                                   
ttACGATTTCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCaaa                                                                                                                                                  >  3:233253/1‑149 (MQ=255)
    gATTTCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACcga                                                                                                                                              <  4:641087/149‑1 (MQ=255)
     aTTTCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACcgac                                                                                                                                             <  3:449855/149‑1 (MQ=255)
      tttCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACcgacg                                                                                                                                            >  4:60282/1‑149 (MQ=255)
        tCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGAt                                                                                                                                          >  4:121020/1‑149 (MQ=255)
               gTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCAt                                                                                                                                   >  2:220668/1‑149 (MQ=255)
                 tGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCAttt                                                                                                                                 <  1:377263/149‑1 (MQ=255)
                   aTCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGa                                                                                                                               >  1:615188/1‑149 (MQ=255)
                       aGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATcc                                                                                                                           >  1:425268/1‑149 (MQ=255)
                                    gCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATcc                                                                                                                           <  2:646142/136‑1 (MQ=255)
                                    gCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATcc                                                                                                                           >  1:646142/1‑136 (MQ=255)
                                               acgCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATc                                                                                                   <  4:664496/149‑1 (MQ=255)
                                                  ccTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTta                                                                                                <  1:595448/149‑1 (MQ=255)
                                                   cTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTtaa                                                                                               <  3:47628/149‑1 (MQ=255)
                                                      tgtTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTtaataa                                                                                            >  3:199052/1‑149 (MQ=255)
                                                           ggCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAAc                                                                                       <  4:227822/149‑1 (MQ=255)
                                                            gCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGAtt                                                                                                                    >  1:267805/1‑119 (MQ=255)
                                                            gCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGAtt                                                                                                                    <  2:267805/119‑1 (MQ=255)
                                                            gCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAAcc                                                                                      >  1:471433/1‑149 (MQ=255)
                                                             cTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCt                                                                                     <  3:527746/149‑1 (MQ=255)
                                                             cTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCt                                                                                     >  4:314112/1‑149 (MQ=255)
                                                              tCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTg                                                                                    >  1:86327/1‑149 (MQ=255)
                                                              tCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGACAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCACGATTTATATCACCTTAATATCTTAATAATCAACCTg                                                                                    >  4:467975/1‑149 (MQ=255)
                                                                      gCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCt                                                                            >  2:500456/1‑149 (MQ=255)
                                                                        aGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACGGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTaa                                                                          >  3:140898/1‑149 (MQ=255)
                                                                        aGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTaa                                                                          >  3:344235/1‑149 (MQ=255)
                                                                         gAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAAc                                                                         >  3:362105/1‑149 (MQ=255)
                                                                         gAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAAc                                                                         <  4:233253/149‑1 (MQ=255)
                                                                         gAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAAc                                                                         <  4:199052/149‑1 (MQ=255)
                                                                          aaTGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACa                                                                        <  2:344976/149‑1 (MQ=255)
                                                                             ggAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGTAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCt                                                                                                  >  1:127399/1‑120 (MQ=255)
                                                                             ggAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGTAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCt                                                                                                  >  1:187605/1‑120 (MQ=255)
                                                                             ggAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGTAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCt                                                                                                  <  2:127399/120‑1 (MQ=255)
                                                                             ggAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGTAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCt                                                                                                  <  2:187605/120‑1 (MQ=255)
                                                                              gAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGAtt                                                                    >  2:230338/1‑149 (MQ=255)
                                                                              gAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGAtt                                                                    <  1:241121/149‑1 (MQ=255)
                                                                               aaaGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTg                                                                              >  4:123970/1‑138 (MQ=255)
                                                                               aaaGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTg                                                                              <  3:123970/138‑1 (MQ=255)
                                                                                     ggCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGaaa                                                             >  1:323022/1‑149 (MQ=255)
                                                                                     ggCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGaaa                                                             <  2:425268/149‑1 (MQ=255)
                                                                                      gCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAAt                                                            >  1:287716/1‑149 (MQ=255)
                                                                                      gCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAAt                                                            <  3:121020/149‑1 (MQ=255)
                                                                                        cGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATgg                                                          <  4:39967/149‑1 (MQ=255)
                                                                                         gCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGa                                                         >  2:319290/1‑149 (MQ=255)
                                                                                          cTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGaa                                                        <  2:615188/149‑1 (MQ=255)
                                                                                             cAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTaa                                                                          <  2:492642/128‑1 (MQ=255)
                                                                                             cAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTaa                                                                          >  1:492642/1‑128 (MQ=255)
                                                                                                  ccGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACttt                                                <  1:220668/149‑1 (MQ=255)
                                                                                                   cGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGACACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTg                                               >  2:565538/1‑149 (MQ=255)
                                                                                                    ggCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAGTCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGc                                              <  2:287716/149‑1 (MQ=255)
                                                                                                     gCTTTGGTCAGCGCGTTGAACGGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGcc                                             <  4:140898/149‑1 (MQ=255)
                                                                                                     gCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGcc                                             <  4:344235/149‑1 (MQ=255)
                                                                                                       tttGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTAtt                                                                               >  3:253920/1‑113 (MQ=255)
                                                                                                       tttGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTAtt                                                                               <  4:253920/113‑1 (MQ=255)
                                                                                                          ggTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCtatta                                        <  1:500456/149‑1 (MQ=255)
                                                                                                            tCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATtata                                      >  4:549371/1‑149 (MQ=255)
                                                                                                            tCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATtata                                      >  4:160654/1‑149 (MQ=255)
                                                                                                                 gcgTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCAt                                                                                                                                   >  2:505873/1‑51 (MQ=255)
                                                                                                                 gcgTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCAt                                                                                                                                   <  1:505873/51‑1 (MQ=255)
                                                                                                                   gTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCa                               >  4:578734/1‑149 (MQ=255)
                                                                                                                       aaCAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCg                           >  2:86429/1‑149 (MQ=255)
                                                                                                                           gggTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGa                                                         >  4:305508/1‑115 (MQ=255)
                                                                                                                           gggTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGa                                                         <  3:305508/115‑1 (MQ=255)
                                                                                                                             gTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCaaaa                     >  1:355550/1‑149 (MQ=255)
                                                                                                                             gTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCaaaa                     >  2:596800/1‑149 (MQ=255)
                                                                                                                                       cGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGa           >  3:230943/1‑149 (MQ=255)
                                                                                                                                        gACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGAc          >  2:675702/1‑149 (MQ=255)
                                                                                                                                             tGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGACAACGa     <  2:323022/149‑1 (MQ=255)
                                                                                                                                               ggCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGACAACGACt   <  1:230338/149‑1 (MQ=255)
                                                                                                                                                gCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGACAACGACtt  >  1:114784/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTACGATTTCAGCCAGTTGATCCAGGCGAGGATCAGGCATTGGTACGACGCCTGTGTTCGGCTCAATGGTGCAGAATGGAAAGTTGGCCGCTTCAATACCGGCTTTGGTCAGCGCGTTGAACAGGGTAGATTTCCCGACGTTGGGCAAACCGACGATACCGCATTTGAATCCCATGATTTAAATCACCTTAATATCTTAATAATCAACCTGTTATTGCTAACAGATTGCAGAAATGGAAATAACTTTGCCTATTATACACGGCACTCGGCAAAAATGCCGCAGACAACGACTT  >  NC_000913/1257638‑1257930

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