Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,610,345 C→G 8.6% A262A (GCC→GCG PP_4081 → type IV/VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,610,3450CG 8.6% 176.2 / 8.4 58A262A (GCC→GCGPP_4081type IV/VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base G (3/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCT  >  NC_002947/4610197‑4610487
                                                                                                                                                     |                                                                                                                                              
tGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCcgccatcgcca                                                                                                                                               >  1:882972/1‑150 (MQ=255)
      aCCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCATAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCgg                                                                                                                                         >  2:744960/1‑150 (MQ=255)
      aCCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCgg                                                                                                                                         >  1:1056298/1‑150 (MQ=255)
       cgaaCTTCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGt                                                                                                                                        >  2:857029/3‑150 (MQ=255)
       ccAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGt                                                                                                                                        <  2:899357/150‑1 (MQ=255)
       ccAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGt                                                                                                                                        >  2:223063/1‑150 (MQ=255)
        cAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTg                                                                                                                                       <  1:30418/150‑1 (MQ=255)
        cAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAACAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTg                                                                                                                                       >  1:615501/1‑150 (MQ=255)
         aaCTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCCACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTga                                                                                                                                      <  1:476538/150‑2 (MQ=255)
               aCGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACg                                                                                                                                <  1:215281/150‑1 (MQ=255)
                                 tGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCGCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGCGGAGCCCGTTGGTCACCACCGCCATCGCCATGGCGGTGGCCTACGCGATTTCCAGCCAACGcc                                                                                                              >  2:270998/1‑150 (MQ=18)
                                         cccACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTg                                                                                                      >  2:949801/1‑150 (MQ=255)
                                          ccACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGa                                                                                                     <  2:1079706/150‑1 (MQ=255)
                                          ccACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGa                                                                                                     <  1:949801/150‑1 (MQ=255)
                                                ccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATcac                                                                                                 <  2:882604/148‑1 (MQ=255)
                                                ccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATcac                                                                                                 >  1:882604/1‑148 (MQ=255)
                                                ccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACaa                                                                                               >  2:742875/1‑150 (MQ=255)
                                                ccGAAGAAGTCTCAGACCCTAACACCAACTTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACaa                                                                                               <  2:615501/150‑1 (MQ=18)
                                                     gaagTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGCTTTACAGCTACCGCCTGCACCTGATCACAACCGAg                                                                                          <  1:869169/150‑1 (MQ=255)
                                                     gaagTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAg                                                                                          <  1:930231/150‑1 (MQ=255)
                                                           tGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTg                                                                                    <  1:178957/150‑1 (MQ=255)
                                                             cGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGa                                                                                  <  2:1056298/150‑1 (MQ=255)
                                                                  ccTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGc                                                                             <  1:842479/150‑1 (MQ=255)
                                                                   cTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCt                                                                            >  2:233617/1‑150 (MQ=255)
                                                                   cTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCAGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCt                                                                            >  2:175464/1‑150 (MQ=255)
                                                                      cacCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGaa                                                                         <  2:893707/150‑1 (MQ=255)
                                                                       acCAACGTCGCCCCCAGAAACTTGAGCTTGGGCCTGATATGGCCGTGGTGGAGCCCGTTCGTCATCGATCGGTGTGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAAc                                                                        <  1:223063/150‑1 (MQ=18)
                                                                          aaCGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTg                                                                                                            >  1:777023/1‑111 (MQ=255)
                                                                          aaCGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTg                                                                                                            <  2:777023/111‑1 (MQ=255)
                                                                                gCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTg                                                               <  1:382852/150‑1 (MQ=255)
                                                                                  ccccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                             >  1:220601/1‑150 (MQ=255)
                                                                                    ccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGa                                                                                  <  2:1049821/127‑1 (MQ=255)
                                                                                    ccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGa                                                                                  >  1:1049821/1‑127 (MQ=255)
                                                                                    ccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgcagc                                                           <  2:1216124/150‑1 (MQ=255)
                                                                                     cAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgcagca                                                          <  1:253642/150‑1 (MQ=255)
                                                                                     cAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgcagca                                                          >  2:1088061/1‑150 (MQ=255)
                                                                                                 cATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCgt                                              <  1:857029/150‑1 (MQ=255)
                                                                                                       ccGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTcact                                        <  2:882972/150‑2 (MQ=255)
                                                                                                          gATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                             <  2:966557/126‑1 (MQ=255)
                                                                                                          gATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                             >  1:966557/1‑126 (MQ=255)
                                                                                                                       tggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                        >  2:108925/1‑150 (MQ=255)
                                                                                                                       tggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                        <  1:175464/150‑1 (MQ=255)
                                                                                                                       tggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                        <  2:801005/150‑1 (MQ=40)
                                                                                                                         gAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGcaccac                                 <  2:1102109/139‑1 (MQ=255)
                                                                                                                         gAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGcaccac                                 >  1:1102109/1‑139 (MQ=255)
                                                                                                                              cGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTc                 <  1:824589/150‑1 (MQ=40)
                                                                                                                              cGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTc                 >  1:194236/1‑150 (MQ=39)
                                                                                                                               gTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcg                >  1:239973/1‑150 (MQ=40)
                                                                                                                                ttGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcgc               >  1:470105/1‑150 (MQ=38)
                                                                                                                                  ggTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCCCTAAGATCAACGAATTCGCGa             >  2:614162/1‑150 (MQ=255)
                                                                                                                                    tCATCTCCGCCATCGCCATGGCGGTGGCNTACGGGATTTACAGCNANCNCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAgg           >  1:880331/1‑150 (MQ=37)
                                                                                                                                      aTCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGtt         >  2:143137/1‑150 (MQ=35)
                                                                                                                                       gccaccGCCATGGCGATGGCGGAGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                        >  1:10379/5‑135 (MQ=11)
                                                                                                                                       gccaccGCCATGGCGATGGCGGAGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                        <  2:10379/131‑1 (MQ=11)
                                                                                                                                        ccaccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCgtgt                                            <  1:242200/111‑1 (MQ=255)
                                                                                                                                        ccaccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCgtgt                                            >  2:242200/4‑114 (MQ=255)
                                                                                                                                          tcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGttttt      <  1:744960/150‑1 (MQ=32)
                                                                                                                                          cccGCCATCGCCATGGCGGTGGCCTACCGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTTAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCAGCACTGAGTTCAACGAATTCGCGAGGttttt      <  1:270998/149‑1 (MQ=12)
                                                                                                                                           ccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAATTCAACGAATTCGCGAGGTTTTTg     >  2:1204422/1‑150 (MQ=255)
                                                                                                                                            cgccatcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGc    <  2:470105/150‑1 (MQ=32)
                                                                                                                                            cgccatcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGc    >  2:250836/1‑150 (MQ=32)
                                                                                                                                            cgccatcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGc    >  2:1003325/1‑150 (MQ=32)
                                                                                                                                              ccatcgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCt  >  1:240814/1‑150 (MQ=32)
                                                                                                                                                   gccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcg                >  2:730065/1‑131 (MQ=33)
                                                                                                                                                   gccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcg                <  1:730065/131‑1 (MQ=33)
                                                                                                                                                    ccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcgc               >  2:724903/1‑131 (MQ=33)
                                                                                                                                                    ccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTcgc               <  1:724903/131‑1 (MQ=32)
                                                                                                                                                     |                                                                                                                                              
TGATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCT  >  NC_002947/4610197‑4610487

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: