Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,990,737 C→T 6.2% Q663* (CAA→TAA)  PP_1780 → mannosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,990,7370CT 6.2% 195.2 / 2.5 64Q663* (CAA→TAA) PP_1780mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GATGTTTGAGCGCAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATGGGGTTGCC  >  NC_002947/1990589‑1990885
                                                                                                                                                    |                                                                                                                                                    
gATGTTTGAGCGCAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGc                                                                                                                                                      <  5:230707/149‑1 (MQ=255)
      tGAGCGCAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGAc                                                                                                                                                 <  5:206033/148‑1 (MQ=255)
       gAGCGCAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGAc                                                                                                                                                 >  6:206033/1‑147 (MQ=255)
         gcgcAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGc                                                                                                                                             <  2:184511/149‑1 (MQ=255)
           gcAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGcc                                                                                                                                            >  4:229399/1‑148 (MQ=255)
           gcAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCa                                                                                                                                           <  2:263033/149‑1 (MQ=255)
             aaaGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCatc                                                                                                                                         >  8:50842/1‑149 (MQ=255)
             aaaGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGAACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCAatc                                                                                                                                         >  6:156552/1‑149 (MQ=255)
                 gCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCatcgatc                                                                                                                                     <  6:335176/149‑1 (MQ=255)
                   ccTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGt                                                                                                                                   >  3:297288/1‑149 (MQ=255)
                    cTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTg                                                                                                                                  <  2:18268/149‑1 (MQ=255)
                    cTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTg                                                                                                                                  <  6:138314/149‑1 (MQ=255)
                     tGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTgg                                                                                                                                 >  1:141566/1‑149 (MQ=255)
                       aCGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTgg                                                                                                                                 <  3:229399/147‑1 (MQ=255)
                           gctgctTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGAc                                                                                                                            >  7:81178/1‑148 (MQ=255)
                                ttGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGc                                                                                                                      >  5:61327/1‑149 (MQ=255)
                                       cTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATgg                                                                                                               <  6:61327/149‑1 (MQ=255)
                                       cTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATgg                                                                                                               <  4:297288/149‑1 (MQ=255)
                                             cGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGAcc                                                                                                         <  6:166844/149‑1 (MQ=255)
                                                         cAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATcg                                                                                             >  3:123081/1‑149 (MQ=255)
                                                             gcCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGATGAAAGCCCTATGGACGACCAAGGCATTATCGCGtt                                                                                         <  1:53790/149‑1 (MQ=255)
                                                             gcCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGtt                                                                                         <  8:81178/149‑1 (MQ=255)
                                                                      gTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGTAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTc                                                                                >  5:144837/1‑149 (MQ=255)
                                                                      gTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGTAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTc                                                                                >  5:150463/1‑149 (MQ=255)
                                                                       tCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCaa                                                                                                       <  6:205792/125‑1 (MQ=255)
                                                                       tCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCaa                                                                                                       >  5:205792/1‑125 (MQ=255)
                                                                          aGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGc                                                                                                    <  2:103683/125‑1 (MQ=255)
                                                                          aGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGACCTATGGACGACCAAGGc                                                                                                    >  1:103683/1‑125 (MQ=255)
                                                                              gTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCcaca                                                                        >  5:304714/1‑149 (MQ=255)
                                                                               tACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATc                                                                                              <  6:52733/126‑1 (MQ=255)
                                                                               tACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATc                                                                                              >  5:52733/1‑126 (MQ=255)
                                                                                 cATCAAGACGTTCACCAACCCACATAACCGCGTAGATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGc                                                                     <  5:156552/149‑1 (MQ=255)
                                                                                  aTCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCt                                                                    >  3:354617/1‑149 (MQ=255)
                                                                                  aTCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCt                                                                    >  2:293242/1‑149 (MQ=255)
                                                                                    cAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTgg                                                                  >  2:187587/1‑149 (MQ=255)
                                                                                        aCGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGtgct                                                              <  2:141566/149‑1 (MQ=255)
                                                                                          gTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGtgctgc                                                            >  7:260445/1‑149 (MQ=255)
                                                                                              cccAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAAc                                                        >  2:268461/1‑149 (MQ=255)
                                                                                                 aaCCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAAcgcg                                                     >  1:11316/1‑149 (MQ=255)
                                                                                                   cccACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCgtg                                                   >  2:355135/1‑149 (MQ=255)
                                                                                                      acaTAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCgtggt                                                 >  3:95466/1‑148 (MQ=255)
                                                                                                        aTAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGaa                                              <  1:187587/149‑1 (MQ=255)
                                                                                                            ccGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTa                                                                                                <  2:252598/95‑1 (MQ=255)
                                                                                                            ccGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTa                                                                                                >  1:252598/1‑95 (MQ=255)
                                                                                                            ccGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCt                                          >  6:364901/1‑149 (MQ=255)
                                                                                                             cgcgTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTTAAGGCtt                                         <  4:95466/149‑1 (MQ=255)
                                                                                                                gTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCaa                                      <  7:50842/149‑1 (MQ=255)
                                                                                                                  aCATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACt                                    >  8:180374/1‑149 (MQ=255)
                                                                                                                  aCATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACt                                    >  5:131636/1‑149 (MQ=255)
                                                                                                                     tGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTAc                                 <  4:123081/149‑1 (MQ=255)
                                                                                                                      gAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTAc                                 >  7:33340/1‑148 (MQ=255)
                                                                                                                         ttGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGtgc                                                               >  8:271623/1‑115 (MQ=255)
                                                                                                                         ttGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGtgc                                                               <  7:271623/115‑1 (MQ=255)
                                                                                                                                    ccTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAgcgc                  <  8:33340/149‑1 (MQ=255)
                                                                                                                                       ctcGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTgg               >  8:220888/1‑149 (MQ=255)
                                                                                                                                        tcGGCTGATCGGTAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTg                <  6:150463/147‑1 (MQ=255)
                                                                                                                                        tcGGCTGATCGGTAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTg                <  6:144837/147‑1 (MQ=255)
                                                                                                                                        tcGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTgg               >  8:294646/1‑148 (MQ=255)
                                                                                                                                        tcGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTgg               >  5:287824/1‑148 (MQ=255)
                                                                                                                                        tcGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTg                <  4:354617/147‑1 (MQ=255)
                                                                                                                                             tGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATgg         >  5:335635/1‑149 (MQ=255)
                                                                                                                                               aTCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCg                                                                              <  7:244577/78‑1 (MQ=255)
                                                                                                                                               aTCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCg                                                                              >  8:244577/1‑78 (MQ=255)
                                                                                                                                               aTCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATggg        <  1:293242/148‑1 (MQ=255)
                                                                                                                                               aTCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATggg        <  5:54616/148‑1 (MQ=255)
                                                                                                                                                    cAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATGGGGTTGcc  <  7:180374/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATGTTTGAGCGCAAAGGCCCTGACGTGCTGCTTGAGGCCTTTGCCGAGGCTTTCTCCAGCGCCGATAACGTCGAGCTGTACATCAAGACGTTCCCCAACCCACATAACCGCGTACATGAATTGATCGCCAACCTCTCGGCTGATCGGCAAGACATGCCATCGATCGTGGTTGACGAAAGCCCTATGGACGACCAAGGCATTATCGCGTTGTATCGTTCCGCCCACACGCTGGTGCTGCCAACGCGTGGTGAAGGCTTCAACTTACCAGCGGCTGAAGCGCTGGCAATGGGGTTGCC  >  NC_002947/1990589‑1990885

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