Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,395,791 A→C 5.3% H9Q (CAT→CAG)  PP_3871 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,395,7910AC 5.3% 250.4 / 4.7 86H9Q (CAT→CAG) PP_3871hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base C (3/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTAGTTGGCCCGAG  >  NC_002947/4395642‑4395938
                                                                                                                                                     |                                                                                                                                                   
aCTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCGGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGa                                                                                                                                                     <  1:33961/150‑1 (MQ=255)
  tGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATg                                                                                                                                                   >  2:780175/1‑150 (MQ=255)
  tGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATg                                                                                                                                                   >  2:666525/1‑150 (MQ=255)
  ngCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATg                                                                                                                                                   >  2:1089044/2‑150 (MQ=255)
     gTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATgcgc                                                                                                                                                <  1:250664/150‑1 (MQ=255)
       tCTGCGCCTCATTCTCCCCGGCTTTGTCATCACCGGTGACCACATCGATCTGCGCCGTCTTTTTTAGCGTGGACGAGCTGACGGCGTCGAGGTCGAGGCCCGCCGAGACCCTGTCGGCCAGGCGCGGGCAGTACTCGCGCACATGCGCGa                                                                                                                                              <  2:942943/150‑1 (MQ=255)
          gcgcCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCTTCGAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCATCTGCGCGagta                                                                                                                                           <  1:1477335/150‑4 (MQ=255)
            gcCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGa                                                                                                                                           >  2:1860211/1‑148 (MQ=255)
            gcCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGa                                                                                                                                           <  1:1860211/148‑1 (MQ=255)
            gcCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACg                                                                                                                                         <  2:994354/150‑1 (MQ=255)
             ccTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACgg                                                                                                                                        <  1:1433061/150‑1 (MQ=255)
              cTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTTGACGAGCGGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACggg                                                                                                                                       <  1:152846/150‑1 (MQ=255)
              cTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACggg                                                                                                                                       <  2:4314/150‑1 (MQ=255)
                     tcGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGAtt                                                                                                                                >  2:1067362/1‑150 (MQ=255)
                     tcGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGAtt                                                                                                                                >  1:596045/1‑150 (MQ=255)
                     tcGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGAtt                                                                                                                                >  1:1292561/1‑150 (MQ=255)
                          ggCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCAt                                                                                                                           <  2:1872079/150‑1 (MQ=255)
                          ggCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCAt                                                                                                                           <  1:666525/150‑1 (MQ=255)
                              ttGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTg                                                                                                                       <  2:1710174/150‑1 (MQ=255)
                               tGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTgg                                                                                                                      <  2:217129/150‑1 (MQ=255)
                                 tcatcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGtt                                                                                                                    >  1:439282/1‑150 (MQ=255)
                                   atcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGGGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTa                                                                                                                  >  2:170561/1‑150 (MQ=255)
                                   atcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATg                                                                                                                          <  2:1454383/142‑1 (MQ=255)
                                   atcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATg                                                                                                                          >  1:1454383/1‑142 (MQ=255)
                                    tcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACgg                                                                                                                                        <  2:1503787/127‑1 (MQ=255)
                                    tcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACgg                                                                                                                                        >  1:1503787/1‑127 (MQ=255)
                                    tcaGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTAc                                                                                                                 >  2:1607733/1‑150 (MQ=255)
                                       gCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTAcc                                                                                                                >  1:1902837/1‑148 (MQ=255)
                                       gCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTAcc                                                                                                                <  2:1902837/148‑1 (MQ=255)
                                            gATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCATCTGCGCGa                                                                                                                                              <  1:1915827/113‑1 (MQ=255)
                                            gATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCATCTGCGCGa                                                                                                                                              >  2:1915827/1‑113 (MQ=255)
                                             aTCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGa                                                                                                        <  2:596045/150‑1 (MQ=255)
                                               cacaTAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAgg                                                                                                      <  2:1761950/150‑1 (MQ=255)
                                                   tAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCGGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTg                                                                                                  >  2:705172/1‑150 (MQ=255)
                                                   tAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTg                                                                                                  <  1:114897/150‑1 (MQ=255)
                                                    aGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGa                                                                                                 <  1:961696/150‑1 (MQ=255)
                                                      gCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATg                                                                                               <  1:780175/150‑1 (MQ=255)
                                                                        aGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACgcggcgg                                                                             <  1:850865/150‑1 (MQ=255)
                                                                           gTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGa                                                                          <  1:1925444/150‑1 (MQ=255)
                                                                             ggACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAg                                                                        <  2:1231288/150‑1 (MQ=255)
                                                                              gACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGCTTTTGATGCTGGTTTACCGGCCCTGAGGCCTGAGGTTCAGCGCAACGCGGCGGCGAAgg                                                                       >  2:119146/1‑150 (MQ=255)
                                                                                       aCGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATg                                                                                               >  2:369716/1‑117 (MQ=255)
                                                                                       aCGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATg                                                                                               <  1:369716/117‑1 (MQ=255)
                                                                                          gCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAgg                                                                       <  2:645682/138‑1 (MQ=255)
                                                                                          gCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAgg                                                                       >  1:645682/1‑138 (MQ=255)
                                                                                               gAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGAt                                                                                                >  1:969459/1‑108 (MQ=255)
                                                                                               gAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGAt                                                                                                <  2:969459/108‑1 (MQ=255)
                                                                                               gAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGAt                                                      <  1:705172/150‑1 (MQ=255)
                                                                                               gAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGAt                                                      <  1:550383/150‑1 (MQ=255)
                                                                                               gAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGAt                                                      <  1:550285/150‑1 (MQ=255)
                                                                                                 ggTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTg                                                    <  2:1292561/150‑1 (MQ=255)
                                                                                                  gTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTgg                                                   >  2:1694510/1‑150 (MQ=255)
                                                                                                      aTGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGGAGTACTCGCGCAGCTGCGCGATGGCGGGGGAGATTTTCATGCTGGTTTCCCGTCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCt                                               >  1:442253/1‑150 (MQ=255)
                                                                                                          cgccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTc                                           >  2:1590759/1‑150 (MQ=255)
                                                                                                          cgccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTc                                           <  1:1607733/150‑1 (MQ=255)
                                                                                                          cgccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTc                                           >  1:894316/1‑150 (MQ=255)
                                                                                                          cgccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAAAGCGATAGCGAACGCTGCTCGGAGTATCGAGTGGACCTGGTc                                           >  2:989962/1‑150 (MQ=255)
                                                                                                           gccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGACGAAGgctgct                                                                  <  2:516206/126‑1 (MQ=255)
                                                                                                           gccCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGACGAAGgctgct                                                                  >  1:516206/1‑126 (MQ=255)
                                                                                                            cccGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGAc                                                 <  1:368882/142‑1 (MQ=255)
                                                                                                            cccGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGAc                                                 >  2:368882/1‑142 (MQ=255)
                                                                                                                  aTCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGt                                   >  1:351511/1‑150 (MQ=255)
                                                                                                                   tCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGAt                                                          >  2:599620/1‑126 (MQ=255)
                                                                                                                   tCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGAt                                                          <  1:599620/126‑1 (MQ=255)
                                                                                                                   tCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGAACTGGTCTTTTTCGTc                                  >  2:431865/1‑150 (MQ=255)
                                                                                                                      tGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGTCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTg                                                    <  2:1130692/129‑1 (MQ=255)
                                                                                                                      tGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTg                                                    >  1:1130692/1‑129 (MQ=255)
                                                                                                                       gTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGAt                                                          >  1:305134/1‑122 (MQ=255)
                                                                                                                       gTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGAt                                                          <  2:305134/122‑1 (MQ=255)
                                                                                                                        tCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGCGGCCTGATGGTCAGCGCAACGCGGCGACGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCa                             >  1:1203631/1‑150 (MQ=255)
                                                                                                                          ggcccggcccGCGCAGAGCGCGCCGAGATGCGCTATGACGGGGGAGATTTTCATGCTGGTTTCCAGTCCCCCACGCCTGATGGTCAGCGCAACGCGGCGGNGNAGGCTGNTCTGAGGATCGANTGGACCTGGTCTTTTTCGTCCTGCNgg                           <  1:1089044/150‑1 (MQ=255)
                                                                                                                           gCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAggg                          >  1:125705/1‑150 (MQ=255)
                                                                                                                                  tCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGc                         >  1:1283484/1‑144 (MQ=255)
                                                                                                                                  tCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGc                         <  2:1283484/144‑1 (MQ=255)
                                                                                                                                   cGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTcc                  <  2:1635978/150‑1 (MQ=255)
                                                                                                                                     ggCAGTACTCGCGCAGATGCGCTAAGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCACCGCGGCGCCGAATGCTGCTTTGAGGAGCGATTGGTCCTGGTCTTTGTCGTCCTGCAGGGCGTCTTCCAt                <  2:606707/150‑1 (MQ=255)
                                                                                                                                     ggCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCCGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCAt                >  1:827277/1‑150 (MQ=255)
                                                                                                                                      gCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGGCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATg               >  1:597610/1‑150 (MQ=255)
                                                                                                                                      gCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACgcggcggc                                                                            >  2:1384585/1‑89 (MQ=255)
                                                                                                                                      gCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACgcggcggc                                                                            <  1:1384585/89‑1 (MQ=255)
                                                                                                                                        aGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTa             >  2:1405658/1‑150 (MQ=255)
                                                                                                                                             tCGCGCATCTGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTAGTTgg        >  2:814641/1‑150 (MQ=255)
                                                                                                                                             tCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGAt                                                                                                <  1:1124301/62‑1 (MQ=255)
                                                                                                                                             tCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGAt                                                                                                >  2:1124301/1‑62 (MQ=255)
                                                                                                                                              cgcgCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTGGTTGGc       <  1:1590759/150‑1 (MQ=255)
                                                                                                                                               gcgcAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTAGTTGGcc      >  2:1250744/1‑150 (MQ=255)
                                                                                                                                                   aGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGACGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTAGTTGGCCCGAg  <  2:1203631/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCAGTAGCGTGGACGAGCTGACGGCGTCGAGGTCGATGCCGCCCGAGATCCTGTCGGCCAGGCTCGGGCAGTACTCGCGCAGATGCGCGATGACGGGGGAGATTTTCATGCTGGTTTACCGGCCCTGAGGCCTGATGGTCAGCGCAACGCGGCGGCGAAGGCTGCTCTGAGGATCGATTGGACCTGGTCTTTTTCGTCCTGCAGGGCGTCTTCCATGTAGTTGGCCCGAG  >  NC_002947/4395642‑4395938

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