Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,184,487 Δ1 bp 100% intergenic (‑21/‑353) ribB ← / → yqiC 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/BMFP family putative fusogenic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,184,4840A.72.5% 184.8 / 125.4 102intergenic (‑18/‑356)ribB/yqiC3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/BMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base . (35/39);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 113 total aligned reads displayed.

CGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTC  >  NC_000913/3184385‑3184581
                                                                                                   |                                                                                                 
cgcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTa                                                                                                   >  1:496934/1‑100 (MQ=255)
cgcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTa                                                                                                   >  1:497037/1‑100 (MQ=255)
 gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaa                                                                                                  >  2:1570513/1‑100 (MQ=255)
  cgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaa                                                                                                 <  2:496935/100‑1 (MQ=255)
   gTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaat                                                                                                <  2:1835644/100‑2 (MQ=255)
     ccTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaatta                                                                                              <  1:1330382/100‑4 (MQ=255)
     ccTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaatta                                                                                              >  2:1676032/1‑97 (MQ=255)
      cTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAACGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaattat                                                                                             >  1:636080/1‑96 (MQ=255)
      cTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaattat                                                                                             >  1:1241892/1‑96 (MQ=255)
      cTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTaaattat                                                                                             >  1:1716695/1‑96 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  2:417156/100‑1 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  1:2110369/100‑1 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  2:1667467/100‑1 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  2:2735438/100‑1 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  2:272135/100‑1 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          >  1:1076292/1‑100 (MQ=255)
        tCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          >  1:89525/1‑100 (MQ=255)
        tCACACAGCGCAGCCAGTGCATTTTAAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATgg                                                                                          <  1:483061/100‑1 (MQ=255)
          aCGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGtt                                                                                        <  1:494454/100‑1 (MQ=255)
             cagcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTAcc                                                                                     >  2:722729/1‑100 (MQ=255)
             cagcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTAcc                                                                                     <  1:1581262/100‑1 (MQ=255)
              agggcagcCAGTGCATTTTGATAACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTAATGGTAAAAAAACCTCACTAAAATTATGGTTAcc                                                                                     <  1:2119552/97‑1 (MQ=255)
              agcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTAcc                                                                                     >  1:765753/1‑100 (MQ=255)
               gcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCa                                                                                    <  1:2104202/100‑1 (MQ=255)
                 gcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGaa                                                                                 <  2:1425877/100‑1 (MQ=255)
                    gcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATc                                                                               >  1:2168122/1‑100 (MQ=255)
                    gcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCa                                                                              >  2:1314483/1‑100 (MQ=255)
                      cAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAg                                                                             >  1:1032283/1‑100 (MQ=255)
                       aGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAggg                                                                           <  2:685344/100‑1 (MQ=255)
                       aGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAggg                                                                           >  2:2556174/1‑100 (MQ=255)
                        gTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGc                                                                          <  2:290866/100‑1 (MQ=255)
                         tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCa                                                                         >  1:824156/1‑100 (MQ=255)
                         tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCa                                                                         <  2:2828927/100‑1 (MQ=255)
                         tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCAAAAAAAAAAAAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGc                                                                          <  1:2512735/99‑1 (MQ=255)
                             ttttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCt                                                                     <  2:215615/100‑1 (MQ=255)
                               ttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTa                                                                   >  1:2381690/1‑100 (MQ=255)
                               ttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTa                                                                   >  2:2257672/1‑100 (MQ=255)
                               ttCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTa                                                                   >  1:1293469/1‑100 (MQ=255)
                                 cAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAg                                                                  <  2:1246479/100‑1 (MQ=255)
                                 cAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAgg                                                                 >  1:350298/1‑100 (MQ=255)
                                  aaCACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAgg                                                                 <  2:710506/100‑1 (MQ=255)
                                  aaCACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGa                                                                >  2:1857186/1‑100 (MQ=255)
                                   acacGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAg                                                               >  1:2252396/1‑100 (MQ=255)
                                     acGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAg                                                               <  2:649660/98‑1 (MQ=255)
                                       gTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGTTTACCAGAATCAGGGCATTCTTAGGAGTgcc                                                           <  2:1960111/100‑3 (MQ=255)
                                       gTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTggc                                                           >  2:1979674/1‑100 (MQ=255)
                                        ttCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTggcg                                                          <  2:1050561/100‑1 (MQ=255)
                                          cGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTggcggc                                                        <  1:1062123/100‑1 (MQ=255)
                                          cGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTggcggc                                                        <  2:2634018/100‑1 (MQ=255)
                                          cGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTggcggc                                                        >  2:2901603/1‑100 (MQ=255)
                                           gAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCCTGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCa                                                       >  1:1554885/1‑100 (MQ=255)
                                             aaGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCata                                                     >  2:1923991/1‑100 (MQ=255)
                                             aaGGCGTACCAAAAGAGGAAAGTAGCGTCTGATCCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCat                                                      >  1:1741493/1‑100 (MQ=255)
                                             aaGGCGTACCAAAAGAGGAAAGTAGCGTCTGATACATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCat                                                      >  1:1741494/1‑100 (MQ=255)
                                              aGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCatat                                                    <  1:2358452/100‑1 (MQ=255)
                                               ggCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCatat                                                    >  2:2097750/1‑100 (MQ=255)
                                               ggCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCAtata                                                   >  2:257521/1‑100 (MQ=255)
                                               ggCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCAtata                                                   >  2:1807354/1‑100 (MQ=255)
                                                gCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCAtata                                                   <  2:613229/100‑1 (MQ=255)
                                                gCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCAtata                                                   >  2:1357869/1‑100 (MQ=255)
                                                  gTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGcc                                                <  2:191350/100‑1 (MQ=255)
                                                    aCCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGcc                                                <  1:1570509/99‑1 (MQ=255)
                                                    aCCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCaa                                              >  2:2881101/1‑100 (MQ=255)
                                                        aaaGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAAt                                           <  1:2011086/100‑1 (MQ=255)
                                                        aaaGAGGAAAGTAGCGTCTGATTCAGGGTAAAAAAACCTCACT‑CAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAtata                                          <  1:2883007/100‑1 (MQ=255)
                                                             ggAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACg                                       <  2:765755/99‑1 (MQ=255)
                                                             ggAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTg                                     >  2:1646621/1‑100 (MQ=255)
                                                               aaaGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAg                                   >  1:2014361/1‑100 (MQ=255)
                                                                 aGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCg                                 >  2:1787464/1‑100 (MQ=255)
                                                                 aGTAGCGTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCg                                 <  1:1345734/100‑1 (MQ=255)
                                                                    aGCGTCTGATTCATGGTAAAAAAACCTCACT‑AACTTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGt                              >  1:1751805/1‑100 (MQ=255)
                                                                     gaggcTGATTCATTGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGt                              <  1:788523/95‑1 (MQ=255)
                                                                       gTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTcc                            >  1:508686/1‑100 (MQ=255)
                                                                       gTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCa                           >  2:1381482/1‑100 (MQ=255)
                                                                       gTCTGATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCa                           <  2:2550113/100‑1 (MQ=255)
                                                                          tGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATg                         >  2:1372343/1‑100 (MQ=255)
                                                                           gATTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGcc                       >  2:188535/1‑100 (MQ=255)
                                                                            atccATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGccc                      >  2:2263867/4‑100 (MQ=255)
                                                                            aTTCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGccc                      <  2:1859563/100‑1 (MQ=255)
                                                                             ttCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCg                     <  2:350300/100‑1 (MQ=255)
                                                                              tCATGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCg                     <  2:476998/99‑1 (MQ=255)
                                                                                aTGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACa                  <  2:2499903/100‑1 (MQ=255)
                                                                                 tGGTAAAAAAACCTCACTAAAATTAGGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACa                  <  1:260580/100‑1 (MQ=255)
                                                                                 tGGTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAg                 >  2:1798841/1‑100 (MQ=255)
                                                                                  ggTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAg                 >  2:2211670/1‑100 (MQ=255)
                                                                                   gTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGaa               >  2:1425505/1‑100 (MQ=255)
                                                                                   gTAAAAAAACCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGaa               <  1:806313/100‑1 (MQ=255)
                                                                                     aaaaaaaCCTCACT‑AAATTATGGTTACAAGAATCAGGGAAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATc             >  1:607385/1‑100 (MQ=255)
                                                                                        aaaaCCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGtt          <  1:2854667/100‑1 (MQ=255)
                                                                                          aaCCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAc        >  2:1635621/1‑100 (MQ=255)
                                                                                          aaCCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAc        <  2:2803796/100‑1 (MQ=255)
                                                                                          aaCCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAc        <  1:35302/100‑1 (MQ=255)
                                                                                           aCCTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTAct       <  1:180063/100‑1 (MQ=255)
                                                                                           aCCTCACT‑AAATTATGGTAACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGTCAGAATCGTTAc        <  1:2134703/99‑1 (MQ=255)
                                                                                            ccTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTActc      <  1:2694502/100‑1 (MQ=255)
                                                                                             cTCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTActct     <  1:1340858/100‑1 (MQ=255)
                                                                                              tCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTActctc    >  2:2653707/1‑100 (MQ=255)
                                                                                              tCACT‑AAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTActctc    >  2:2213471/1‑100 (MQ=255)
                                                                                                 cTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc  >  1:2455619/1‑100 (MQ=255)
                                                                                                 cTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACtctctc  >  2:946813/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CGCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTC  >  NC_000913/3184385‑3184581

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