Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,710,769 G→A 8.0% V273V (GTC→GTT gpsA ← glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,710,7690GA 8.0% 264.6 / 8.8 96V273V (GTC→GTTgpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base A (6/2);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGAC  >  NC_002947/4710625‑4710909
                                                                                                                                                |                                                                                                                                            
cTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGc                                                                                                                                         >  2:670910/1‑150 (MQ=255)
cTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGc                                                                                                                                         <  2:488502/150‑1 (MQ=255)
cTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGc                                                                                                                                         >  2:117979/1‑150 (MQ=255)
 tCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGcc                                                                                                                                        >  1:897170/1‑150 (MQ=255)
 ataaTCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGcc                                                                                                                                        >  1:761001/4‑150 (MQ=255)
   gATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGCctg                                                                                                                                      >  2:96851/1‑149 (MQ=255)
    aTCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGCctgg                                                                                                                                     <  2:897170/150‑3 (MQ=255)
     tCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCTGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTCAACCGCctggg                                                                                                                                    <  1:96851/150‑4 (MQ=255)
         cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTTCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAg                                                                                                                                >  1:100091/1‑150 (MQ=255)
         cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGAcc                                                                                                                                           <  2:936825/139‑1 (MQ=255)
         cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGAcc                                                                                                                                           >  1:936825/1‑139 (MQ=255)
         cTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAg                                                                                                                                >  2:1471/1‑150 (MQ=255)
            gTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCt                                                                                                                             >  2:675492/1‑150 (MQ=255)
              tCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCttt                                                                                                                           <  2:965698/150‑1 (MQ=255)
                aGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTgg                                                                                                                         >  2:1288907/1‑150 (MQ=255)
                      cgcCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGc                                                                                                                   <  1:1471/150‑1 (MQ=255)
                       gcCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTgg                                                                                                                         <  2:1293048/143‑1 (MQ=255)
                       gcCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTgg                                                                                                                         >  1:1293048/1‑143 (MQ=255)
                       gcCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGcc                                                                                                                  >  1:438335/1‑150 (MQ=255)
                       gcCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGcc                                                                                                                  >  1:817828/1‑150 (MQ=255)
                          cTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCa                                                                                                               >  2:896961/1‑150 (MQ=255)
                           ttCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTc                                                                                                                                   <  1:300223/129‑1 (MQ=255)
                           ttCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTc                                                                                                                                   >  2:300223/1‑129 (MQ=255)
                            tCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGt                                                                                                             <  1:675492/150‑1 (MQ=255)
                               aaCAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCg                                                                                                          >  2:641842/1‑150 (MQ=255)
                                 cAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCGCTTCGCCCAGGCGGTTGACCGCCTCTTTCATGCTTTGGCCCTGCCCCGGTGCGTa                                                                                                        <  2:216798/150‑1 (MQ=255)
                                    gATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGcc                                                                                                     >  1:615526/1‑150 (MQ=255)
                                      tGGCATGCAGGCCAGCCACCAGCGGCATGTACACCCGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCaa                                                                                                   <  2:761001/150‑1 (MQ=255)
                                           tGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGa                                                                                              <  1:767956/150‑1 (MQ=255)
                                           tGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGa                                                                                              >  1:1233121/1‑150 (MQ=255)
                                                   agccaCCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGGGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGATTGGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCg                                                                                      >  2:1046423/1‑150 (MQ=255)
                                                       accaGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAAcc                                                                                                 >  2:475877/1‑135 (MQ=255)
                                                       accaGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAAcc                                                                                                 <  1:475877/135‑1 (MQ=255)
                                                       accaGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGctc                                                                                  >  1:730519/1‑150 (MQ=255)
                                                          aGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTgg                                                                                                                         >  1:71371/1‑108 (MQ=255)
                                                          aGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTgg                                                                                                                         <  2:71371/108‑1 (MQ=255)
                                                          aGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTg                                                                               <  1:641842/150‑1 (MQ=255)
                                                          aGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTg                                                                               <  2:438335/150‑1 (MQ=255)
                                                            cGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTggg                                                                             >  1:1373318/1‑150 (MQ=255)
                                                               cATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCga                                                                          >  2:1350624/1‑150 (MQ=255)
                                                               cATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCga                                                                          <  1:670910/150‑1 (MQ=255)
                                                               cATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCga                                                                          <  2:680570/150‑1 (MQ=255)
                                                                      acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGc                                                                                                                   >  2:1277927/1‑102 (MQ=255)
                                                                      acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGc                                                                                                                   <  1:1277927/102‑1 (MQ=255)
                                                                      acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                   <  1:1350624/150‑1 (MQ=255)
                                                                      acCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                   >  1:998581/1‑150 (MQ=255)
                                                                       ccTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAgg                                                                  <  1:829898/150‑1 (MQ=255)
                                                                       ccTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAgg                                                                  >  1:588511/1‑150 (MQ=255)
                                                                        cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGcc                                                                                                                  >  1:8154/1‑101 (MQ=255)
                                                                        cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGcc                                                                                                                  <  2:8154/101‑1 (MQ=255)
                                                                        cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                   <  1:670786/148‑1 (MQ=255)
                                                                        cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAg                                                                   >  2:670786/1‑148 (MQ=255)
                                                                        cTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGt                                                                 >  1:571867/1‑150 (MQ=255)
                                                                         tGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGc                                                                                                           <  1:762871/107‑1 (MQ=255)
                                                                         tGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGc                                                                                                           >  2:762871/1‑107 (MQ=255)
                                                                         tGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTg                                                                >  2:1163527/1‑150 (MQ=255)
                                                                          gCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCATCTATAGAGTTGCGGCTCTTGGGCGAGGAGCAGGTGa                                                               <  1:1046423/150‑1 (MQ=255)
                                                                              ttGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGAt                                                           >  2:545663/1‑150 (MQ=255)
                                                                                gCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGCCCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCa                                                         <  1:838803/150‑1 (MQ=255)
                                                                                 cTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAg                                                        <  2:914689/150‑1 (MQ=255)
                                                                                 cTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAg                                                        <  1:545663/150‑1 (MQ=255)
                                                                                       cTTTGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCAcc                                                  <  2:588511/150‑1 (MQ=255)
                                                                                                 aGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCt                                                                                                                     >  2:192742/1‑73 (MQ=255)
                                                                                                 aGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCt                                                                                                                     <  1:192742/73‑1 (MQ=255)
                                                                                                     ccTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGc                                    <  2:100091/150‑1 (MQ=255)
                                                                                                      cTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGAt                                                           >  1:210580/1‑126 (MQ=255)
                                                                                                      cTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGAt                                                           <  2:210580/126‑1 (MQ=255)
                                                                                                      cTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGATGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGcc                                   >  1:1356517/1‑150 (MQ=255)
                                                                                                        tGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACGGACACCGGCCAGGCCGa                                 <  2:730519/150‑1 (MQ=255)
                                                                                                           gCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCAcc                                                  <  2:482717/130‑1 (MQ=255)
                                                                                                           gCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCAcc                                                  >  1:482717/1‑130 (MQ=255)
                                                                                                              tgtTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACg                           >  1:989236/1‑150 (MQ=255)
                                                                                                               gtTGACCCCCTCGGCCACTTCGCCCAGGCGGTTAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGTGCGAGGAGCAGGTGACGCTCAGGTTACCTACACCGGCCAGGCCGAGCAACgc                          >  2:242685/1‑149 (MQ=255)
                                                                                                                ttGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTc                         >  1:680745/1‑150 (MQ=255)
                                                                                                                   aCCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTTTTCCAGGCTTTGACCCTGCCCCAATGCGTCGCCAACCTGATAGTTGCGTCTCTTTTGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATg                      >  1:194139/1‑150 (MQ=255)
                                                                                                                   aCCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAg                                <  1:177723/140‑1 (MQ=255)
                                                                                                                   aCCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAg                                >  2:177723/1‑140 (MQ=255)
                                                                                                                   aCCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATg                      <  2:680745/150‑1 (MQ=255)
                                                                                                                    cccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgag                                                                         >  2:1390908/1‑98 (MQ=255)
                                                                                                                    cccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATgg                     <  1:1288907/150‑1 (MQ=255)
                                                                                                                    cccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCCTCCAGGCTTTGGCACTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgag                                                                         <  1:1390908/98‑1 (MQ=255)
                                                                                                                     ccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCgg                                                                                     >  2:331305/1‑85 (MQ=255)
                                                                                                                     ccccTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCgg                                                                                     <  1:331305/85‑1 (MQ=255)
                                                                                                                         tCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTg                <  1:816138/150‑1 (MQ=255)
                                                                                                                           ggCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGgc              <  2:998581/150‑1 (MQ=255)
                                                                                                                              cACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgag                                                                         >  2:1313837/1‑88 (MQ=255)
                                                                                                                              cACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCgag                                                                         <  1:1313837/88‑1 (MQ=255)
                                                                                                                                 ttCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTg        <  2:571867/150‑1 (MQ=255)
                                                                                                                                   cGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGc      <  2:817828/150‑1 (MQ=255)
                                                                                                                                      agaGGCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGGGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGa   >  1:578356/3‑150 (MQ=255)
                                                                                                                                      agaGGCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGa   >  1:1111530/3‑150 (MQ=255)
                                                                                                                                       gaGGCGGTCAACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGAc  >  2:476653/2‑150 (MQ=255)
                                                                                                                                       cAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGAc  >  1:643132/1‑150 (MQ=255)
                                                                                                                                       cAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGAc  <  2:1233121/150‑1 (MQ=255)
                                                                                                                                         ggCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGc                                    >  2:759156/1‑114 (MQ=255)
                                                                                                                                         ggCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGc                                    <  1:759156/114‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CTCGATCACCTGGTTCAGCGTGCGCCCTTCGAACAGGATGGCATGCAGGCCAGCCACCAGCGGCATGTACACCTGCACTTGCTGGGCCTTGGTCTTGAGCACCTTGAGCGTGTTGACCCCCTCGGCCACTTCGCCCAGGCGGTTGACCGCCTCTTCCAGGCTTTGGCCCTGCCCCAGTGCGTAGCCAACCTGATAGTTGCGGCTCTTGGGCGAGGAGCAGGTGACGATCAGGTCACCGACACCGGCCAGGCCGAGGAACGTCATGGGGTTGGCGCCCTGGCTGAC  >  NC_002947/4710625‑4710909

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