Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,542:1 +T 10.3% coding (1099/1116 nt) PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5421.T10.3% 318.0 / 31.2 87coding (1099/1116 nt)PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base . (37/41);  new base T (4/5);  total (41/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 = 7.94e-01

GAAGGGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGG  >  NC_002947/1261397‑1261688
                                                                                                                                                  |                                                                                                                                                  
gAAGGGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCg                                                                                                                                                >  1:782969/1‑150 (MQ=255)
gAAGGGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCg                                                                                                                                                >  1:1465077/1‑150 (MQ=255)
   gggCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGcga                                                                                                                                             <  2:1465077/150‑3 (MQ=255)
    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:1926106/1‑150 (MQ=255)
           aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATAc                                                                                                                                     <  1:356358/150‑1 (MQ=255)
           aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATAc                                                                                                                                     <  1:1164364/150‑1 (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                                                                                    >  2:1485135/1‑142 (MQ=255)
                                                                    gCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATgcgc                                                                                    <  1:1485135/142‑1 (MQ=255)
                                                                    gCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGACATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCt                                                                            <  1:1318468/150‑1 (MQ=255)
                                                                        tCCGTCGAAGGCGGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTCCCGGGTCAATTGGCCAGCAGGTA‑GCCGGTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCTAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                        <  2:244857/150‑1 (MQ=255)
                                                                           ttCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTATGCCGCAGCGGC‑TACCTGttgg                                                                                                                              >  2:651212/1‑91 (MQ=25)
                                                                           ttCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTATGCCGCAGCGGC‑TACCTGttgg                                                                                                                              <  1:651212/93‑3 (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:186732/150‑1 (MQ=255)
                                                                                             aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                   >  2:373729/1‑150 (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                                    <  1:409066/150‑1 (MQ=255)
                                                                                                            gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                    <  2:814593/150‑1 (MQ=255)
                                                                                                            gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCGGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                    >  2:673886/1‑150 (MQ=255)
                                                                                                                cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                >  1:498570/1‑150 (MQ=255)
                                                                                                                cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                >  1:1626938/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                    <  1:1457625/150‑1 (MQ=255)
                                                                                                                            cGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGa                    <  2:583901/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                                                                                             >  1:1677621/1‑61 (MQ=17)
                                                                                                                                            cAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                             <  2:1677621/61‑1 (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)
                                                                                                                                                  |                                                                                                                                                  
GAAGGGCTTCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGG  >  NC_002947/1261397‑1261688

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