Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,034,969 G→A 5.4% A335T (GCG→ACG)  PP_0892 → phospholipase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,034,9690GA 5.4% 220.9 / 2.2 76A335T (GCG→ACG) PP_0892phospholipase family protein
Reads supporting (aligned to +/- strand):  ref base G (40/31);  new base A (2/2);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCG  >  NC_002947/1034820‑1035117
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                          <  2:139734/150‑1 (MQ=255)
                                            gCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                          >  2:1857667/1‑150 (MQ=255)
                                            gCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                          <  1:1135058/150‑1 (MQ=255)
                                                  gTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTc                                                                                                    <  1:1366548/150‑1 (MQ=255)
                                                  gTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTc                                                                                                    <  1:199173/150‑1 (MQ=255)
                                                   tAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAATCGCTCGACCCGCCGTTGACAGCTGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCg                                                                                                   <  1:705410/150‑1 (MQ=255)
                                                     gacgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGCCAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCACAGCGAAGGAGTcaaa                                                                                                 <  2:325181/150‑4 (MQ=255)
                                                        gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                              >  1:1264855/1‑150 (MQ=255)
                                                        gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                              >  2:1901745/1‑150 (MQ=255)
                                                        gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                              >  1:1542867/1‑150 (MQ=255)
                                                        gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTc                                                                                              >  1:1027062/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:1857667/150‑1 (MQ=255)
                                                                                aaCTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtgg                                                                      >  1:1724227/1‑150 (MQ=255)
                                                                                  cTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTtggtg                                                                    <  2:861853/150‑1 (MQ=255)
                                                                                       aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTGGCTACAGCGGTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTAGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                               >  1:1789701/1‑150 (MQ=255)
                                                                                       aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTg                                                                                      >  1:754774/1‑127 (MQ=255)
                                                                                       aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTg                                                                                      <  2:754774/127‑1 (MQ=255)
                                                                                       aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                               <  1:1335002/150‑1 (MQ=255)
                                                                                       aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCa                                                               <  1:1769019/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                                        >  1:1631487/1‑150 (MQ=255)
                                                                                                              cTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggc                                        >  2:1641069/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                   <  2:1664380/144‑1 (MQ=255)
                                                                                                                                      ccaccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc                   >  1:1664380/4‑147 (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:1626698/147‑1 (MQ=255)
                                                                                                                                               ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc          >  1:1626744/1‑147 (MQ=255)
                                                                                                                                               ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc          >  1:1626698/1‑147 (MQ=255)
                                                                                                                                               ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc          <  2:1626744/147‑1 (MQ=255)
                                                                                                                                               ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAc       <  1:1641069/150‑1 (MQ=255)
                                                                                                                                               ttGACAGCGGCGGTGGTGGCCAGCTTCCAGCAAGAGGTGGCACAAAGCGAAGAATTCGAAATCGCCCACTGGAATAAGAAGGTATACATATTTATGGTTAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGATAACGTCCTCGAc       >  1:1876699/1‑150 (MQ=255)
                                                                                                                                               ttGACAGCGGCGATGGTGGCCCGCTTCTAGCGCGACTTCCCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATATGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCGCGAc       <  1:388698/150‑1 (MQ=255)
                                                                                                                                                    aGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAccgccg  >  1:1665382/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGAAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACCGCCG  >  NC_002947/1034820‑1035117

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: