Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,195,904 A→C 6.9% C284G (TGC→GGC)  rfaE ← fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,195,9040AC 6.9% 312.5 / 12.0 101C284G (TGC→GGC) rfaEfused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/48);  new base C (4/3);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 101 total aligned reads displayed.

GCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCT  >  NC_000913/3195756‑3196050
                                                                                                                                                    |                                                                                                                                                  
gCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAg                                                                                                                                                   <  1:1355301/150‑1 (MQ=255)
gCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAg                                                                                                                                                   >  2:250187/1‑150 (MQ=255)
  tACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGc                                                                                                                                                 <  2:1586689/150‑1 (MQ=255)
      gCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCCCGTCCCCGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCAGCATTGGCAAAGAAGCAGGcttct                                                                                                                                             >  1:821842/1‑150 (MQ=255)
      gCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCCTGcttct                                                                                                                                             >  1:360660/1‑150 (MQ=255)
      gCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGcttct                                                                                                                                             >  2:878199/1‑150 (MQ=255)
       ccAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGcttctt                                                                                                                                            >  1:1073813/1‑150 (MQ=255)
          gCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCa                                                                                                                                         <  2:37862/150‑1 (MQ=255)
           cTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCGAAGAAGCCTGcttctttgcc                                                                                                                                        <  2:1390391/150‑5 (MQ=255)
           cTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCCTGcttcttttcc                                                                                                                                        <  2:1234658/150‑5 (MQ=255)
           cTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCCTGcttctttgcc                                                                                                                                        <  2:360660/150‑5 (MQ=255)
            ttCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGc                                                                                                                                       <  1:225226/150‑1 (MQ=255)
            ttCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGc                                                                                                                                       <  2:1318036/150‑1 (MQ=255)
                gTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCTTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAAt                                                                                                                                   <  2:830113/150‑1 (MQ=255)
                gTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAAt                                                                                                                                   >  2:1254498/1‑150 (MQ=255)
                gTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCAATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAAt                                                                                                                                   <  2:1208064/150‑1 (MQ=255)
                  tCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTa                                                                                                                                 <  2:1048291/150‑1 (MQ=255)
                   ccTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACTACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTAc                                                                                                                                >  2:752413/1‑150 (MQ=255)
                       ttCGGTCATCACGCAAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGc                                                                                                                            <  1:92020/150‑1 (MQ=255)
                        tCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCt                                                                                                                           <  1:688880/150‑1 (MQ=255)
                              atcaCGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGc                                                                                                                     <  1:462966/150‑1 (MQ=255)
                               tcaCGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCg                                                                                                                    <  2:962938/150‑1 (MQ=255)
                                  cGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTg                                                                                                                 >  1:1099751/1‑150 (MQ=255)
                                   gCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTgc                                                                                                                >  1:1219207/1‑150 (MQ=255)
                                    ccAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTgcc                                                                                                               >  1:541773/1‑150 (MQ=255)
                                     cAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTgccg                                                                                                              <  1:108098/150‑1 (MQ=255)
                                        aGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCa                                                                                                           >  2:1010661/1‑150 (MQ=255)
                                        aGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCa                                                                                                           <  2:1028796/150‑1 (MQ=255)
                                            aaatgcTGTACGTGGACGTGCAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGAc                                                                                                       >  2:1030751/6‑150 (MQ=255)
                                             gTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACg                                                                                                      >  2:985570/1‑150 (MQ=255)
                                             gTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACg                                                                                                      >  2:89606/1‑150 (MQ=255)
                                               aTCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGcc                                                                                                    <  2:1219207/150‑1 (MQ=255)
                                               aTCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGcc                                                                                                    >  1:1059090/1‑150 (MQ=255)
                                               aTCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGcc                                                                                                    <  2:1012525/150‑1 (MQ=255)
                                                tCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGTCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCa                                                                                                   >  2:928668/1‑150 (MQ=255)
                                                tCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCa                                                                                                   >  2:922765/1‑150 (MQ=255)
                                                tCTGAACTTCAACGTACAACATTTTAAACCTCGATCGCACAAACCTAGGAGGTTCAAAGTTTTCAGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCa                                                                                                   >  2:1360862/1‑150 (MQ=255)
                                                   gCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATc                                                                                                <  2:741279/150‑1 (MQ=255)
                                                    cACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCa                                                                                               >  1:83174/1‑150 (MQ=255)
                                                     aCGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCAc                                                                                              >  1:1107550/1‑150 (MQ=255)
                                                       gTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCg                                                                                            <  2:1048482/150‑1 (MQ=255)
                                                        tCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCgt                                                                                           >  2:635211/1‑150 (MQ=255)
                                                        tCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCgt                                                                                           >  2:1661231/1‑150 (MQ=255)
                                                                cAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCg                                                                                   >  1:1138590/1‑150 (MQ=255)
                                                                 aGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGc                                                                                  >  1:955845/1‑150 (MQ=255)
                                                                 aGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGc                                                                                  <  2:1465669/150‑1 (MQ=255)
                                                                 aGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGc                                                                                  <  2:83174/150‑1 (MQ=255)
                                                                  gCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGcc                                                                                 >  1:1509472/1‑149 (MQ=255)
                                                                  gCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCa                                                                                 >  1:1280697/1‑150 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                               <  1:250187/150‑1 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                               <  2:1073813/150‑1 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                               >  1:136754/1‑150 (MQ=255)
                                                                      tttCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCgg                                                                             <  1:1082524/150‑1 (MQ=255)
                                                                         ccAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTaa                                                                          >  1:425518/1‑150 (MQ=255)
                                                                             gtCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTc                                                                      <  1:966571/149‑1 (MQ=255)
                                                                             cTCGATCGGCGAAACCGGCTAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCTAATTACCCGCTGCCAGCGTTGCCGCCACGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTc                                                                      <  1:444128/150‑1 (MQ=255)
                                                                               cGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCAt                                                                    <  2:1570918/150‑1 (MQ=255)
                                                                               cGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCAt                                                                    >  2:1540968/1‑150 (MQ=255)
                                                                               cGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCAt                                                                    >  2:873195/1‑150 (MQ=255)
                                                                                gATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATa                                                                   >  2:1668663/1‑150 (MQ=255)
                                                                                       gAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTcc                                                            <  1:1248191/150‑1 (MQ=255)
                                                                                        aaaCCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCt                                                           <  2:781055/150‑1 (MQ=255)
                                                                                           ccGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTgcg                                                        >  2:880894/1‑150 (MQ=255)
                                                                                               ggAGGTTCCCAGTTTGCTGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTg                                                    <  1:1302516/150‑1 (MQ=255)
                                                                                               ggAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTgcg                                                        <  2:691935/146‑1 (MQ=255)
                                                                                               ggAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTgcg                                                        >  1:691935/1‑146 (MQ=255)
                                                                                                       ccAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCa                                            <  2:1059090/150‑1 (MQ=255)
                                                                                                         aGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCata                                          >  1:863408/1‑150 (MQ=255)
                                                                                                            ttGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGc                                       <  1:878199/150‑1 (MQ=255)
                                                                                                             tGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCa                                      >  2:1547379/1‑150 (MQ=255)
                                                                                                               ccGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAgc                                    >  1:1175246/1‑150 (MQ=255)
                                                                                                                   ccaccaCGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAgcggcg                                <  2:1280697/150‑1 (MQ=255)
                                                                                                                   ccaccaCGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAgcggcg                                <  2:541773/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                               <  1:129704/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                               <  2:216804/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                               <  1:135728/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                               <  2:207894/150‑1 (MQ=255)
                                                                                                                           ccAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGCCGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTAccc                        >  2:1128193/1‑150 (MQ=255)
                                                                                                                           agcgccgccGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTAccc                        <  1:1192098/147‑1 (MQ=255)
                                                                                                                                cgccgcATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTg                   >  2:98480/1‑150 (MQ=255)
                                                                                                                                 gccgcATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgc                  >  1:1615685/1‑150 (MQ=255)
                                                                                                                                   cgcATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgcag                >  1:825790/1‑150 (MQ=255)
                                                                                                                                    gcATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgcagc               <  1:1360862/150‑1 (MQ=255)
                                                                                                                                     cATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgcagca              <  1:635211/150‑1 (MQ=255)
                                                                                                                                      aTTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGcagcag             >  1:174684/1‑150 (MQ=255)
                                                                                                                                       cgctggaagaagCCTGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCagcagc            >  1:1411158/7‑150 (MQ=255)
                                                                                                                                        tGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCg           <  1:182728/150‑1 (MQ=255)
                                                                                                                                        tGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCg           <  1:176868/150‑1 (MQ=255)
                                                                                                                                        tGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGCCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCg           >  1:268413/1‑150 (MQ=255)
                                                                                                                                         ctggaagaagCCTGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGa          >  1:210126/5‑150 (MQ=255)
                                                                                                                                         ctggaagaagCCTGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGa          >  2:210986/5‑150 (MQ=255)
                                                                                                                                          gCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGAc         <  1:89606/150‑1 (MQ=255)
                                                                                                                                          gCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGCACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGAATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGAc         >  2:709240/1‑150 (MQ=255)
                                                                                                                                           cAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACa        >  2:832021/1‑150 (MQ=255)
                                                                                                                                           cAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACa        >  1:1236100/1‑150 (MQ=255)
                                                                                                                                              agaagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgca                 <  1:460556/138‑1 (MQ=255)
                                                                                                                                              agaagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgca                 >  2:460556/1‑138 (MQ=255)
                                                                                                                                                aagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATAccc   <  2:855092/150‑1 (MQ=255)
                                                                                                                                                aagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATAccc   <  2:425518/150‑1 (MQ=255)
                                                                                                                                                 agCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGCCGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCt  >  1:1548490/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCT  >  NC_000913/3195756‑3196050

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