Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 327,213 G→C 6.7% S12* (TCA→TGA)  betA ← choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913327,2130GC 6.7% 378.3 / 14.9 119S12* (TCA→TGA) betAcholine dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (54/57);  new base C (3/5);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

Only 100 of 123 total aligned reads displayed.

GTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTG  >  NC_000913/327068‑327362
                                                                                                                                                 |                                                                                                                                                     
gtaccctttACCCTGTAGCGGGAAAGCCACGGCCTCCGGCATCTGGGGGGGGAAGGCAAAGCGATAGGCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGcc                                                                                                                                                   <  2:1115603/144‑1 (MQ=255)
    cggttACCCTGTATGGGGAATGCCAGGGCAGCGGGCATGTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCTCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCTGCCGGCa                                                                                                                                               <  2:375015/147‑1 (MQ=255)
     gTTTACCCTGTAGCGGGAATTCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAc                                                                                                                                              <  2:1334043/150‑1 (MQ=255)
      tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                             <  1:1192259/150‑1 (MQ=255)
      tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaac                                                                                                                                             <  1:576021/150‑3 (MQ=255)
      tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaac                                                                                                                                             <  1:144437/150‑3 (MQ=255)
      tttACCCTGTAGCGGGAATGCCAGCGCAGCGGGCATCTGGGTGCGGAAGGCAACGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                             <  1:1065178/150‑1 (MQ=255)
      cttACCCTGTAGCGGGAACGCCCGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaac                                                                                                                                             <  1:832717/149‑3 (MQ=255)
       ttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCa                                                                                                                                            >  2:473529/1‑150 (MQ=255)
        tACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGCCGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaa                                                                                                                                           >  1:448956/1‑150 (MQ=255)
         aCCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaat                                                                                                                                          <  1:403516/150‑1 (MQ=255)
         aCCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaacgtt                                                                                                                                          <  1:525837/150‑6 (MQ=255)
          cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaata                                                                                                                                         >  2:1693477/1‑150 (MQ=255)
          cccTGAAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCGTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCAGCGCCGGCACCaata                                                                                                                                         >  2:1516448/1‑150 (MQ=255)
           ccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa                                                                                                                                        >  1:505068/1‑150 (MQ=255)
           ccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa                                                                                                                                        >  1:1425061/1‑150 (MQ=255)
           ccTGTAGCGGGAATGCCAGGGCAGCGGGCATCGGGGTGCGGAAGTCAACGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa                                                                                                                                        <  1:596510/150‑1 (MQ=255)
                  cGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTa                                                                                                                                 >  2:1501382/1‑150 (MQ=255)
                   gggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCCGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGAAGTAg                                                                                                                                >  1:284181/1‑150 (MQ=255)
                    ggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGt                                                                                                                               >  2:326559/1‑150 (MQ=255)
                    ggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGATTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCTGCACCAATAATGATGTAGt                                                                                                                               >  2:1177434/1‑150 (MQ=255)
                       aTGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCaa                                                                                                                            <  1:196284/150‑1 (MQ=255)
                          ccAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAtt                                                                                                                         <  1:925837/150‑1 (MQ=255)
                           cAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTg                                                                                                                        <  2:1620227/150‑1 (MQ=255)
                           cAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCAGGCACCAATAATGATGTAGTCAAATTg                                                                                                                        >  1:199643/1‑150 (MQ=255)
                            aGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGc                                                                                                                       <  2:1229832/150‑1 (MQ=255)
                             gggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCa                                                                                                                      <  2:1512664/150‑1 (MQ=255)
                              ggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCaa                                                                                                                     >  1:1629160/1‑150 (MQ=255)
                                cAGCGGGCATCGGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTGCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAAt                                                                                                                   <  2:1478974/150‑1 (MQ=255)
                                 aGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATa                                                                                                                  >  2:1009405/1‑150 (MQ=255)
                                   cGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATaaa                                                                                                                >  1:1689177/1‑150 (MQ=255)
                                     ggCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAAcc                                                                                                              <  1:1126737/150‑1 (MQ=255)
                                       cATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCCCCAATAATGATGTAGTCAAATTGCAAATAAAcctt                                                                                                            >  2:450401/1‑149 (MQ=255)
                                               tGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAg                                                                                                    <  1:1282292/150‑1 (MQ=255)
                                               tGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAg                                                                                                    >  1:956286/1‑150 (MQ=255)
                                                  ggAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAat                                                                                                 >  1:899502/1‑150 (MQ=255)
                                                  ggAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAAGAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAat                                                                                                 >  2:208346/1‑150 (MQ=255)
                                                   gAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAata                                                                                                <  2:480438/150‑1 (MQ=255)
                                                    aaGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAatat                                                                                               <  1:1405188/150‑1 (MQ=255)
                                                     cgTCAAAGCGAGAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGTGGGATACGGTTCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATg                                                                                              <  2:830811/149‑1 (MQ=255)
                                                     aGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATg                                                                                              >  2:1339526/1‑150 (MQ=255)
                                                     aGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATg                                                                                              >  2:1346225/1‑150 (MQ=255)
                                                      gTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATgg                                                                                             >  1:1122404/1‑150 (MQ=255)
                                                         aaaGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACt                                                                                          >  1:765008/1‑150 (MQ=255)
                                                          aaGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTg                                                                                         >  2:594390/1‑150 (MQ=255)
                                                               agaGTCCGGGCCGCCCGCTTCACGCAGCAGCACGGAGGTCCTCGTCTCTGCAGTCAGACGGGAAGCGAGAACGTTGCCTGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGCAtt                                                                                    <  2:88909/148‑1 (MQ=255)
                                                               aTAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAtt                                                                                    >  2:619667/1‑150 (MQ=255)
                                                                 aGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTa                                                                                  >  2:504195/1‑150 (MQ=255)
                                                                    ccGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCTGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATTATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGcc                                                                               <  1:1497992/150‑1 (MQ=255)
                                                                    ccGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGGTTGAATATGGCCTGGAATTTAGcc                                                                               >  1:692783/1‑150 (MQ=255)
                                                                     cGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCa                                                                              <  1:924832/150‑1 (MQ=255)
                                                                        gccgccCGCTTCAATCAGCAGCACGGAGGTATTCGGATCTTCAGGCAGACTTGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAATGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCAtct                                                                           <  1:843697/150‑1 (MQ=255)
                                                                        gccgccCGCTTCAAGCAGCAGCACGGAGGTATTCGGAGCTTCAGGCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCAtct                                                                           <  2:1600485/150‑1 (MQ=255)
                                                                         ccgcccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCAtctc                                                                          <  2:1629160/150‑1 (MQ=255)
                                                                          cgcccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCa                                                                         >  2:210614/1‑150 (MQ=255)
                                                                            cccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAAc                                                                       <  1:1112087/150‑1 (MQ=255)
                                                                             ccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAAAGATGTAGTCAAATTGCAAATAAACCTCCGGGTTAGAATATGGACTGTAATTTATCCATCGCAAcc                                                                      >  2:1525151/1‑150 (MQ=255)
                                                                                cTTCAAGCAGCACCACGGAGGTATTCGGAACTTCAGGCAGCCGGGGAGCGAGAACGTTGCCGGCTGAGCCGGCACCAAAAATGATGAAGTCAAATTGCAAATAAACCACCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtgg                                                                   <  2:1630262/150‑1 (MQ=255)
                                                                                  tCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggat                                                                 <  1:473529/150‑1 (MQ=255)
                                                                                   cAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                                <  2:28925/150‑1 (MQ=255)
                                                                                   cAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                                <  2:29963/150‑1 (MQ=255)
                                                                                     agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                              >  1:1275775/1‑150 (MQ=255)
                                                                                     agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                              <  1:991271/150‑1 (MQ=255)
                                                                                       cagcagcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCACCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACt                                                            <  1:245751/150‑1 (MQ=255)
                                                                                        agcagcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAATGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACtt                                                           >  1:1540897/1‑150 (MQ=255)
                                                                                             caCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCt                                                      >  1:582132/1‑150 (MQ=255)
                                                                                              aCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTg                                                     <  2:899502/150‑1 (MQ=255)
                                                                                              aCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTg                                                     <  2:534537/150‑1 (MQ=255)
                                                                                                ggAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTggg                                                   <  2:1565392/150‑1 (MQ=255)
                                                                                                  aGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGgtg                                                 <  2:765008/150‑1 (MQ=255)
                                                                                                    gTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTa                                               <  1:1339526/150‑1 (MQ=255)
                                                                                                     tATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTaa                                              >  2:1100025/1‑150 (MQ=255)
                                                                                                      tttCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAAc                                             <  2:1275775/149‑1 (MQ=255)
                                                                                                         cGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCCTCTCTACCTGGATGGACTTCACCTGGGTGTAActct                                          >  2:323043/1‑150 (MQ=255)
                                                                                                          ggATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTg                                         <  1:1278807/150‑1 (MQ=255)
                                                                                                             tCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTTACTCTGGAg                                      >  1:637126/1‑150 (MQ=255)
                                                                                                              cTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGTTGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGc                                     <  1:1693477/150‑1 (MQ=255)
                                                                                                              cTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGc                                     >  2:1267689/1‑150 (MQ=255)
                                                                                                              cTTCAGTCAGACGGGTAGCGAGAACATTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTgg                                        >  2:1457457/1‑147 (MQ=255)
                                                                                                               ttCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCg                                    >  2:1596287/1‑150 (MQ=255)
                                                                                                                  agtcagCCGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTAAAATTGCAAATAAAACTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtca                                 <  1:877865/150‑1 (MQ=255)
                                                                                                                    tcagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGTTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtcatc                               <  1:1232501/150‑1 (MQ=255)
                                                                                                                    tcagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtcatc                               <  1:1009405/150‑1 (MQ=255)
                                                                                                                     cagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCg                                    >  2:1045257/1‑144 (MQ=255)
                                                                                                                     cagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCg                                    <  1:1045257/144‑1 (MQ=255)
                                                                                                                      agACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCAc                             <  1:735957/150‑1 (MQ=255)
                                                                                                                          gggTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGTCTGGAATTTAGCCATCTCACCCTGGATGGACGTCACCTGGGTGTACCTCTGGAGCGTCATCACGCCg                         >  1:1175642/1‑150 (MQ=255)
                                                                                                                                cGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcg                   <  2:1187547/150‑1 (MQ=255)
                                                                                                                                cGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGAAGCGTCATCACGCCGTTCTcg                   <  2:448956/150‑1 (MQ=255)
                                                                                                                                   tattggtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGa                >  1:1596991/7‑150 (MQ=255)
                                                                                                                                    aaCGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGAc               <  1:1549203/150‑1 (MQ=255)
                                                                                                                                     ttggtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTTGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGAcc              >  1:136206/5‑150 (MQ=255)
                                                                                                                                     ttggtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGGcc              >  1:1035713/5‑150 (MQ=255)
                                                                                                                                        ttGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCTCGCCGTTCTCGCGACCAAt           >  1:245306/1‑150 (MQ=255)
                                                                                                                                         tGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATg          >  1:246590/1‑150 (MQ=255)
                                                                                                                                          gCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGc         >  2:1021740/1‑150 (MQ=255)
                                                                                                                                           ccGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCGGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACTCCGTTCTCGCGCCCAATGcc        >  2:238344/1‑150 (MQ=255)
                                                                                                                                              gCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGa     >  1:308744/1‑150 (MQ=255)
                                                                                                                                                 gAGCCGGCACCAATAATGATGTAGTCAAATTGCATATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAgtg  <  1:930667/150‑1 (MQ=255)
                                                                                                                                                 gAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAgtg  >  2:1424364/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTG  >  NC_000913/327068‑327362

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