Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,554 Δ1 bp 11.1% coding (1087/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,5540A.11.1% 227.7 / 23.8 81coding (1087/1116 nt)PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base . (4/5);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

TCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTGGCCGG  >  NC_002947/1261405‑1261703
                                                                                                                                                      |                                                                                                                                                     
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                                                                                                                          >  2:520558/1‑150 (MQ=255)
                             gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGa                                                                                                                          <  1:457539/150‑1 (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                                                                                                            >  1:186732/1‑150 (MQ=255)
                                           agcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                            <  2:1707640/150‑1 (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                                                                                                                                             >  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                                                   <  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:498570/1‑150 (MQ=255)
                                                                                                        cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                               >  1:1626938/1‑150 (MQ=255)
                                                                                                        cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGAACATGCCAGCGGTGGCAACGCGCCGGACGCGACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                               <  1:1180372/150‑1 (MQ=255)
                                                                                                         aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                            <  2:142055/72‑1 (MQ=255)
                                                                                                         aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                            >  1:142055/1‑72 (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                                                                                                            <  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)
                                                                                                                                                      aTACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTggccgg  >  2:10526/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCCAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTGGCCGG  >  NC_002947/1261405‑1261703

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