Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,803,118 C→T 6.8% A246T (GCA→ACA)  PP_5081 ← type IV pili biogenesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,803,1180CT 6.8% 170.0 / 2.8 59A246T (GCA→ACA) PP_5081type IV pili biogenesis protein
Reads supporting (aligned to +/- strand):  ref base C (35/20);  new base T (2/2);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CTCGTGCCATACCCCCGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCGGCAGGTGCGCC  >  NC_002947/5802969‑5803260
                                                                                                                                                     |                                                                                                                                              
cTCGTGCCATACCCCCGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTgc                                                                                                                                                <  2:831351/150‑1 (MQ=255)
           tacccGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCcgcgcg                                                                                                                                     >  1:948254/3‑150 (MQ=255)
            ccccGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCgcgcgc                                                                                                                                    >  1:382172/1‑150 (MQ=255)
              ccGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCa                                                                                                                                  >  2:923558/1‑150 (MQ=255)
                          gtcgcgCACGCACTGCGCCGTGACGTTGGTCTATATGCATGACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCa                                                                                                                      <  2:185220/148‑1 (MQ=255)
                             gCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCAttt                                                                                                                   >  1:795389/1‑150 (MQ=255)
                                   gCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCTGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCt                                                                                                             >  2:937808/1‑150 (MQ=255)
                                       ttCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTgcgc                                                                                                         >  1:658639/1‑150 (MQ=255)
                                           aCCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCaa                                                                                                      <  1:134191/149‑1 (MQ=255)
                                           aCCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCaa                                                                                                      >  2:134191/1‑149 (MQ=255)
                                            cttGACGTTGGACTATATGCATCACTACCCCATGATGACGCCGTACACGGTCGCCCGGCCGTAGGCGATACAGCTTGCCAGCCGCCAAAAACAGTGCGTCGCGTGCAGCCACGCGCGCAAGGCTGCCTACCATTTGAAGCTGCGCCCAct                                                                                                    >  2:156206/3‑149 (MQ=255)
                                               gACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAAc                                                                                                     <  1:900517/146‑1 (MQ=255)
                                               gACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAAc                                                                                                     >  2:900517/1‑146 (MQ=255)
                                               gACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCa                                                                                                 >  1:454835/1‑150 (MQ=255)
                                                  gTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGc                                                                                              >  1:563404/1‑150 (MQ=255)
                                                       tCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGAcc                                                                                         >  1:560866/1‑150 (MQ=255)
                                                            atGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCCTGCCGCCCCGCCCCCCAGGCTGCCGACCATTTCCACCTTCTCCAACGGCGCGCCGACAAGccc                                                                                    >  2:197907/1‑150 (MQ=255)
                                                              gCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCgg                                                                                  >  2:14175/1‑150 (MQ=255)
                                                              gCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCgg                                                                                  >  1:1061445/1‑150 (MQ=255)
                                                                aTCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGgcg                                                                                >  2:1018507/1‑150 (MQ=255)
                                                                aTCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGgcg                                                                                >  1:902738/1‑150 (MQ=255)
                                                                     tACGCCCTGATCACGCCCTACGCGATCGCCCGGCCGTAGGCGATAGAGCTTGCCAGCCGCCATCAGAAGTGCATCGCGTGCCGCCCGGCGCGCCAGGCTGCCGCCCATTTCCACCTGCGCCACGCGCCGTCCGACCAACCCGGCGCTGGc                                                                           >  2:482737/1‑150 (MQ=255)
                                                                        ccccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCAc                                                                        >  1:382370/1‑150 (MQ=255)
                                                                         cccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCTATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAAc                                                                                                     >  1:630306/1‑120 (MQ=255)
                                                                         cccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAAc                                                                                                     <  2:630306/120‑1 (MQ=255)
                                                                         cccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACCCCGAAGAGCCAGGCGCTGGCCACt                                                                       >  2:821719/1‑150 (MQ=255)
                                                                          ccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCGGACCATTTCCATCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTc                                                                      <  2:948254/150‑1 (MQ=255)
                                                                          ccTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTc                                                                      >  1:307813/1‑150 (MQ=255)
                                                                            tGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTcgc                                                                    <  2:1061445/150‑1 (MQ=255)
                                                                              aTCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCgcgc                                                                  <  1:144159/150‑1 (MQ=255)
                                                                                  cGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGACCGGCGCTGGCCACTCGCGCCGCa                                                              >  2:893753/1‑150 (MQ=255)
                                                                                   gCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAg                                                             >  2:435848/1‑150 (MQ=255)
                                                                                    cccTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGc                                                            >  1:44114/1‑150 (MQ=255)
                                                                                       tACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCggg                                                         >  1:1189678/1‑150 (MQ=255)
                                                                                        aCCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGgc                                                        >  1:814331/1‑150 (MQ=255)
                                                                                           cGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGTCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCa                                                     >  1:873113/1‑150 (MQ=255)
                                                                                           cGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGTCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCa                                                     >  1:872966/1‑150 (MQ=255)
                                                                                                   cGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTc                                             <  2:563404/150‑1 (MQ=255)
                                                                                                     ggCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTc                                             >  1:504860/1‑148 (MQ=255)
                                                                                                     ggCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTc                                             <  2:504860/148‑1 (MQ=255)
                                                                                                     ggCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCa                                           <  1:768093/150‑1 (MQ=255)
                                                                                                      gCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGACCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCAc                                          <  1:893753/150‑1 (MQ=255)
                                                                                                           cgcgATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAg                                     <  2:307813/150‑1 (MQ=255)
                                                                                                                tAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCa                                <  2:1189678/150‑1 (MQ=255)
                                                                                                                       tGCCAGCCGCCATCAGCAGTGCCTCGCGTGTCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTg                         <  2:873113/150‑1 (MQ=255)
                                                                                                                       tGCCAGCCGCCATCAGCAGTGCCTCGCGTGTCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTg                         <  2:872966/150‑1 (MQ=255)
                                                                                                                        gCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGc                        <  2:607806/150‑1 (MQ=255)
                                                                                                                         ccAGCCGCAATAAGAAGTGCCTCACGTGCCGCTCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCg                       <  2:560866/150‑1 (MQ=255)
                                                                                                                          cAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCgg                      <  1:923558/150‑1 (MQ=255)
                                                                                                                            gccgccATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGcc                    <  1:213417/150‑1 (MQ=255)
                                                                                                                             ccgccATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCa                                      >  1:25360/1‑131 (MQ=255)
                                                                                                                             ccgccATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCa                                      <  2:25360/131‑1 (MQ=255)
                                                                                                                                ccATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCGGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTg                >  1:800964/1‑150 (MQ=255)
                                                                                                                                  aTCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCATCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGcc              >  1:379034/1‑150 (MQ=255)
                                                                                                                                   tCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCg             >  2:78063/1‑150 (MQ=255)
                                                                                                                                   tCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCg             <  2:454835/150‑1 (MQ=255)
                                                                                                                                       cagTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCCACCACAGCCATTGCGGCCCTTGCCGGCAg         <  1:937808/150‑1 (MQ=255)
                                                                                                                                            ccTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCGGCAGGTgcg    >  1:923655/1‑150 (MQ=255)
                                                                                                                                              tCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCGGCAGGTGCGcc  >  2:793929/1‑150 (MQ=255)
                                                                                                                                              tCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGACGGCAGGTGCGcc  >  1:666821/1‑150 (MQ=255)
                                                                                                                                                    gccgccCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGc                              >  2:600585/1‑116 (MQ=255)
                                                                                                                                                    gccgccCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGc                              <  1:600585/116‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CTCGTGCCATACCCCCGCCATGAACAGCCGCTCACGCACTTCGACCTGACGTTGGTCTATATGCATCACTACCCCCTGATCACGCCCTACCCGGTCGCCCGGGCGTACGCGATAGAGCTTGCCAGCCGCCATCAGCAGTGCCTCGCGTGCCGCCCCGCGCGCCAGGCTGCCGACCATTTCCAGCTGCGCCAACGGCACGCCGACCAGCCCGGCGCTGGCCACTCGCGCCGCAGCGGGCGCAAACGGGTCCACTTCAGGCGCAGCCATTGCGGCCCTTGCCGGCAGGTGCGCC  >  NC_002947/5802969‑5803260

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: