Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,548 T→A 10.7% S365C (AGC→TGC)  PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5480TA10.7% 225.9 / 18.0 85S365C (AGC→TGC) PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base A (4/5);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCA  >  NC_002947/1261401‑1261694
                                                                                                                                                    |                                                                                                                                                  
ggCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTTCCTTCGCAATTGGCCAACAGGTA‑GCCGCTGc                                                                                                                                                  <  1:1483372/150‑1 (MQ=255)
   ttCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCgac                                                                                                                                               >  2:1318468/1‑148 (MQ=255)
   ttCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGc                                                                                                                                               >  1:1567288/1‑150 (MQ=255)
   ttCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGc                                                                                                                                               >  1:1707640/1‑150 (MQ=255)
    tCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCa                                                                                                                                              <  1:1649548/150‑1 (MQ=255)
       aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATAc                                                                                                                                           <  1:1164364/150‑1 (MQ=255)
       aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATAc                                                                                                                                           <  1:356358/150‑1 (MQ=255)
       aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATAc                                                                                                                                           >  1:1926106/1‑150 (MQ=255)
        gCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCTGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGAATAcc                                                                                                                                          >  2:833505/1‑150 (MQ=255)
         ccGCGCTTGGCAAACCAGTCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCt                                                                                                                                         <  1:1882906/150‑1 (MQ=255)
         ccGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCCACAGGTA‑GCCGCTGCGGCATACCt                                                                                                                                         <  2:169451/150‑1 (MQ=255)
             gcTTGGCAAACCAGGCATTGGCGAGGTCGGTCGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAGCGGC‑TACCTGttg                                                                                                                                     <  1:9976/150‑2 (MQ=255)
             gcTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAGCGGC‑TACCTGttg                                                                                                                                     <  1:1393751/150‑2 (MQ=255)
                  gCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCgg                                                                                                                                >  1:1037049/1‑150 (MQ=255)
                          ggCATTGGCGAGGTCGGTAACTTCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCCGTGCTTCTACAGACCACTTGACGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATgg                                                                                                                        <  1:1038542/150‑1 (MQ=255)
                           gCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GACGCTGCGGCTTACCTGGTCGCCGGCCACATggg                                                                                                                       >  2:144903/1‑150 (MQ=255)
                           gCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGGGCTTCTCCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATggg                                                                                                                       <  2:1538748/150‑1 (MQ=255)
                             aTTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATGAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCg                                                                                                                     >  2:1038328/1‑150 (MQ=255)
                               gggCGAGGTCGGTAGGCGCACGCATGAGGGGCGGCGGTCCGTGGACCAGCCCGCCGTTACAGAGCACTTGCCGGCTGCAATGGTCTTACCGGATCAATTGGCCAACAGGGA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                   <  2:990513/149‑1 (MQ=255)
                                 gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGa                                                                                                                 <  1:457539/150‑1 (MQ=255)
                                 gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGa                                                                                                                 >  2:520558/1‑150 (MQ=255)
                                   gAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGAcc                                                                                                               >  1:583901/1‑150 (MQ=255)
                                           ttgcagcaGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCCCATGGGCGATGACCATGCCAGc                                                                                                       >  2:903137/3‑150 (MQ=255)
                                             gcagcaGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCgg                                                                                                     >  2:1396560/1‑150 (MQ=255)
                                               ggcagGCGTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   <  2:1239369/149‑1 (MQ=255)
                                               agcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   <  2:1707640/150‑1 (MQ=255)
                                               agcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   >  1:186732/1‑150 (MQ=255)
                                                gcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTgg                                                                                                  <  1:1723414/150‑1 (MQ=255)
                                                gcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTgg                                                                                                  <  1:520558/150‑1 (MQ=255)
                                                 cagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGc                                                                                                 >  2:1280535/1‑150 (MQ=255)
                                                    gCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCAt                                                                                              <  1:266531/150‑1 (MQ=255)
                                                     cATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATg                                                                                             >  2:1602428/1‑150 (MQ=255)
                                                            gCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGa                                                                                      <  1:1038328/150‑1 (MQ=255)
                                                            gCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGACGACCATGCCAGCGGTGGCCATGCGCCGGa                                                                                      <  2:1369943/150‑1 (MQ=255)
                                                               ggCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCCGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGc                                                                                   <  1:256642/150‑1 (MQ=255)
                                                                gCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATgcgc                                                                                          <  1:1485135/142‑1 (MQ=255)
                                                                gCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATgcgc                                                                                          >  2:1485135/1‑142 (MQ=255)
                                                                gCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGACATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCt                                                                                  <  1:1318468/150‑1 (MQ=255)
                                                                    tCCGTCGAAGGCGGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTCCCGGGTCAATTGGCCAGCAGGTA‑GCCGGTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCTAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                              <  2:244857/150‑1 (MQ=255)
                                                                       ttCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTATGCCGCAGCGGC‑TACCTGttgg                                                                                                                                    <  1:651212/93‑3 (MQ=25)
                                                                       ttCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTATGCCGCAGCGGC‑TACCTGttgg                                                                                                                                    >  2:651212/1‑91 (MQ=25)
                                                                         cGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATa                                                                         <  2:1567288/150‑1 (MQ=255)
                                                                          gAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAc                                                                        >  2:1457625/1‑150 (MQ=255)
                                                                              gCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                    >  2:70734/1‑150 (MQ=255)
                                                                                     ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGGTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCGTACAGCAGCCCCgg                                                             >  2:590968/1‑150 (MQ=255)
                                                                                     ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACGGCCGCCCCgg                                                             >  2:687272/1‑150 (MQ=255)
                                                                                     ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                             >  2:875645/1‑150 (MQ=255)
                                                                                        tACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctg                                                          >  1:304666/1‑150 (MQ=255)
                                                                                         aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                         >  2:373729/1‑150 (MQ=255)
                                                                                         aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                         <  2:186732/150‑1 (MQ=255)
                                                                                         aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCAATTTTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                         <  1:33931/150‑1 (MQ=255)
                                                                                              ccACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                    <  1:1715340/150‑1 (MQ=255)
                                                                                               cACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTg                                                   >  1:1280870/1‑150 (MQ=255)
                                                                                                  ttGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTgcgc                                                >  2:257924/1‑150 (MQ=255)
                                                                                                   tGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCa                                               <  1:373729/150‑1 (MQ=255)
                                                                                                    gCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAg                                              <  2:1280870/150‑1 (MQ=255)
                                                                                                        gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                          <  2:814593/150‑1 (MQ=255)
                                                                                                        gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                          <  1:409066/150‑1 (MQ=255)
                                                                                                        gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCGGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                          >  2:673886/1‑150 (MQ=255)
                                                                                                            cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      >  1:1626938/1‑150 (MQ=255)
                                                                                                            cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      >  1:498570/1‑150 (MQ=255)
                                                                                                            cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGAACATGCCAGCGGTGGCAACGCGCCGGACGCGACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      <  1:1180372/150‑1 (MQ=255)
                                                                                                             aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                   >  1:142055/1‑72 (MQ=255)
                                                                                                             aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                   <  2:142055/72‑1 (MQ=255)
                                                                                                             aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCAAGAGCGTCGTGTCCGGCTGGCAGTTGCGCAGGGGGGTGAc                                     >  2:152571/1‑150 (MQ=255)
                                                                                                                ggTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCg                                  <  2:304666/150‑1 (MQ=255)
                                                                                                                  tttgACCGGATGAATTGGGCAGCGGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGCCCTTGCCAGCGGTGGCCATGCGCCAGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGt                                <  2:1004376/150‑1 (MQ=255)
                                                                                                                     tACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCAc                             <  1:1236243/150‑1 (MQ=255)
                                                                                                                       ccGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTg                           >  2:357617/1‑150 (MQ=255)
                                                                                                                        cGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGa                          <  2:583901/150‑1 (MQ=255)
                                                                                                                        cGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGa                          <  1:1457625/150‑1 (MQ=255)
                                                                                                                          taTCAATTGGCCAACAGGTA‑GCCGCTGCGCCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                        <  1:777866/149‑1 (MQ=255)
                                                                                                                           aTCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTg                       >  1:711389/1‑150 (MQ=255)
                                                                                                                               aTTGGCCAACAGGTA‑GCCGCTGCGGTATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCAATGAGTGGGTc                   >  1:546591/1‑150 (MQ=255)
                                                                                                                                 tGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGa                 <  1:379773/150‑1 (MQ=255)
                                                                                                                                    ccAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGa              >  1:991297/1‑150 (MQ=255)
                                                                                                                                     gaaCAGGTATGCCGCAGCGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGAc             >  2:805609/2‑150 (MQ=255)
                                                                                                                                     cAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGCGGGTCGATGAc             >  2:1914488/1‑150 (MQ=255)
                                                                                                                                       aCAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACtt           <  2:1626938/150‑1 (MQ=255)
                                                                                                                                        cAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   <  2:1677621/61‑1 (MQ=17)
                                                                                                                                        cAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   >  1:1677621/1‑61 (MQ=17)
                                                                                                                                        cAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACt                            <  1:970133/132‑1 (MQ=255)
                                                                                                                                        cAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACt                            >  2:970133/1‑132 (MQ=255)
                                                                                                                                          ggTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCgg        <  1:257924/150‑1 (MQ=255)
                                                                                                                                               gCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGAt    <  2:1898058/150‑1 (MQ=255)
                                                                                                                                                 cGCNGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTTAGTGGGTCGATGACTTCGGCGATCa  >  2:1187220/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCA  >  NC_002947/1261401‑1261694

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