Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 339,228 C→G 6.8% M12I (ATG→ATC)  PP_0281 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947339,2280CG 6.8% 222.3 / 7.4 74M12I (ATG→ATC) PP_0281amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (28/41);  new base G (3/2);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

ACGCCGCGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGC  >  NC_002947/339079‑339374
                                                                                                                                                     |                                                                                                                                                  
aCGCCGCGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGTCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCcatc                                                                                                                                                    >  1:128042/1‑150 (MQ=255)
 cgccgcGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCcatca                                                                                                                                                   <  1:304183/150‑1 (MQ=255)
     gcgAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGc                                                                                                                                               <  1:235123/150‑1 (MQ=255)
     gcgAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGc                                                                                                                                               >  2:1333960/1‑150 (MQ=255)
        aaCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATGATGGCCgg                                                                                                                                            <  1:1310184/150‑1 (MQ=255)
         aCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGa                                                                                                                                           <  1:7465/150‑1 (MQ=255)
         aCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGa                                                                                                                                           >  2:1298158/1‑150 (MQ=255)
         aCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGa                                                                                                                                           <  2:656295/150‑1 (MQ=255)
         aCCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGa                                                                                                                                           >  2:519577/1‑150 (MQ=255)
               gTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAAt                                                                                                                                     <  1:1333960/150‑1 (MQ=255)
                  gAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCg                                                                                                                                  >  1:357415/1‑150 (MQ=255)
                   aGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGt                                                                                                                                 <  1:1447282/150‑1 (MQ=255)
                    gTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTg                                                                                                                                <  2:201852/150‑1 (MQ=255)
                     tAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGc                                                                                                                               <  2:1814675/150‑1 (MQ=255)
                     tAGAAGCCGCCAAGCCATTGCAGCGGCTGGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGc                                                                                                                               >  1:1287403/1‑150 (MQ=255)
                         aGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGt                                                                                                                           >  1:373553/1‑150 (MQ=255)
                         aGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGGCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGt                                                                                                                           <  2:724366/150‑1 (MQ=255)
                          gccgccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTc                                                                                                                          >  2:179378/1‑150 (MQ=255)
                           ccgccAAGCCATTGCAGCGGCTTGAACGGGGAAGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCg                                                                                                                         >  1:883818/1‑150 (MQ=255)
                              ccAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGTCCGGACCGAATCCGTGCAGGTCGata                                                                                                                      <  1:1298158/150‑1 (MQ=255)
                                 aGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTc                                                                                                                   <  2:105699/150‑1 (MQ=255)
                                 aGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTc                                                                                                                   <  2:1284254/150‑1 (MQ=255)
                                     aTTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATgg                                                                                                               >  1:153555/1‑150 (MQ=255)
                                        gCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCa                                                                                                            <  2:991588/150‑1 (MQ=255)
                                         cAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCACATTACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAg                                                                                                           <  1:1645537/150‑1 (MQ=255)
                                                ttGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGACATTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAttt                                                                                                    >  1:426906/1‑150 (MQ=255)
                                                ttGAACGGGGAGGTCTTGGCCAGGCTGCGTAGGAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAttt                                                                                                    >  1:1483974/1‑150 (MQ=255)
                                                   aaCGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAg                                                                                                 <  2:153555/150‑1 (MQ=255)
                                                    aCGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAgg                                                                                                <  1:1474058/150‑1 (MQ=255)
                                                    aCGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAgg                                                                                                <  1:760047/150‑1 (MQ=255)
                                                     cGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCGGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCcatca                                                                                                                                                   >  2:1037068/1‑98 (MQ=255)
                                                     cGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCGGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCcatca                                                                                                                                                   <  1:1037068/98‑1 (MQ=255)
                                                      ggggAAGCCTCGGCCAGGGCGCCGAGCAGGGCCAGGACCAGCCCCAGCAGCAGGGCGCAAAGCGCCAGTTTTGCGGTCAGCAGGGTGCCGGCCATCATGGTCGGACCGAATCCGGGCAGGGCGATATTCATGGGCGGGTCATTTTAggaa                                                                                              <  2:652736/150‑3 (MQ=255)
                                                       gggAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGAc                                                                                             <  2:128042/150‑1 (MQ=255)
                                                         gAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATAATCCTGGGCCGGTAATTCTAaggcttt                                                                                           >  2:1355546/1‑146 (MQ=255)
                                                          aGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTgg                                                                                          <  1:632632/150‑1 (MQ=255)
                                                                tGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGGGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGcgcc                                                                                    <  1:108008/150‑1 (MQ=255)
                                                                     aGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGCTAAACATGGGCATGGCATTTTAGGGACTGGCGCGCCGCAAg                                                                               >  2:152226/1‑150 (MQ=255)
                                                                     aGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAt                                                                                                      <  2:1836707/127‑1 (MQ=255)
                                                                     aGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAt                                                                                                      >  1:1836707/1‑127 (MQ=255)
                                                                          gcCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGt                                                                                                                                 <  2:1736672/95‑1 (MQ=255)
                                                                          gcCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGt                                                                                                                                 >  1:1736672/1‑95 (MQ=255)
                                                                            cGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGgcgcgc                                                                        <  1:600895/150‑1 (MQ=255)
                                                                              aGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCg                                                                      <  2:373553/150‑1 (MQ=255)
                                                                                     ccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAggg                                                               <  1:472012/150‑1 (MQ=255)
                                                                                     ccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAggg                                                               >  1:395925/1‑150 (MQ=255)
                                                                                     ccAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCAGGCCATCATGGCCGTCACGAATCCGTTCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAggg                                                               >  1:1671040/1‑150 (MQ=255)
                                                                                        ggACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAtt                                                                                                     <  2:1545536/109‑1 (MQ=255)
                                                                                        ggACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCAtt                                                                                                     >  1:1545536/1‑109 (MQ=255)
                                                                                         gACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTgg                                                                                          <  1:902879/119‑1 (MQ=255)
                                                                                         gACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTgg                                                                                          >  2:902879/1‑119 (MQ=255)
                                                                                               ccccAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTc                                                     >  1:1236740/1‑150 (MQ=255)
                                                                                                  cagcagCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGCACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGAt                                                  <  1:519577/150‑1 (MQ=255)
                                                                                                        caAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGaa                                            >  2:820545/1‑150 (MQ=255)
                                                                                                         aaGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTa                                                                                                  <  2:596168/95‑1 (MQ=255)
                                                                                                         aaGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTa                                                                                                  >  1:596168/1‑95 (MQ=255)
                                                                                                          aGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGATCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACg                                          <  1:152226/150‑1 (MQ=255)
                                                                                                              gctttgcgcCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGaa                                      <  1:1648308/145‑1 (MQ=255)
                                                                                                              gcAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGaa                                      <  1:1450518/150‑1 (MQ=255)
                                                                                                                 aaGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTa                                   >  2:1668000/1‑150 (MQ=255)
                                                                                                                  aGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTAc                                  <  1:680760/150‑1 (MQ=255)
                                                                                                                      ccAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGt                              <  2:519617/150‑1 (MQ=255)
                                                                                                                          ttttACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGtt                          <  2:1610936/150‑1 (MQ=255)
                                                                                                                               cGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGa                        >  1:587889/1‑147 (MQ=255)
                                                                                                                               cGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGa                        <  2:587889/147‑1 (MQ=255)
                                                                                                                               cGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGAttt                     <  1:1668000/150‑1 (MQ=255)
                                                                                                                                ggTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGGACTTGTCGTTGAtttt                    >  2:1903634/1‑150 (MQ=255)
                                                                                                                                  tcatcaGGGTGCCGGCCATCATGCCAAGTGCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCt                  <  1:941884/150‑1 (MQ=255)
                                                                                                                                  gcatgaGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCCATGCTGAACGGGAAGTACTTGTCGTTGATTTTCt                  >  2:179343/2‑150 (MQ=255)
                                                                                                                                   catcaGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCtt                 <  2:1287403/150‑1 (MQ=255)
                                                                                                                                       aGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTgtac             >  2:1580714/1‑150 (MQ=255)
                                                                                                                                          gTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTgtacgta          <  2:395925/150‑1 (MQ=255)
                                                                                                                                           gtccGGCCATGATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAgg                                                                                                >  2:785409/3‑63 (MQ=255)
                                                                                                                                           gtccGGCCATGATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAgg                                                                                                <  1:785409/61‑1 (MQ=255)
                                                                                                                                              cGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGt      <  1:91793/150‑1 (MQ=255)
                                                                                                                                               ggCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGTCGCGCCGCAAGGGCGCGCCGGGCAGGGCGGATCATTCGATGCTGCACGCGAAGAACTTGTGGTTGATTTTCTTGTACGTACCGTc     >  1:576230/1‑150 (MQ=255)
                                                                                                                                                gCCATGATGGCCGGCCCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCg    >  1:1139296/1‑150 (MQ=255)
                                                                                                                                                 ccATGATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCgg   >  1:1195824/1‑150 (MQ=255)
                                                                                                                                                  catgatGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGc  >  2:1210660/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACGCCGCGAACCAGGGTCGAGTAGAAGCCGCCAAGCCATTGCAGCGGCTTGAACGGGGAGGTCTTGGCCAGGGCGCCGAGCAGGCCCAGGACCAGCCCCAGCAGCAAGGCGCAAAGCGCCAGTTTTACGGTCATCAGGGTGCCGGCCATCATGGCCGGACCGAATCCGTGCAGGTCGATATTCATGGGCAGGTCATTTTAGGGACTGGCGCGCCGCAAGGGCGCGCCGGTCAGGGCGGATCATTCGATGCTGAACGGGAAGTACTTGTCGTTGATTTTCTTGTACGTACCGTCGGC  >  NC_002947/339079‑339374

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