Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,554,974 C→T 6.8% intergenic (‑31/+140) PP_4886 ← / ← PP_4887 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,554,9740CT 6.8% 301.6 / 8.2 102intergenic (‑31/+140)PP_4886/PP_4887hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (48/47);  new base T (2/5);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 108 total aligned reads displayed.

CTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGATGCTGCTTAAAGCTTC  >  NC_002947/5554830‑5555123
                                                                                                                                                |                                                                                                                                                     
cTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGTTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTa                                                                                                                                                  <  1:1565157/150‑1 (MQ=255)
 ttATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTac                                                                                                                                                 >  1:1163053/1‑150 (MQ=255)
 ttATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTAACTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTac                                                                                                                                                 >  1:865182/1‑150 (MQ=255)
    tCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCTGCTGATAACGACGATGTTCAGAAATACACATCTTAACTTTGCTCCCGGAGTAAAACGGCGGCCGGTaccga                                                                                                                                              >  1:1737946/1‑150 (MQ=255)
    tCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTaccga                                                                                                                                              <  2:1702955/150‑1 (MQ=255)
       aGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTaccgaccg                                                                                                                                           >  2:1881441/1‑150 (MQ=255)
         ttcTTCTCCAGTGATTTCCGCGATATGAACATAAGCATCGACAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAACCGGCGGCCGGTACCGACCGcc                                                                                                                                         <  2:1810539/148‑1 (MQ=255)
           cTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgt                                                                                                                                       <  2:801008/150‑3 (MQ=255)
           cTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCCCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgt                                                                                                                                       <  1:1495248/150‑3 (MQ=255)
            ttCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGCTAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                                                                      <  1:1974419/150‑4 (MQ=255)
            ttCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATTCAAGCGCCTGTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTc                                                                                                                                      >  1:1022058/1‑150 (MQ=255)
            ttCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTc                                                                                                                                      >  1:229545/1‑150 (MQ=255)
             tCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCa                                                                                                                                     <  2:108318/150‑1 (MQ=255)
               gCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAgg                                                                                                                                   <  2:1591872/150‑1 (MQ=255)
                ccAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGcc                                                                                                                                         <  2:311158/143‑1 (MQ=255)
                ccAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGcc                                                                                                                                         >  1:311158/1‑143 (MQ=255)
                ccAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGc                                                                                                                                  >  2:683394/1‑150 (MQ=255)
                ccAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGc                                                                                                                                  >  1:1709674/1‑150 (MQ=255)
                  aGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                                                                      >  1:493018/1‑141 (MQ=255)
                  aGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                                                                      <  2:493018/144‑4 (MQ=255)
                     gATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAAt                                                                                                                             >  2:1011191/1‑150 (MQ=255)
                         tCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATacaa                                                                                                                         >  1:972145/1‑150 (MQ=255)
                         tCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATacaa                                                                                                                         <  1:1011191/150‑1 (MQ=255)
                               atatGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTaa                                                                                                                   <  2:582248/150‑1 (MQ=255)
                                 atGAACATAAGCATCGAAAGCCACGCTTACCTCTCCCTGGTAAAAGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGACGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAAcc                                                                                                                 >  2:1178934/1‑150 (MQ=255)
                                      cATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTAt                                                                                                            <  2:619691/150‑1 (MQ=255)
                                       aTAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATg                                                                                                           >  2:2240368/1‑150 (MQ=255)
                                          aGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCa                                                                                                        <  2:229545/150‑1 (MQ=255)
                                            cgtCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                                                                      <  2:758279/148‑1 (MQ=255)
                                              tCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGccc                                                                                                    <  1:1848965/150‑1 (MQ=255)
                                                 aaaGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTg                                                                                                 <  2:1163053/150‑1 (MQ=255)
                                                    gCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATacaa                                                                                                                         <  2:1741850/123‑1 (MQ=255)
                                                    gCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATacaa                                                                                                                         >  1:1741850/1‑123 (MQ=255)
                                                              ccTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGc                                                                                    <  1:2036449/150‑1 (MQ=255)
                                                                tctcCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGccc                                                                                  >  1:1060912/1‑150 (MQ=255)
                                                                tctcCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGccc                                                                                  <  1:2240368/150‑1 (MQ=255)
                                                                  tcCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTa                                                                                >  1:1074950/1‑150 (MQ=255)
                                                                  tcCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTa                                                                                <  2:180857/150‑1 (MQ=255)
                                                                      tGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCaa                                                                            <  2:1992475/150‑1 (MQ=255)
                                                                      tGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCaa                                                                            <  1:1881441/150‑1 (MQ=255)
                                                                          aaaGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGaa                                                                        >  1:2060191/1‑150 (MQ=255)
                                                                              gTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              >  1:1830854/1‑108 (MQ=255)
                                                                              gTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              <  2:1830854/108‑1 (MQ=255)
                                                                              gTTCGGCTGATAACAAAGCTGTTCATACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATAtt                                                                    >  2:1419935/1‑150 (MQ=255)
                                                                               ttCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              >  2:1080360/1‑107 (MQ=255)
                                                                               ttCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              <  1:1080360/107‑1 (MQ=255)
                                                                                 cGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 <  2:1022058/150‑1 (MQ=255)
                                                                                 cGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAAAGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 >  1:52606/1‑150 (MQ=255)
                                                                                   gCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAatat                                                                     <  1:1011084/144‑1 (MQ=255)
                                                                                   gCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAatat                                                                     >  2:1011084/1‑144 (MQ=255)
                                                                                   gCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGa                                                               >  1:1027320/1‑150 (MQ=255)
                                                                                      gATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAAt                                                            <  2:865182/150‑1 (MQ=255)
                                                                                         gacgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCg                                                                                                                                                     >  1:2254412/2‑58 (MQ=255)
                                                                                         gacgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCg                                                                                                                                                     <  2:2254412/57‑1 (MQ=255)
                                                                                         aaCGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCCAAGTCCCTATCAAGGAATATTTGGGAAAtctc                                                         <  1:1270193/150‑1 (MQ=255)
                                                                                          acgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTCAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTcc                                                        <  2:1737946/150‑1 (MQ=255)
                                                                                          acgacgATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTcc                                                        <  1:683394/150‑1 (MQ=255)
                                                                                           cgacgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAATGCCCTATCAAGGAATATTTGGGAAAtctc                                                         <  2:1206694/148‑1 (MQ=255)
                                                                                            gacggagTTCAAACGTCGTCATCCTGACGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTa                                                      >  2:93483/7‑150 (MQ=255)
                                                                                            gacgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACa                                                                                            <  2:145641/112‑1 (MQ=255)
                                                                                            gacgaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACa                                                                                            >  1:145641/1‑112 (MQ=255)
                                                                                               gaTGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAg                                                   >  1:546242/1‑150 (MQ=255)
                                                                                                  tgtCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCa                                                                                                                               <  2:694914/69‑1 (MQ=255)
                                                                                                  gTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCa                                                                                                                               >  1:694914/1‑71 (MQ=255)
                                                                                                   ttCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGACTc                                               <  1:2405228/150‑1 (MQ=255)
                                                                                                    tCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTcg                                              >  1:1626658/1‑150 (MQ=255)
                                                                                                     cAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                                                                      >  1:2162154/1‑58 (MQ=37)
                                                                                                     cAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGTCGGTACCGgccgccgtt                                                                                                                                      <  2:2162154/61‑4 (MQ=37)
                                                                                                        aCATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACGGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTcgacga                                          >  1:1761298/1‑150 (MQ=255)
                                                                                                         cATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACgat                                         <  2:2123754/150‑1 (MQ=255)
                                                                                                         cATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACgat                                         <  1:1116255/150‑1 (MQ=255)
                                                                                                              tCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGa                                                                         <  1:1053646/113‑1 (MQ=255)
                                                                                                              tCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGa                                                                         >  2:1053646/1‑113 (MQ=255)
                                                                                                                   cTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAg                                                                           >  1:2339978/1‑106 (MQ=255)
                                                                                                                   cTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAg                                                                           <  2:2339978/106‑1 (MQ=255)
                                                                                                                   cTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 <  1:1101061/116‑1 (MQ=255)
                                                                                                                   cTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 >  2:1101061/1‑116 (MQ=255)
                                                                                                                    tGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGAACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGAttt                              >  1:1038140/1‑150 (MQ=255)
                                                                                                                     gAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTg                             <  2:1074950/150‑1 (MQ=255)
                                                                                                                              cccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTcgacga                                          <  2:2263733/128‑1 (MQ=255)
                                                                                                                              cccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTcgacga                                          >  1:2263733/1‑128 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              <  1:2342122/59‑1 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCt                                                                                                              >  2:2342122/1‑59 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTcc                                                        >  2:2311463/1‑113 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCg                   >  1:712930/1‑150 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCg                   >  2:2338618/1‑150 (MQ=255)
                                                                                                                               ccGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCg                   >  1:1723272/1‑150 (MQ=255)
                                                                                                                                cGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGa                  >  1:1064931/1‑150 (MQ=255)
                                                                                                                                 gAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                                                                      <  2:2009521/65‑1 (MQ=255)
                                                                                                                                 gAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGc                                                                                                      >  1:2009521/1‑65 (MQ=255)
                                                                                                                                 gAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGCCGATTTGGTCTGTTTCGAt                 <  2:712930/150‑1 (MQ=255)
                                                                                                                                 gAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGAt                 <  2:546242/150‑1 (MQ=255)
                                                                                                                                    gcctgaaggcggtCGGTACCGGCCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGAtgct              >  2:729548/8‑150 (MQ=255)
                                                                                                                                    gTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGAtgct              >  2:2093193/1‑150 (MQ=255)
                                                                                                                                    gTAAAACGGCGGCCGGTACCGACCGCCTTCACGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGTAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGAtgct              <  2:1060912/150‑1 (MQ=255)
                                                                                                                                     cctgaaggcggtCGGTACCGGCCGACTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGAtgctg             >  1:1771810/7‑150 (MQ=255)
                                                                                                                                       aaaCGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTAACGTTGCGACGATGATGACGATTTGGTCAGTTTCGATgctgct           >  1:1833515/1‑150 (MQ=255)
                                                                                                                                             cggcCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGATGCTGCTTAAAGc     <  1:1419935/150‑1 (MQ=255)
                                                                                                                                               gcCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 <  2:743215/88‑1 (MQ=255)
                                                                                                                                               gcCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTgg                                                                 >  1:743215/1‑88 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTTATCCAGTGCTTCGCCAGTGATTTCCGCGATATGAACATAAGCATCGAAAGCCATGCTTACCTCTCCCTGGTAAAGGTTCGGCTGATAACGACGATGTTCAGACATCGTCATCCTGAGGTTGCTCCCGAAGTAAAACGGCGGCCGGTACCGACCGCCTTCAGGCTCAATACAACACTAACCCCTATGTCAGCCCTTGCAACAGCAAACGCCCTATCAAGGAATATTTGGGAAATCTCCTACAGATTCGACGATGATGACGATTTGGTCAGTTTCGATGCTGCTTAAAGCTTC  >  NC_002947/5554830‑5555123

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