Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,286,411 A→G 6.8% V33A (GTC→GCC)  PP_3757 ← chemotaxis protein CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,286,4110AG 6.8% 205.9 / 3.9 73V33A (GTC→GCC) PP_3757chemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base G (3/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCTTCAGGTCGGGCCGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTA  >  NC_002947/4286262‑4286555
                                                                                                                                                     |                                                                                                                                                
ccTTCAGGTCGGGCCGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGa                                                                                                                                                  >  2:1651135/1‑150 (MQ=255)
     aGGTCGGGCCGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTgg                                                                                                                                             >  2:1511365/1‑150 (MQ=255)
     aGGTCGGGCCGCAACTTCACCGCGGGCTCGAAAAACAACACGCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCTGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTgg                                                                                                                                             <  2:1096183/150‑1 (MQ=255)
        tCGGGCCGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTAc                                                                                                                                          <  1:1330347/150‑1 (MQ=255)
              cGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCa                                                                                                                                    <  2:30692/150‑1 (MQ=255)
               gCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCAc                                                                                                                                   <  2:1428475/150‑1 (MQ=255)
               gCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCAc                                                                                                                                   >  1:931566/1‑150 (MQ=255)
                      ctccGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAg                                                                                                                            <  1:661102/148‑1 (MQ=255)
                                        atCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCAt                                                                                                          >  1:1210203/1‑150 (MQ=255)
                                         cccACTCCCGACTCCCCCCCCGGCAAACGGAAGCCGATCACCATGATTTCCAGGTGCGGTTTGCTGGCCAAGGTGGCAAACGCCTGCTCTCCATCCTCAGCCACCAGGACGTGGTACCCCTCTCCCGCCAGGTAATCGCGCAAAATCATg                                                                                                         <  1:423883/149‑1 (MQ=255)
                                          ccACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGc                                                                                                        <  1:1588543/150‑1 (MQ=255)
                                                ccGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCg                                                                                                  >  1:263623/1‑150 (MQ=255)
                                                   agaTCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCg                                                                                               <  1:1029392/150‑1 (MQ=255)
                                                   agaTCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCg                                                                                               >  1:57523/1‑150 (MQ=255)
                                                       cccccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTc                                                                                           >  2:707575/1‑150 (MQ=255)
                                                         cccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATg                                                                                                         >  1:1219919/1‑134 (MQ=255)
                                                         cccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATg                                                                                                         <  2:1219919/134‑1 (MQ=255)
                                                         cccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCAt                                                                                         >  2:378226/1‑150 (MQ=255)
                                                         cccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCAt                                                                                         >  2:216740/1‑150 (MQ=255)
                                                          ccccGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATc                                                                                        >  2:1710037/1‑150 (MQ=255)
                                                               tcAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTcga                                                                                   >  1:1715457/2‑150 (MQ=255)
                                                                 aaaCGGAAGTCCGTCACCATCAGGTCCAGGTTCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTcgacg                                                                                 >  2:1269237/1‑150 (MQ=255)
                                                                    cGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGAcc                                                                              <  2:1862132/150‑1 (MQ=255)
                                                                      gAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCg                                                                                               >  2:1775453/1‑131 (MQ=255)
                                                                      gAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGCACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCg                                                                                               <  1:1775453/131‑1 (MQ=255)
                                                                          tCCGTCACCGTCAGGTAAAAGTGCGGTTTGCTGGGCAGGATGGCAAACGCCTGCTCGCCATTCGCATCCACCAGGACCTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCAAATCGCCGGCTGATCCTCGACGACCAGAAcc                                                                        >  2:21414/1‑150 (MQ=255)
                                                                          tCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGGCGTGGTACCCCTCACCCGCCa                                                                                                                             >  2:1829409/1‑97 (MQ=255)
                                                                          tCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGGCGTGGTACCCCTCACCCGCCa                                                                                                                             >  2:1829421/1‑97 (MQ=255)
                                                                          tCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGGCGTGGTACCCCTCACCCGCCa                                                                                                                             <  1:1829409/97‑1 (MQ=255)
                                                                          tCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGGCGTGGTACCCCTCACCCGCCa                                                                                                                             <  1:1829421/97‑1 (MQ=255)
                                                                             gTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACcaca                                                                     >  2:1226619/1‑150 (MQ=255)
                                                                                aCCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCt                                                                  <  1:707575/150‑1 (MQ=255)
                                                                                  cATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTtg                                                                <  1:1651135/150‑1 (MQ=255)
                                                                                   aTCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTtgt                                                               >  2:523133/1‑150 (MQ=255)
                                                                                   aTCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTtgt                                                               <  1:1438599/150‑1 (MQ=255)
                                                                                   aTCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTtgt                                                               <  1:1430578/150‑1 (MQ=255)
                                                                                   aTCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTtgt                                                               <  1:216740/150‑1 (MQ=255)
                                                                                        gTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCCGGTAATCGCGCAAAATCATGCGAATCGCCTGGCCATCCTCGACGACCAGAACCACATCTTGTGCATc                                                          >  1:99276/1‑150 (MQ=255)
                                                                                           cAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGGCGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCAAATCAAGTGCATCTTc                                                       >  2:986861/1‑150 (MQ=255)
                                                                                            aGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCgacgac                                                                               >  2:728920/1‑125 (MQ=255)
                                                                                            aGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCgacgac                                                                               <  1:728920/125‑1 (MQ=255)
                                                                                                 cGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCAt                                                 <  1:523133/150‑1 (MQ=255)
                                                                                                        cTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCGTCTTCACTCATGTTAACg                                          >  2:1319834/1‑150 (MQ=255)
                                                                                                         tGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGc                                         >  1:1091589/1‑150 (MQ=255)
                                                                                                         tGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGAAGGCCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGc                                         >  1:311071/1‑150 (MQ=255)
                                                                                                           gCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCt                                       >  2:1721206/1‑150 (MQ=255)
                                                                                                            ccAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCtt                                      <  1:1511365/150‑1 (MQ=255)
                                                                                                             cAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTg                                     >  1:1357009/1‑150 (MQ=255)
                                                                                                                 aTGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGa                                    <  1:370220/147‑1 (MQ=255)
                                                                                                                 aTGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGa                                    >  2:370220/1‑147 (MQ=255)
                                                                                                                   ggCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCGTCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATAc                               <  2:263623/150‑1 (MQ=255)
                                                                                                                      aaaCGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTa                            >  1:1043613/1‑150 (MQ=255)
                                                                                                                      aaaCGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTa                            >  2:1734977/1‑150 (MQ=255)
                                                                                                                      aaaCGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTa                            <  1:1269237/150‑1 (MQ=255)
                                                                                                                          gCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCTAGGTAATCGCGCAAAATCATGCGAATCGCCGGGTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCGTGAAATACTTACAGt                        >  2:468328/1‑150 (MQ=255)
                                                                                                                           ccTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGtt                       <  2:57523/150‑1 (MQ=255)
                                                                                                                           ccTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACAAATCTTAACGCCTTGAAATACTTACAGtt                       <  1:731278/150‑1 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGaa                                                                          <  1:659368/97‑1 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGaa                                                                          >  2:659368/1‑97 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAAc                                                                         >  1:602601/1‑98 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAAc                                                                         <  2:602601/98‑1 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGTATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTgg                     >  1:1622368/1‑150 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTgg                     >  1:1268162/1‑150 (MQ=255)
                                                                                                                             tGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTgg                     <  1:1721206/150‑1 (MQ=255)
                                                                                                                                 cGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGTCGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCAt                 <  2:931566/150‑1 (MQ=255)
                                                                                                                                  gCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATc                >  1:373139/1‑150 (MQ=255)
                                                                                                                                   tcgtCCTCAGCCTCCGGGACGTGGGTCCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCa               <  2:1357009/147‑1 (MQ=255)
                                                                                                                                   ccATCCTCAGCCACCAGGACGTGGTACCCGTCACCGGCGTGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCGTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCa               <  2:311071/150‑1 (MQ=255)
                                                                                                                                         tCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGc                                                             >  1:980235/1‑98 (MQ=255)
                                                                                                                                         tCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGc                                                             <  2:980235/98‑1 (MQ=255)
                                                                                                                                         tCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGAcccc         >  1:201712/1‑150 (MQ=255)
                                                                                                                                          cAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCa        <  1:561015/150‑1 (MQ=255)
                                                                                                                                            gCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAgg      <  2:1715457/150‑1 (MQ=255)
                                                                                                                                                ccaGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTa  >  2:475440/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCTTCAGGTCGGGCCGCAACTTCACCGCGGGCTCGGCAATATCCACTCCCGAGATCCCCCCCGGCAAACGGAAGTCCGTCACCATCAGGTCCAGGTGCGGTTTGCTGGCCAGGATGGCAAACGCCTGCTCGCCATCCTCAGCCACCAGGACGTGGTACCCCTCACCCGCCAGGTAATCGCGCAAAATCATGCGAATCGCCGGTTCATCCTCGACGACCAGAACCACATCTTGTGCATCTTCACTCATGTTAACGCCTTGAAATACTTACAGTTGGCCATCAGACCCCAGGCTTA  >  NC_002947/4286262‑4286555

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