Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,259 G→C 7.8% A21G (GCT→GGT) ‡ yfiE ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,2590GC 7.8% 299.6 / 13.9 91A21G (GCT→GGT) ‡yfiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (43/40);  new base C (4/3);  total (47/44)
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.97e-01

TTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCA  >  NC_000913/2715111‑2715404
                                                                                                                                                    |                                                                                                                                                  
ttCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAg                                                                                                                                                    >  1:358431/1‑149 (MQ=255)
   aTAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCAACGTCGATTGTGTACAGCACAATTTTTGCGAAGCTc                                                                                                                                                >  1:24447/1‑149 (MQ=255)
          tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGGATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaa                                                                                                                                         <  2:708619/149‑1 (MQ=255)
          tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCaaaaa                                                                                                                                         <  1:172047/149‑1 (MQ=255)
           gCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAg                                                                                                                                        >  3:29360/1‑149 (MQ=255)
             ggCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaattg                                                                                                                                      <  3:255091/149‑4 (MQ=255)
             ggCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGa                                                                                                                                      <  4:383217/149‑1 (MQ=255)
              gcagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCACAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                                     <  4:192909/149‑1 (MQ=255)
                agcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAcc                                                                                                                                   <  1:667920/149‑1 (MQ=255)
                   aaCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctc                                                                                                                                <  1:830914/149‑1 (MQ=255)
                    aCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctct                                                                                                                               >  2:87834/1‑149 (MQ=255)
                      ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                             >  1:87762/1‑149 (MQ=255)
                                cgcgGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGtt                                                                                                                   <  2:690906/149‑1 (MQ=255)
                                       aaGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                                                                                                            >  1:695295/1‑149 (MQ=255)
                                        aGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTCAATATGAAAAGTCACCGTCGATTATGTACAGCACAATTTTTGAGAAGCTCACAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGt                                                                                                           <  3:742601/149‑1 (MQ=255)
                                          gCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaa                                                                                                         >  4:458691/1‑149 (MQ=255)
                                          gCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaa                                                                                                         >  1:719004/1‑149 (MQ=255)
                                          gCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaa                                                                                                         >  3:337665/1‑149 (MQ=255)
                                            acaTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaata                                                                                                       >  2:595149/1‑149 (MQ=255)
                                              aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaa                                                                                                     >  1:669838/1‑149 (MQ=255)
                                                   tCGGCCAATTTTCTCAAATAACTGGGCTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                                <  4:198527/149‑1 (MQ=255)
                                                   tCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGgc                                                                                                <  4:337665/149‑1 (MQ=255)
                                                      gCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCa                                                                                             >  4:706674/1‑149 (MQ=255)
                                                         gaTTTTCTCACATAACGGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTATACAGCACAATTTTTGCGAAGCACGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAGCGGCGCAGag                                                                                          <  2:289318/148‑2 (MQ=255)
                                                          aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATc                                                                                         <  1:469790/149‑1 (MQ=255)
                                                           ttttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAAAGGCGCAGATcc                                                                                        >  4:449269/1‑149 (MQ=255)
                                                             ttCTCAAATAACTGGAATGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAAATTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCAt                                                                                      >  4:219620/1‑149 (MQ=255)
                                                               ctcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATaa                                                                                    <  2:872257/149‑1 (MQ=255)
                                                                tcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                   >  3:296395/1‑149 (MQ=255)
                                                                      aaCTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCaaa                                                                             <  2:87762/149‑1 (MQ=255)
                                                                         tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                          >  2:813499/1‑149 (MQ=255)
                                                                          ggACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTa                                                                         <  1:810035/149‑1 (MQ=255)
                                                                          ggACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTa                                                                         >  3:216065/1‑149 (MQ=255)
                                                                           gACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTAt                                                                        <  3:458691/149‑1 (MQ=255)
                                                                            gcTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATc                                                                       <  2:708564/148‑1 (MQ=255)
                                                                            aCTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATc                                                                       <  2:887287/149‑1 (MQ=255)
                                                                              tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGa                                                                     >  1:274848/1‑149 (MQ=255)
                                                                                aaaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaa                                                                   >  3:443625/1‑149 (MQ=255)
                                                                                aaaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaa                                                                   >  1:315353/1‑149 (MQ=255)
                                                                                aaaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTCCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTAGCGaaa                                                                   <  2:190120/149‑1 (MQ=255)
                                                                                 aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                  >  4:80311/1‑149 (MQ=255)
                                                                                 aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCATCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                  <  4:571819/149‑1 (MQ=255)
                                                                                   aCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                        <  1:730041/125‑1 (MQ=255)
                                                                                   aCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc                                                                                        >  2:730041/1‑125 (MQ=255)
                                                                                    cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                               <  4:678162/149‑1 (MQ=255)
                                                                                    cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                               <  4:29360/149‑1 (MQ=255)
                                                                                    cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                               <  2:785890/149‑1 (MQ=255)
                                                                                     tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                              >  4:743555/1‑149 (MQ=255)
                                                                                     tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                              >  1:244438/1‑149 (MQ=255)
                                                                                     tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                              >  3:520296/1‑149 (MQ=255)
                                                                                        tGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGaa                                                           >  2:217514/1‑149 (MQ=255)
                                                                                         tcTCAAGCTTCATAATAAGAAAAGTCACCCTCCATTTGTTACATCCCAATTTTTGCGAAACTCCCAAAAAGGAACCCTTTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAAt                                                          <  1:631071/148‑1 (MQ=255)
                                                                                          cTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATc                                                         <  1:135736/149‑1 (MQ=255)
                                                                                           tCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCt                                                        >  3:653152/1‑149 (MQ=255)
                                                                                               gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                                    <  1:295955/149‑1 (MQ=255)
                                                                                               gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                                    <  3:706674/149‑1 (MQ=255)
                                                                                               gctgctGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGaa                                                    <  2:263859/149‑1 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAATCCCAAACCTATCGAAAATAGCGAATCTAGAAtata                                                <  2:151003/149‑1 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAAtata                                                >  2:867068/1‑149 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAAtata                                                >  2:47674/1‑149 (MQ=255)
                                                                                                   ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAAtata                                                >  2:62967/1‑149 (MQ=255)
                                                                                                    tgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaa                                               >  1:464533/1‑149 (MQ=255)
                                                                                                       atatGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                            <  1:595947/149‑1 (MQ=255)
                                                                                                         atGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAca                                          >  1:177100/1‑149 (MQ=255)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                       >  2:213955/1‑134 (MQ=255)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                       <  1:213955/134‑1 (MQ=255)
                                                                                                           gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcatt                                        <  2:358431/149‑1 (MQ=255)
                                                                                                                      gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                             >  2:635911/1‑149 (MQ=255)
                                                                                                                      gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                             <  1:595149/149‑1 (MQ=255)
                                                                                                                       tCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaaa                            >  4:130327/1‑149 (MQ=255)
                                                                                                                         gATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATg                          >  1:795099/1‑149 (MQ=255)
                                                                                                                            tgtgTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTc                       <  1:783565/149‑1 (MQ=255)
                                                                                                                            tgtgTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTc                       <  1:867068/149‑1 (MQ=255)
                                                                                                                                tACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg                   >  4:131159/1‑149 (MQ=255)
                                                                                                                                   aGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg                >  4:609044/1‑149 (MQ=255)
                                                                                                                                   aGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg                <  4:578826/149‑1 (MQ=255)
                                                                                                                                   aGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg                >  2:179818/1‑149 (MQ=255)
                                                                                                                                    gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg               <  2:669838/149‑1 (MQ=255)
                                                                                                                                     cacaATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg              >  1:711866/1‑149 (MQ=255)
                                                                                                                                      acaATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGt             <  3:743555/149‑1 (MQ=255)
                                                                                                                                        aaTTTTTGCGACGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTAc           >  3:570518/1‑149 (MQ=255)
                                                                                                                                        aaTTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTAc           <  2:719004/149‑1 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttt                                  <  2:333653/124‑1 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc                      <  2:322399/136‑1 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc                      >  1:322399/1‑136 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg         >  1:683005/1‑149 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg         >  3:570731/1‑149 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg         >  3:644766/1‑149 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttt                                  >  1:333653/1‑124 (MQ=255)
                                                                                                                                          tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg         >  3:427926/1‑149 (MQ=255)
                                                                                                                                               gCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTAc    <  1:87834/149‑1 (MQ=255)
                                                                                                                                                 gAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCa  >  2:169593/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCA  >  NC_000913/2715111‑2715404

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