Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,034,963 T→C 5.4% L333L (TTG→CTG)  PP_0892 → phospholipase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,034,9630TC 5.4% 215.6 / 2.2 75L333L (TTG→CTG) PP_0892phospholipase family protein
Reads supporting (aligned to +/- strand):  ref base T (40/29);  new base C (2/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

GCGTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAC  >  NC_002947/1034817‑1035112
                                                                                                                                                  |                                                                                                                                                     
gCGTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGa                                                                                                                                                    <  2:27728/150‑1 (MQ=255)
  gTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACa                                                                                                                                                  >  2:923530/1‑150 (MQ=255)
  gTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACa                                                                                                                                                  >  2:835155/1‑150 (MQ=255)
  gTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACGTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACa                                                                                                                                                  >  2:999189/1‑150 (MQ=255)
   tACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAg                                                                                                                                                 <  1:509546/150‑1 (MQ=255)
    aCGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAgc                                                                                                                                                <  1:1825400/150‑1 (MQ=255)
     cGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAgcg                                                                                                                                               >  2:1769019/1‑150 (MQ=255)
             cGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCggtggt                                                                                                                                       <  2:1734433/150‑1 (MQ=255)
              gAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGgtggtg                                                                                                                                      >  2:765257/1‑150 (MQ=255)
                gTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCTAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGc                                                                                                                                    <  1:651109/150‑1 (MQ=255)
                   ccAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAg                                                                                                                                 >  2:785420/1‑150 (MQ=255)
                    cAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGc                                                                                                                                >  1:1547447/1‑150 (MQ=255)
                      gCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGGTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCtt                                                                                                                              >  2:1456284/1‑150 (MQ=255)
                       ccGTGCTTTCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTc                                                                                                                             <  1:26033/150‑1 (MQ=255)
                       ccGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTc                                                                                                                             <  2:105216/150‑1 (MQ=255)
                         gTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGa                                                                                                                           <  1:437803/150‑1 (MQ=255)
                             ttCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCATCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgcgc                                                                                                                       <  2:1166124/150‑1 (MQ=255)
                                 tGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACt                                                                                                                   >  2:839374/1‑150 (MQ=255)
                                    aCCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCg                                                                                                                >  2:1335002/1‑150 (MQ=255)
                                      cTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGcc                                                                                                              >  1:1072242/1‑150 (MQ=255)
                                          aGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                          <  1:590598/150‑1 (MQ=255)
                                            aTGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCAAACCTCGAAGCGCTCGCCCCGCCTTACACACCAACTATCGGATCCATCTTCGAGCGCGACTTCGCCCAGAGc                                                                                                        >  2:1366548/1‑150 (MQ=255)
                                              ggCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCga                                                                                                      >  1:199841/1‑150 (MQ=255)
                                               gCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                     <  1:1135058/150‑1 (MQ=255)
                                               gCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                     >  2:1857667/1‑150 (MQ=255)
                                               gCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                     <  2:139734/150‑1 (MQ=255)
                                                     gTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTc                                                                                               <  1:199173/150‑1 (MQ=255)
                                                     gTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTc                                                                                               <  1:1366548/150‑1 (MQ=255)
                                                      tAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAATCGCTCGACCCGCCGTTGACAGCTGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCg                                                                                              <  1:705410/150‑1 (MQ=255)
                                                        gacgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGCCAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCACAGCGAAGGAGTcaaa                                                                                            <  2:325181/150‑4 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGtt                                                                                                                                                      >  1:1838576/1‑89 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGtt                                                                                                                                                      <  2:1838576/89‑1 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                         >  1:1027062/1‑150 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                         >  1:1264855/1‑150 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                         >  2:1901745/1‑150 (MQ=255)
                                                           gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                         >  1:1542867/1‑150 (MQ=255)
                                                                ggTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACATGAAAGTGGGCCTGCACCCCACCCTCGAAGCGCCCCACCCGCCCTCGACAGCCGCGGTGCTGGCCACCTTCGAGCGCGACTTCGCCCCGAGCGAACACGTCCACCCCGcccc                                                                                    >  2:65722/1‑149 (MQ=255)
                                                                   cAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTg                                                                                 <  2:1224790/150‑1 (MQ=255)
                                                                           gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATgcgc                                                                         >  2:1515282/1‑150 (MQ=255)
                                                                            ttCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATgcgcg                                                                        >  2:305166/1‑150 (MQ=255)
                                                                                tGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGGGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCt                                                                    >  2:859625/1‑150 (MQ=255)
                                                                                 gCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCtt                                                                   <  1:785420/150‑1 (MQ=255)
                                                                                  cAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtg                                                                  >  1:390978/1‑150 (MQ=255)
                                                                                   aaCTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtgg                                                                 >  1:1724227/1‑150 (MQ=255)
                                                                                   aaCTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtgg                                                                 <  1:1857667/150‑1 (MQ=255)
                                                                                     cTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtggtg                                                               <  2:861853/150‑1 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTGGCTACAGCGGTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTAGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                          >  1:1789701/1‑150 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTg                                                                                 <  2:754774/127‑1 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTg                                                                                 >  1:754774/1‑127 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                          <  1:1769019/150‑1 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                          <  1:1335002/150‑1 (MQ=255)
                                                                                          aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                          >  2:852711/1‑150 (MQ=255)
                                                                                           ctaCTGGGACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAg                                                         >  2:1714599/3‑150 (MQ=255)
                                                                                             aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCAACTGGCATGCGCGCCCTTGGGGGCGCAggg                                                       >  2:851298/1‑150 (MQ=255)
                                                                                               tGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTg                                                     >  2:356858/1‑150 (MQ=255)
                                                                                                 gAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGaa                                                   >  1:263916/1‑150 (MQ=255)
                                                                                                         gcTTCAACCTTGAAGCCGTCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGa                                           >  1:1096858/1‑150 (MQ=255)
                                                                                                                ccTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGgtggtg                                                                                                                                      >  1:1209509/1‑52 (MQ=255)
                                                                                                                ccTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGgtggtg                                                                                                                                      <  2:1209509/52‑1 (MQ=255)
                                                                                                                ccTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgggg                                    >  2:1674285/1‑150 (MQ=255)
                                                                                                                 cTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggc                                   >  2:1641069/1‑150 (MQ=255)
                                                                                                                 cTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggc                                   >  1:1631487/1‑150 (MQ=255)
                                                                                                                  ttGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggct                                  >  1:274968/1‑150 (MQ=255)
                                                                                                                       gCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTc                             <  2:263916/150‑1 (MQ=255)
                                                                                                                          aaCATCGAAGCGCTCGAACCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                          >  2:1155983/1‑150 (MQ=255)
                                                                                                                           aCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATc                         <  1:999189/150‑1 (MQ=255)
                                                                                                                                gacgcgcTCGACGCGGCGGTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                    <  2:1542867/147‑1 (MQ=255)
                                                                                                                                gAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                    <  1:765257/150‑1 (MQ=255)
                                                                                                                                 aaGCGCTCGACCCGCCGTTGACAGCGGGCCAGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtgg                   >  1:51850/1‑150 (MQ=255)
                                                                                                                                 aaGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtgg                   <  1:305166/150‑1 (MQ=255)
                                                                                                                                      cTNGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGNGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc              >  1:1286883/1‑150 (MQ=255)
                                                                                                                                         ccaccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc              >  1:1664380/4‑147 (MQ=255)
                                                                                                                                         ccaccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc              <  2:1664380/144‑1 (MQ=255)
                                                                                                                                          caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                    <  2:444338/138‑1 (MQ=255)
                                                                                                                                          caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                    >  1:444338/3‑140 (MQ=255)
                                                                                                                                            ccgccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTc        >  1:1429314/1‑150 (MQ=255)
                                                                                                                                               ccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc     <  1:356858/150‑1 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc     <  2:1626744/147‑1 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc     <  2:1626698/147‑1 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc     >  1:1626744/1‑147 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc     >  1:1626698/1‑147 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAc  <  1:1641069/150‑1 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGGTGGTGGCCAGCTTCCAGCAAGAGGTGGCACAAAGCGAAGAATTCGAAATCGCCCACTGGAATAAGAAGGTATACATATTTATGGTTAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGATAACGTCCTCGAc  >  1:1876699/1‑150 (MQ=255)
                                                                                                                                                  ttGACAGCGGCGATGGTGGCCCGCTTCTAGCGCGACTTCCCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATATGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCGCGAc  <  1:388698/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAC  >  NC_002947/1034817‑1035112

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: