Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,471,732 T→A 5.6% intergenic (‑499/‑395) PP_3961 ← / → PP_3962 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,471,7320TA 5.6% 185.2 / 4.2 71intergenic (‑499/‑395)PP_3961/PP_3962transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base A (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-02

GGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCACCCAATTTTCAC  >  NC_002947/4471584‑4471879
                                                                                                                                                    |                                                                                                                                                   
ggATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCAtt                                                                                                                                                    <  1:775871/150‑1 (MQ=255)
 gATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTg                                                                                                                                                   >  2:463137/1‑150 (MQ=255)
 gATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTg                                                                                                                                                   <  1:828448/150‑1 (MQ=255)
 gATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTg                                                                                                                                                   >  2:828448/1‑150 (MQ=255)
  aTTACGTATCGCTATTGCTGTGGAAGTGAAAGTGGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAAAGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTgg                                                                                                                                                  >  2:444728/1‑150 (MQ=255)
    tACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGc                                                                                                                                                <  1:1181496/150‑1 (MQ=255)
       gAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCgg                                                                                                                                             <  2:240363/150‑1 (MQ=255)
              tATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTggg                                                                                                                                                 <  1:72400/139‑1 (MQ=255)
              tATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTggg                                                                                                                                                 >  2:72400/1‑139 (MQ=255)
                   ccGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCtt                                                                                                                                 >  2:984508/1‑150 (MQ=255)
                         aaGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAtt                                                                                                                           <  2:1196102/150‑1 (MQ=255)
                           gTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTg                                                                                                                         >  2:383416/1‑150 (MQ=255)
                           gTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGGTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTg                                                                                                                         >  2:441996/1‑150 (MQ=255)
                                  tCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCaa                                                                                                                  >  1:316103/1‑150 (MQ=255)
                                     gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAttt                                                                                                                          >  2:1065777/1‑139 (MQ=255)
                                     gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAttt                                                                                                                          <  1:1065777/139‑1 (MQ=255)
                                         aGCCACGGTAAACCACGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCaa                                                                                                           <  1:488346/150‑1 (MQ=255)
                                          gCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACAATTTGACACCAAGATTCAAt                                                                                                          >  1:639998/1‑150 (MQ=255)
                                            cACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGa                                                                                                        >  1:251404/1‑150 (MQ=255)
                                             aCGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGAc                                                                                                       >  1:118642/1‑150 (MQ=255)
                                                     ccAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCg                                                                                               >  1:1147535/1‑150 (MQ=255)
                                                       aaGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTc                                                                                             <  1:497834/150‑1 (MQ=255)
                                                       aaGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTc                                                                                             <  1:70155/150‑1 (MQ=255)
                                                           ccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCat                                                                                         >  2:167776/1‑150 (MQ=255)
                                                                 cTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCg                                                                                   >  1:1059599/1‑150 (MQ=255)
                                                                     cGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCa                                                                                                            <  1:420271/121‑1 (MQ=255)
                                                                     cGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCa                                                                                                            >  2:420271/1‑121 (MQ=255)
                                                                               ggAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGCCGCATATGTCGCAGATCACGCGCAt                                                                     >  1:1196245/1‑150 (MQ=255)
                                                                                    gTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATggg                                                                  >  2:128819/1‑148 (MQ=255)
                                                                                    gTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATggg                                                                  <  1:128819/148‑1 (MQ=255)
                                                                                     tGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTg                                                               <  1:463137/150‑1 (MQ=255)
                                                                                       gATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGc                                                             <  1:1101177/150‑1 (MQ=255)
                                                                                            gTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACGCCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                                                        >  1:1084771/1‑150 (MQ=255)
                                                                                              tAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCa                                                                                 <  2:230179/123‑1 (MQ=255)
                                                                                              tAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCa                                                                                 >  1:230179/1‑123 (MQ=255)
                                                                                                aGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTg                                                    >  1:787000/1‑150 (MQ=255)
                                                                                                aGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTg                                                    <  1:167776/150‑1 (MQ=255)
                                                                                                 gAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg                                                   <  2:118642/150‑1 (MQ=255)
                                                                                                  aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGc                                                  <  2:1147535/150‑1 (MQ=255)
                                                                                                  aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACAATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGc                                                  <  2:639998/150‑1 (MQ=255)
                                                                                                       aggaaaGCCCTGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGACACCCCGCATATGTGTCACAACGCGTGCAGGGTCTGGCCCTGCGTTGGCGAAcc                                             >  1:569658/4‑150 (MQ=255)
                                                                                                          aaaGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTg                                          >  1:1103260/1‑150 (MQ=255)
                                                                                                           aaGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGc                                         <  2:1059599/150‑1 (MQ=255)
                                                                                                            aGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                                                        >  1:1104165/1‑134 (MQ=255)
                                                                                                            aGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                                                        <  2:1104165/134‑1 (MQ=255)
                                                                                                              ccTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGAt                                      <  2:787000/150‑1 (MQ=255)
                                                                                                              ccTAGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAGTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGAt                                      <  1:924150/150‑1 (MQ=255)
                                                                                                                ttGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGa                                    >  1:380907/1‑150 (MQ=255)
                                                                                                                   aataatCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACAt                                 >  2:102245/1‑150 (MQ=255)
                                                                                                                    ataatCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACAtt                                <  2:143477/150‑1 (MQ=255)
                                                                                                                       atCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTg                                                    <  2:1188503/127‑1 (MQ=255)
                                                                                                                       atCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTg                                                    >  1:1188503/1‑127 (MQ=255)
                                                                                                                        tCAGGGCTTTCCTCATGCCTCTGGACCATTGGGCCGGCGCGCGCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAg                            >  1:26440/1‑150 (MQ=255)
                                                                                                                         cAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAgg                           <  1:441996/150‑1 (MQ=255)
                                                                                                                          aGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAggg                          >  1:90954/1‑150 (MQ=255)
                                                                                                                            ggCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGAc                        <  1:102245/150‑1 (MQ=255)
                                                                                                                              cTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACtt                      <  2:1103260/150‑1 (MQ=255)
                                                                                                                                  ccTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTtttct                  <  2:380907/150‑1 (MQ=255)
                                                                                                                                   cTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATGCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTCGACCTCCGTTGGCGCACCCTGCGATGACATTGTAGGGACGtttctt                 <  2:979312/150‑1 (MQ=255)
                                                                                                                                    tCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAAt                                                                                                          <  2:132857/60‑1 (MQ=255)
                                                                                                                                    tCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAAt                                                                                                          >  1:132857/1‑60 (MQ=255)
                                                                                                                                     cATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTtttctttc               <  2:251404/150‑1 (MQ=255)
                                                                                                                                       tGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCAc             <  1:383416/150‑1 (MQ=255)
                                                                                                                                        gCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCAcc            >  2:993801/1‑150 (MQ=255)
                                                                                                                                        gCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCAcc            >  1:117203/1‑150 (MQ=255)
                                                                                                                                         gcgccggCCCAATGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAgg                           <  2:1047116/129‑1 (MQ=255)
                                                                                                                                         gcgccggCCCAATGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAgg                           >  1:1047116/6‑134 (MQ=255)
                                                                                                                                           gccggCCCAATGGGCCAGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGc                                                                                           >  2:819949/4‑68 (MQ=255)
                                                                                                                                           gccggCCCAATGGGCCAGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGc                                                                                           <  1:819949/65‑1 (MQ=255)
                                                                                                                                            ctGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGCCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCACCCAAt        <  2:1196245/150‑1 (MQ=255)
                                                                                                                                                 ccATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCACCCAATTTTCa   <  2:316103/150‑1 (MQ=255)
                                                                                                                                                  cATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGAGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCACCCAATTTTCAc  >  1:698917/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACATTGTAGGGACTTTTCTTTCACCCAATTTTCAC  >  NC_002947/4471584‑4471879

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: