Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,000,100 A→G 6.5% C140R (TGC→CGC)  PP_2623 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,000,1000AG 6.5% 170.3 / 2.0 62C140R (TGC→CGC) PP_2623hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/24);  new base G (2/2);  total (35/27)
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.57e-01

GGACATGGCGCCAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTGGAACT  >  NC_002947/2999951‑3000243
                                                                                                                                                     |                                                                                                                                               
ggACATGGCGCCAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCa                                                                                                                                                 <  2:1807667/150‑1 (MQ=255)
 gACATGGCGCCAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAg                                                                                                                                                >  1:158615/1‑150 (MQ=255)
        cgcCAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGt                                                                                                                                         <  1:973215/150‑1 (MQ=255)
           cAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTcag                                                                                                                                      <  2:933180/150‑1 (MQ=255)
           cAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTcag                                                                                                                                      >  2:1664577/1‑150 (MQ=255)
            aaCGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCGGGTAACGTcagc                                                                                                                                     >  2:1478271/1‑150 (MQ=255)
                   gCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCAttt                                                                                                                              >  2:520846/1‑150 (MQ=255)
                     aaTGTAAATTAGCAACTTGATGGCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTc                                                                                                                            <  2:745307/149‑1 (MQ=255)
                      aTGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGTCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCg                                                                                                                           >  1:1799505/1‑150 (MQ=255)
                      aTGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCg                                                                                                                           >  1:932590/1‑150 (MQ=255)
                        gTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGaa                                                                                                                         <  1:572032/150‑1 (MQ=255)
                           aCTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCagcagc                                                                                                                                  <  2:1456779/138‑1 (MQ=255)
                           aCTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCagcagc                                                                                                                                  >  1:1456779/1‑138 (MQ=255)
                             ttAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAgtcgt                                                                                                                    >  1:1859946/1‑150 (MQ=255)
                              tAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTGACTTCAGCAGCATTTTCGAAgtcgtc                                                                                                                   <  1:989700/150‑1 (MQ=255)
                                  aaCTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTcag                                                                                                                                      <  2:189314/127‑1 (MQ=255)
                                  aaCTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTcag                                                                                                                                      >  1:189314/1‑127 (MQ=255)
                                   aCTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCa                                                                                                              <  1:1110778/150‑1 (MQ=255)
                                         tGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCct                                                                                                        <  1:1406511/150‑1 (MQ=255)
                                           tCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCctct                                                                                                      >  1:1688735/1‑150 (MQ=255)
                                           tCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCctct                                                                                                      >  1:710727/1‑150 (MQ=255)
                                              ggggAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTa                                                                                                   <  2:1375198/150‑1 (MQ=255)
                                                     agGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACAt                                                                                            <  2:1165310/150‑1 (MQ=255)
                                                       gggCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCGGGTAACGTcagc                                                                                                                                     <  2:777370/107‑1 (MQ=255)
                                                       gggCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCGGGTAACGTcagc                                                                                                                                     >  1:777370/1‑107 (MQ=255)
                                                                 tGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCaa                                                                                <  1:520846/150‑1 (MQ=255)
                                                                        gCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGAt                                                                         <  2:1206172/150‑1 (MQ=255)
                                                                            aCCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTtc                                                                     >  1:1803412/1‑150 (MQ=255)
                                                                             ccATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTtct                                                                    >  2:134896/1‑150 (MQ=255)
                                                                              cATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAAc                                                                                                                                           >  2:1651049/1‑78 (MQ=255)
                                                                              cATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAAc                                                                                                                                           <  1:1651049/78‑1 (MQ=255)
                                                                               aTGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTCCTGACATGCAAGATTTCAAGCTTGATGTTCTcc                                                                  >  1:1353835/1‑150 (MQ=255)
                                                                               aTGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTcc                                                                  >  1:1699795/1‑150 (MQ=255)
                                                                                  gCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTg                                                               >  1:1803161/1‑150 (MQ=255)
                                                                                   ccGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTGATGCTTGATGTTCTCCCTGa                                                              >  1:1531311/1‑150 (MQ=255)
                                                                                   ccGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGa                                                              >  2:3690/1‑150 (MQ=255)
                                                                                    cGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTACATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGACGGTTTCCCTGaa                                                             >  1:1239848/1‑150 (MQ=255)
                                                                                     gCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAg                                                            <  1:759413/150‑1 (MQ=255)
                                                                                         ccGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATATTCTCCCTGAAGTACg                                                        >  2:640291/1‑150 (MQ=255)
                                                                                           ggCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCGGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTg                                                      <  1:1478271/150‑1 (MQ=255)
                                                                                           ggCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTg                                                      <  2:1688735/150‑1 (MQ=255)
                                                                                                  ccccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCCGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGa                                               <  1:1565802/150‑1 (MQ=255)
                                                                                                  ccccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCctc                                                                                                       <  2:1182176/94‑1 (MQ=255)
                                                                                                  ccccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCctc                                                                                                       >  1:1182176/1‑94 (MQ=255)
                                                                                                   cccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGAc                                              >  2:1132973/1‑150 (MQ=255)
                                                                                                   cccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGAc                                              >  2:1653591/1‑150 (MQ=255)
                                                                                                   cccGAATTGGCCATAACCCGCGGTGTAAACATCACGGCAAATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAAACGACGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTAGAc                                              >  1:627527/1‑150 (MQ=255)
                                                                                                    ccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCa                                                                                                                                                 <  1:701077/50‑1 (MQ=255)
                                                                                                    ccGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCa                                                                                                                                                 >  2:701077/1‑50 (MQ=255)
                                                                                                         ttGGCGATAAGCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGGCAGCAGCGTCTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATATTCTCCCTGAATTCCGCGCGGTCGACCAGAAg                                        >  2:552385/1‑150 (MQ=255)
                                                                                                         ttGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAg                                        >  1:1464739/1‑150 (MQ=255)
                                                                                                              cATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCa                                   <  1:699569/150‑1 (MQ=255)
                                                                                                                    ccGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGt                                                  <  2:1661411/129‑1 (MQ=255)
                                                                                                                    ccGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGt                                                  >  1:1661411/1‑129 (MQ=255)
                                                                                                                    ccGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGCTCGACCAGAAGTTTCAGGCTgc                             >  1:1464917/1‑150 (MQ=255)
                                                                                                                          tgtAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCatg                       >  1:1889846/1‑150 (MQ=255)
                                                                                                                              aaCATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTAACAGGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCAtgatga                   <  1:230925/150‑1 (MQ=255)
                                                                                                                                    cTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCATCCCTCATGATAGGGAACACTACAATCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAg             <  1:134896/150‑1 (MQ=255)
                                                                                                                                     ttGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAgc            >  1:1812986/1‑150 (MQ=255)
                                                                                                                                         aGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGCTTCAAgctgcc        >  2:513205/1‑149 (MQ=255)
                                                                                                                                          gATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGctgctg       <  1:1664577/150‑1 (MQ=255)
                                                                                                                                           aTACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTgg      >  2:532513/1‑150 (MQ=255)
                                                                                                                                             aCCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTGGaa    <  1:1132973/150‑1 (MQ=255)
                                                                                                                                              ccGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTGGAAc   <  1:1653591/150‑1 (MQ=255)
                                                                                                                                               cGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTGGAACt  >  2:1293004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGACATGGCGCCAACGGAGGCCATGTAACTTAGCAACTTGATGTCAGGGGAAGAGGGGCCGAAGTTGTAGTTGCCTACCATGGCCGCAACCGGCTCTCCCCCGAATTGGCCATAACCCGCGGTGTAAACATGCTTGTAGATACCGCTGCAGGTAACGTCAGCAGCATTTTCGAAGTCGTCGAGCAGGTCCTCTTTACTGACATGCAAGATTTCAAGCTTGATGTTCTCCCTGAAGTCCGTGCGGTCGACCAGAAGTTTCAGGCTGCGCCATGATGATTCAAGCTGCTGGAACT  >  NC_002947/2999951‑3000243

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