Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,110 T→A 11.1% intergenic (+32/‑82) PP_5353 → / → PP_5354 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,1100TA11.1% 206.2 / 15.2 72intergenic (+32/‑82)PP_5353/PP_5354hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base A (3/5);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTC  >  NC_002947/6103965‑6104255
                                                                                                                                                 |                                                                                                                                                 
tCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTAGTAGCCGGTATCCGTCAAAACAGATAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGccc                                                                                                                                               >  2:785943/1‑150 (MQ=255)
  gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                                                                                                             >  2:1201843/1‑150 (MQ=255)
   aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                                                                                                            >  1:173855/1‑150 (MQ=255)
   aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCggg                                                                                                                                            <  1:436766/150‑1 (MQ=255)
      aaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATTGTTCGTTCCTAAGTCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCGGGGtt                                                                                                                                         <  1:774624/150‑1 (MQ=255)
        ttCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCTAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGtttt                                                                                                                                       >  2:312555/1‑150 (MQ=255)
         tCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGttttt                                                                                                                                      <  1:1201843/150‑1 (MQ=255)
           gcgcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACAAGCCTTGGTTGATAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCGGGGttctgtg                                                                                                                                    <  1:1432848/150‑6 (MQ=255)
             gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                                                                                                                  >  1:1424922/1‑150 (MQ=255)
              cTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCACCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttt                                                                                                                                 >  1:1622306/1‑150 (MQ=255)
               ttCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctttt                                                                                                                                <  1:1145761/150‑1 (MQ=255)
               ttCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCACCCAAGACGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctttt                                                                                                                                <  2:1622306/150‑1 (MQ=255)
                    ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaa                                                                                                                           >  2:743776/1‑150 (MQ=255)
                       gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGt                                                                                                                        <  2:1146241/150‑1 (MQ=255)
                         gCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTgg                                                                                                                      <  1:1079972/150‑1 (MQ=255)
                                 ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                                                                                                           >  1:1027901/1‑121 (MQ=255)
                                 ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                                                                                                           <  2:1027901/121‑1 (MQ=255)
                                     gCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                                                                                                             <  1:1616540/147‑1 (MQ=255)
                                     gCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                                                                                                             >  2:1616540/1‑147 (MQ=255)
                                       gAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGtata                                                                                                        >  1:799404/1‑150 (MQ=255)
                                               gCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAg                                                                                                <  2:1424922/150‑1 (MQ=255)
                                                   aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                            >  1:761440/1‑150 (MQ=255)
                                                   aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                            <  2:427404/150‑1 (MQ=255)
                                                      gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACATCCGTCCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTc                                                                                         <  1:115339/150‑1 (MQ=255)
                                                       ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCa                                                                                        >  2:839276/1‑150 (MQ=255)
                                                       ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCggg                                                                                                                                            >  2:1162196/1‑98 (MQ=255)
                                                       ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCggg                                                                                                                                            <  1:1162196/98‑1 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                                                                                                            >  2:708276/1‑123 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                                                                                                            <  1:708276/123‑1 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGta                                                                                                          <  2:971399/125‑1 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGta                                                                                                          >  1:971399/1‑125 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCa                                                                                        <  2:1523566/143‑1 (MQ=255)
                                                              ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCa                                                                                        >  1:1523566/1‑143 (MQ=255)
                                                                   gCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCagag                                                                            <  2:260343/150‑1 (MQ=255)
                                                                       gaACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTGTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAg                                                                        <  1:670463/150‑1 (MQ=255)
                                                                        aaCCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGt                                                                       <  1:743776/150‑1 (MQ=255)
                                                                        aaCCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGt                                                                       >  1:182340/1‑150 (MQ=255)
                                                                           cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                                                                                                           <  2:1558433/79‑1 (MQ=255)
                                                                           cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                                                                                                           >  1:1558433/1‑79 (MQ=255)
                                                                           cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTCaa                                                                    <  1:1642618/150‑1 (MQ=255)
                                                                           cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttct                                                                                                                                       >  1:1138927/1‑81 (MQ=38)
                                                                           cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttct                                                                                                                                       <  2:1138927/83‑3 (MQ=38)
                                                                               cTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGt                                                                <  2:799404/150‑1 (MQ=255)
                                                                                ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                            >  2:911405/1‑121 (MQ=255)
                                                                                ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                            <  1:911405/121‑1 (MQ=255)
                                                                                 tGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTgg                                                              <  1:839276/150‑1 (MQ=255)
                                                                                            cGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATATAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATc                                                   <  2:1568034/150‑1 (MQ=255)
                                                                                              gTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGa                                                 >  1:580335/1‑150 (MQ=255)
                                                                                               tATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaa                                                <  2:173855/150‑1 (MQ=255)
                                                                                                         aaCAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaac                                               <  1:1513565/141‑1 (MQ=255)
                                                                                                         aaCAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaac                                               >  2:1513565/1‑141 (MQ=255)
                                                                                                                aaTCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTGTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTACGTGAAGCGaa                               >  2:822251/1‑150 (MQ=255)
                                                                                                                 aTCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAAt                              >  2:1081946/1‑150 (MQ=255)
                                                                                                                     cTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGc                          <  2:952686/150‑1 (MQ=255)
                                                                                                                     cTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGc                          <  1:1081946/150‑1 (MQ=255)
                                                                                                                      tGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGcc                         >  2:925523/1‑150 (MQ=255)
                                                                                                                       gCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCg                        <  2:580335/150‑1 (MQ=255)
                                                                                                                             tAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                                                                                                 >  1:914247/1‑55 (MQ=255)
                                                                                                                             tAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                                                                                                 <  2:914247/55‑1 (MQ=255)
                                                                                                                             tAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcgc                                                                                                  >  2:1479601/1‑70 (MQ=255)
                                                                                                                             tAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcgc                                                                                                  <  1:1479601/70‑1 (MQ=255)
                                                                                                                               aaCACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTAGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAAGCGTGGAAGAAATGATCGAAAAACTTCATAAAGCGAATGAGCCGGTaaagat                >  2:14744/1‑145 (MQ=255)
                                                                                                                                 gcaaaaaCCCCGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGc              >  1:1399549/6‑150 (MQ=255)
                                                                                                                                    aGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGaa                                                                         >  1:1044939/1‑88 (MQ=255)
                                                                                                                                    aGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGaa                                                                         <  2:1044939/88‑1 (MQ=255)
                                                                                                                                    aGAACCCCGGCGTTGCCCGGGGTTCTTGCTTCTGCTAAAGTTTGCTACGAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGcccc           <  2:761440/150‑1 (MQ=255)
                                                                                                                                     gAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCt          <  1:312555/150‑1 (MQ=255)
                                                                                                                                      aaCCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTg         >  2:996724/1‑150 (MQ=255)
                                                                                                                                          ccGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGa                                                     <  2:581692/102‑1 (MQ=255)
                                                                                                                                          ccGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGa                                                     >  1:581692/1‑102 (MQ=255)
                                                                                                                                             ggCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTc  <  2:1561441/150‑1 (MQ=255)
                                                                                                                                             ggCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTc  <  1:996724/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTC  >  NC_002947/6103965‑6104255

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