Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,195,905 G→T 7.0% A283A (GCC→GCA 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,9050GT 7.0% 335.8 / 8.6 99A283A (GCC→GCArfaEfused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base G (45/47);  new base T (4/3);  total (49/50)
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.30e-01

Only 100 of 104 total aligned reads displayed.

GCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTC  >  NC_000913/3195756‑3196054
                                                                                                                                                     |                                                                                                                                                     
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)
            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:1028796/150‑1 (MQ=255)
                                        aGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCa                                                                                                               >  2:1010661/1‑150 (MQ=255)
                                             gTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACg                                                                                                          >  2:89606/1‑150 (MQ=255)
                                             gTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACg                                                                                                          >  2:985570/1‑150 (MQ=255)
                                               aTCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGcc                                                                                                        >  1:1059090/1‑150 (MQ=255)
                                               aTCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGcc                                                                                                        <  2:1219207/150‑1 (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                                                                                      <  2:1465669/150‑1 (MQ=255)
                                                                 aGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGc                                                                                      <  2:83174/150‑1 (MQ=255)
                                                                 aGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGc                                                                                      >  1:955845/1‑150 (MQ=255)
                                                                  gCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGcc                                                                                     >  1:1509472/1‑149 (MQ=255)
                                                                  gCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCa                                                                                     >  1:1280697/1‑150 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                                   <  1:250187/150‑1 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                                   >  1:136754/1‑150 (MQ=255)
                                                                    aTTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCAcc                                                                                   <  2:1073813/150‑1 (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:1540968/1‑150 (MQ=255)
                                                                               cGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCAt                                                                        >  2:873195/1‑150 (MQ=255)
                                                                               cGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCAt                                                                        <  2:1570918/150‑1 (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                                                            >  1:691935/1‑146 (MQ=255)
                                                                                               ggAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTgcg                                                            <  2:691935/146‑1 (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:207894/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                                   <  2:216804/150‑1 (MQ=255)
                                                                                                                    caccacGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGgcgc                                   <  1:135728/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:176868/150‑1 (MQ=255)
                                                                                                                                        tGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGCCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCg               >  1:268413/1‑150 (MQ=255)
                                                                                                                                         ctggaagaagCCTGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGa              >  2:210986/5‑150 (MQ=255)
                                                                                                                                         ctggaagaagCCTGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGa              >  1:210126/5‑150 (MQ=255)
                                                                                                                                          gCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGAc             <  1:89606/150‑1 (MQ=255)
                                                                                                                                          gCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGCACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGAATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGAc             >  2:709240/1‑150 (MQ=255)
                                                                                                                                           cAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACa            >  1:1236100/1‑150 (MQ=255)
                                                                                                                                           cAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACa            >  2:832021/1‑150 (MQ=255)
                                                                                                                                              agaagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgca                     >  2:460556/1‑138 (MQ=255)
                                                                                                                                              agaagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTgca                     <  1:460556/138‑1 (MQ=255)
                                                                                                                                                aagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATAccc       <  2:425518/150‑1 (MQ=255)
                                                                                                                                                aagCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATAccc       <  2:855092/150‑1 (MQ=255)
                                                                                                                                                 agCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGCCGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCt      >  1:1548490/1‑150 (MQ=255)
                                                                                                                                                     ggCTTCTTCCAGCGAATTACCCGCTGTCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCACACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTc  >  2:903090/1‑150 (MQ=255)
                                                                                                                                                     ggCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTc  <  2:955845/150‑1 (MQ=255)
                                                                                                                                                     ggCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTc  >  2:1038101/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCTTCTTCCAGCGAATTACCCGCTGCCAGCGTTGCCGCCAGGACGCCAATCACCGTGTCGCCCGCACCGGTAACGTCATACACTTCCTGCGCTTGGGTTGGCATATGCAGCGGCGCTTTACCCGGTTGCAGCAGCGACATACCCTGTTC  >  NC_000913/3195756‑3196054

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: