Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,933,707 G→T 9.2% R30R (CGC→CGA)  minD ← MinC‑MinD‑MinE system ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,933,7070GT 9.2% 195.9 / 6.7 65R30R (CGC→CGA) minDMinC‑MinD‑MinE system ATPase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base T (3/3);  total (32/33)
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.25e-01

GCTTGTCCTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATA  >  NC_002947/1933562‑1933851
                                                                                                                                                 |                                                                                                                                                
gcTTGTCCTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCa                                                                                                                                              >  2:842451/1‑150 (MQ=255)
       cTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAg                                                                                                                                             <  2:1048991/144‑1 (MQ=255)
       cTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAg                                                                                                                                             >  1:1048991/1‑144 (MQ=255)
       cTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAg                                                                                                                                       >  2:1165489/1‑150 (MQ=255)
        ttGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAgg                                                                                                                                      >  1:1351431/1‑150 (MQ=255)
         tGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCTCGCACTGCGCGGc                                                                                                                                     <  2:31447/150‑1 (MQ=255)
            tCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCgg                                                                                                                                  >  2:632513/1‑150 (MQ=255)
            tCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCgg                                                                                                                                  <  2:1289494/150‑1 (MQ=255)
                 gCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCTCGCACTGCGCGGccacaa                                                                                                                                <  1:143006/147‑5 (MQ=255)
                 gCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCTCGCACTGCGCGGccacaa                                                                                                                                >  2:143006/1‑143 (MQ=255)
                 gCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCa                                                                                                                             >  2:1220599/1‑150 (MQ=255)
                  ccTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGCCGGTCTTGTGGCCGCGCAGTGCGAGGCGGGTACaaa                                                                                                                            >  2:688751/1‑150 (MQ=255)
                              tGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTgg                                                                                                                <  1:1220599/150‑1 (MQ=255)
                              tGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTgg                                                                                                                <  2:574002/150‑1 (MQ=255)
                               ggCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGt                                                                                                               >  2:886715/1‑150 (MQ=255)
                                gCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTa                                                                                                              <  1:842451/150‑1 (MQ=255)
                                  ttCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGt                                                                                                            >  2:564216/1‑150 (MQ=255)
                                   tCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                           >  2:819241/1‑150 (MQ=255)
                                     gCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGt                                                                                                         <  1:1165489/150‑1 (MQ=255)
                                         ttGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTg                                                                                                     <  2:319878/150‑1 (MQ=255)
                                          ttaCCACGTTGCCGATGTCCTACACCACGCCGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGc                                                                                                    <  2:167938/148‑1 (MQ=255)
                                          tGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGc                                                                                                    <  1:1199386/150‑1 (MQ=255)
                                             ccACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCca                                                                                                 <  2:952961/150‑1 (MQ=255)
                                                     aCGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACAcccccc                                                                                         >  2:422180/1‑150 (MQ=255)
                                                       gAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCtt                                                                                       >  1:257586/1‑150 (MQ=255)
                                                             gTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCcaca                                                                                               <  2:471371/136‑1 (MQ=255)
                                                             gTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCcaca                                                                                               >  1:471371/1‑136 (MQ=255)
                                                                 accacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTa                                                                             <  2:257586/150‑1 (MQ=255)
                                                                   cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAAc                                                                           >  2:640812/1‑150 (MQ=255)
                                                                   cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAAc                                                                           >  2:436075/1‑150 (MQ=255)
                                                                   cacGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAAc                                                                           <  1:559945/150‑1 (MQ=255)
                                                                      gcggcgTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCAc                                                                        <  1:819241/150‑1 (MQ=255)
                                                                        ggcgTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACga                                                                      <  1:632513/150‑1 (MQ=255)
                                                                               gCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCtt                                                               <  1:436075/150‑1 (MQ=255)
                                                                                 aGCCCATGATCAGGTGGAGGTTACTCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATTGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                             <  1:688751/150‑1 (MQ=255)
                                                                                 aGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTgg                                                             >  2:804044/1‑150 (MQ=255)
                                                                                   cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGat                                                                                                                                       <  1:49946/74‑2 (MQ=255)
                                                                                   cccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAg                                                                                                                                       >  2:49946/1‑74 (MQ=255)
                                                                                    ccATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCa                                                          >  2:909637/1‑150 (MQ=255)
                                                                                     cATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCaa                                                         <  1:886715/150‑1 (MQ=255)
                                                                                          tCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGATGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGgtgt                                                    <  2:1043249/150‑1 (MQ=255)
                                                                                              gTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTg                                                                                                     >  2:434624/1‑97 (MQ=255)
                                                                                              gTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTg                                                                                                     <  1:434624/97‑1 (MQ=255)
                                                                                               tCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAc                                               >  2:1255893/1‑150 (MQ=255)
                                                                                                 gAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                           <  1:450515/88‑1 (MQ=255)
                                                                                                 gAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                           >  2:450515/1‑88 (MQ=255)
                                                                                                  aGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAcccc                                            <  1:422180/150‑1 (MQ=255)
                                                                                                   ggTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCt                                           >  1:894142/1‑150 (MQ=255)
                                                                                                            ggCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGt                                  <  1:804044/150‑1 (MQ=255)
                                                                                                              cccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCg                                <  2:324158/150‑1 (MQ=255)
                                                                                                               ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                           <  2:651721/74‑1 (MQ=255)
                                                                                                               ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTg                                                                                                           >  1:651721/1‑74 (MQ=255)
                                                                                                               ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAg                                                        <  2:136470/125‑1 (MQ=255)
                                                                                                               ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAg                                                        >  1:136470/1‑125 (MQ=255)
                                                                                                               ccACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCgg                               >  1:766943/1‑150 (MQ=255)
                                                                                                                            gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGaa                                                                                >  1:1361485/1‑88 (MQ=255)
                                                                                                                            gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGaa                                                                                <  2:1361485/88‑1 (MQ=255)
                                                                                                                            gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAACACTAGTCGGGACACTTATACCt                  >  2:441282/1‑150 (MQ=255)
                                                                                                                            gACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCt                  <  1:640812/150‑1 (MQ=255)
                                                                                                                                    cGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCaaaaaa          <  2:1351431/150‑1 (MQ=255)
                                                                                                                                         ttGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGaa     >  2:428116/1‑150 (MQ=255)
                                                                                                                                         ttGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGaa     >  1:579836/1‑150 (MQ=255)
                                                                                                                                         taccggCCTCGCACTGCGCGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAc                                               <  1:455469/104‑1 (MQ=255)
                                                                                                                                         taccggCCTCGCACTGCGCGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCAc                                               >  2:455469/5‑108 (MQ=255)
                                                                                                                                            cggCCTCGCACTGCGCGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATa  >  2:356681/2‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTTGTCCTTGATCAGGGCCTGTTGCAGGTTGGCTTCGCCGTTGACCACGTTGACGAAGTCGTACACCACGCGGCGTTCGCAGCCCATGATCAGGTCGAGGTTACGCAGGCCCACGTCGAAGTCGACAATGACGGTCTTGTGGCCGCGCAGTGCGAGGCCGGTACCAATGGCGGCGCTGGTAGTGGTCTTGCCCACACCCCCCTTGCCGGAAGTAACCACGAGAATCTTGGCCAAGGTGTTTCACCCCTAAAAATAGTCGGGACACTTGTCCCTGCAAAAAACTGAAATA  >  NC_002947/1933562‑1933851

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: ‑