Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,529,039 C→T G402D (GGT→GAT)  sgcC ← putative PTS system EIIC permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,0390CT94.7% 1506.9 / 50.6 470G402D (GGT→GAT) sgcCputative PTS system EIIC permease component
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base T (259/185);  total (273/197)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 471 total aligned reads displayed.

CCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCA  >  NC_000913/4528891‑4529186
                                                                                                                                                    |                                                                                                                                                   
 ccAATAACCTCTTTTAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGa                                                                                                                                                  <  6:376811/151‑1 (MQ=255)
 ccAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGa                                                                                                                                                  <  8:1587597/151‑1 (MQ=255)
   aaTAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAc                                                                                                                                                <  8:1436478/151‑1 (MQ=255)
    aTAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAAc                                                                                                                                                >  6:914929/1‑150 (MQ=255)
      aaCCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACat                                                                                                                                              <  6:1658531/150‑1 (MQ=255)
        ccTCTTTCAGCCAACTCATTGAAAGTACTCCTTTTTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGATCATAtt                                                                                                                                           >  1:392315/1‑151 (MQ=255)
        ccTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAtt                                                                                                                                           >  8:770511/1‑151 (MQ=255)
        ccTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAtt                                                                                                                                           >  3:50359/1‑151 (MQ=255)
        ccTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATAtt                                                                                                                                           >  3:1294807/1‑151 (MQ=255)
         ctctTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTg                                                                                                                                          >  6:514091/1‑151 (MQ=255)
            tttCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGCCAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGacc                                                                                                                                       >  5:1203079/1‑151 (MQ=255)
             ttCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGCCGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGccc                                                                                                                                       >  5:690566/1‑150 (MQ=255)
             ttCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGccc                                                                                                                                       <  1:511972/150‑1 (MQ=255)
              tCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTAAGAAATCGACCAATCGAACATATTGCCCAc                                                                                                                                     >  2:211388/1‑151 (MQ=255)
                   cAACTCATTGAAAGTACTCCTTATTCCGATATTTCACGTTTACGTAAGTACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCGACggcc                                                                                                                                 >  4:282066/1‑148 (MQ=255)
                        cATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCt                                                                                                                           >  8:842360/1‑151 (MQ=255)
                         aTTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGACAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                          >  7:1597261/1‑151 (MQ=255)
                           tGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                          >  1:829045/1‑149 (MQ=255)
                            gAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTgag                                                                                                                        <  6:1448773/150‑1 (MQ=255)
                            gAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAAATCAGAAATCGACCAATCGAGCATATTGCCCACGGACAACGCTgaga                                                                                                                       <  1:607282/151‑1 (MQ=255)
                             aaaGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTg                                                                                                                          >  3:151655/1‑147 (MQ=255)
                             aaaGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACCACAACGCCAATAATCCCAAGCGCCATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACACCGCTg                                                                                                                          >  6:62150/1‑147 (MQ=255)
                              aaGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAtt                                                                                                                     <  3:1168125/151‑1 (MQ=255)
                               aGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAttt                                                                                                                    >  5:278950/1‑151 (MQ=255)
                               aGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAttt                                                                                                                    >  8:1688787/1‑151 (MQ=255)
                                  aCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTgcg                                                                                                                 <  3:597242/151‑1 (MQ=255)
                                  aCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTgcg                                                                                                                 >  2:306117/1‑151 (MQ=255)
                                  aCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTgc                                                                                                                  <  5:780991/150‑1 (MQ=255)
                                     ccTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCAc                                                                                                              <  2:538241/151‑1 (MQ=255)
                                     ccTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCAc                                                                                                              <  1:1266481/151‑1 (MQ=255)
                                        tATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACttt                                                                                                           >  8:359025/1‑151 (MQ=255)
                                         aTTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACttt                                                                                                           <  2:469187/150‑1 (MQ=255)
                                          ttCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACCACAACGAAAATAATCACAAGCGACATCAATTCAGAAATAGACAAATCGAAAATATTGGCCACGGACAACGCTGAGATTTGATCACTTTc                                                                                                          >  6:804856/1‑150 (MQ=255)
                                           tCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCAc                                                                                                        >  7:1107775/1‑151 (MQ=255)
                                             cGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCg                                                                                                      <  4:1411638/151‑1 (MQ=255)
                                             cGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCg                                                                                                      >  5:833944/1‑151 (MQ=255)
                                              gATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCgg                                                                                                     >  3:417306/1‑151 (MQ=255)
                                                tAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGc                                                                                                    >  6:856446/1‑150 (MQ=255)
                                                 aaTTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCa                                                                                                   >  5:474490/1‑150 (MQ=255)
                                                  aTTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCaa                                                                                                  >  6:981286/1‑150 (MQ=255)
                                                      aCGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTg                                                                                             <  1:1099790/151‑1 (MQ=255)
                                                         tttACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGaa                                                                                           <  3:1004104/150‑1 (MQ=255)
                                                           tACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAAc                                                                                         <  1:981375/150‑1 (MQ=255)
                                                               tAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCt                                                                                    >  6:1064879/1‑151 (MQ=255)
                                                                aaGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACCACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCtt                                                                                   >  4:358018/1‑151 (MQ=255)
                                                                aaGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCtt                                                                                   >  2:399478/1‑151 (MQ=255)
                                                                aaGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCt                                                                                    <  5:609002/150‑1 (MQ=255)
                                                                    aCAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAg                                                                               <  8:677641/151‑1 (MQ=255)
                                                                     cAGCCCGACGCTGGAGACAATACCGACCGCCACGACAACGCCAATAATCCGACCCGACATCAATTCAGAAATCGACCAATCTAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCTCCTTTGAgc                                                                              <  7:297414/151‑1 (MQ=255)
                                                                     cAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAgc                                                                              >  5:1320016/1‑151 (MQ=255)
                                                                      aGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAgcg                                                                             <  8:254080/151‑1 (MQ=255)
                                                                       gCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAgcgc                                                                            >  5:1488237/1‑151 (MQ=255)
                                                                         acgacgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGcc                                                                           >  5:967160/1‑150 (MQ=255)
                                                                          cgacgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCAt                                                                         <  6:226990/151‑1 (MQ=255)
                                                                           gagtcTGGCGTCAATACCGACCTCCACGACAAGTCGAATAATCCCAAGCGTCATCAATTCTGAATTCGTCCATTCGTTCATATTGCCCAATGACAAGGCTGAGATTTGCGCACTTTCACCTGCAAAACTGAAACCGGCTTTGAGGGCGAtt                                                                        <  1:1537769/147‑1 (MQ=255)
                                                                             cgCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTc                                                                      <  7:504937/151‑1 (MQ=255)
                                                                               cTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCAGAGAATTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTc                                                                      >  2:649228/1‑149 (MQ=255)
                                                                                tGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGt                                                                   >  5:1132902/1‑151 (MQ=255)
                                                                                  gCGACAATACCGACCGCCACGACAACGCCAATAATCCGAAGCGACCTAAATTCAGAAATCGACCAATCTAACATATTGCCAACGGACAACGCTGAGATTTGCGCACTTACACCGGCAAAACTGAAACCGCCTTTGAGAGCCATTTCAGTg                                                                  >  8:552722/1‑150 (MQ=255)
                                                                                   cGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTc                                                                      >  5:1421552/1‑145 (MQ=255)
                                                                                    gACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCACTCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACCGAAACCGCCTTTGAGCGCCATTTCAGTGAAg                                                               <  2:1278075/151‑1 (MQ=255)
                                                                                     aCAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTTAGATTAGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGt                                                              <  2:1413798/151‑1 (MQ=255)
                                                                                     aCAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGt                                                              >  6:263218/1‑151 (MQ=255)
                                                                                      cAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGt                                                              >  4:1617325/1‑150 (MQ=255)
                                                                                       aaTACCGACCGCCACGACAACGCCAATAATACCAAGCGACATCAAATCAGAAATCGACGAAACGAACATATAGCCCACCCACTACGCTGAGATTTGCGCAGATTCACCCGCAAATCTGAACCCGCCATTGAGCGCCATTTCACTGAAgaa                                                             >  4:786955/1‑148 (MQ=255)
                                                                                          accgaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGgc                                                         <  1:692096/151‑1 (MQ=255)
                                                                                           ccgaccgCCACGACAAAGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGgcg                                                        <  5:1591860/151‑1 (MQ=255)
                                                                                               ccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGaa                                                    <  8:500004/151‑1 (MQ=255)
                                                                                                cgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAAc                                                   <  3:982432/151‑1 (MQ=255)
                                                                                                     cgtcAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGATCATATTGGCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAATACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGt                                              <  6:590446/148‑1 (MQ=255)
                                                                                                     cGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGt                                              <  7:1373849/151‑1 (MQ=255)
                                                                                                     cGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGt                                              <  4:726209/151‑1 (MQ=255)
                                                                                                       aCAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTgg                                            <  3:791076/151‑1 (MQ=255)
                                                                                                         aaCGCCAATAATGCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGAGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCg                                          <  1:860791/151‑1 (MQ=255)
                                                                                                         aaCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCg                                          <  7:823618/151‑1 (MQ=255)
                                                                                                           cGCCAATACTCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCACGAACTGGGTGGCGAt                                        <  7:691081/151‑1 (MQ=255)
                                                                                                           cGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGAt                                        <  3:205163/151‑1 (MQ=255)
                                                                                                               aataatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATcagc                                    >  7:405200/1‑151 (MQ=255)
                                                                                                                  aatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                 <  1:1026207/151‑1 (MQ=255)
                                                                                                                  aatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                 <  2:770821/151‑1 (MQ=255)
                                                                                                                  aatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAgg                                 >  6:122913/1‑151 (MQ=255)
                                                                                                                    tCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGAc                               >  3:747678/1‑151 (MQ=255)
                                                                                                                    tCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGAc                               >  4:1207127/1‑151 (MQ=255)
                                                                                                                     cccAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACg                              <  3:1369109/151‑1 (MQ=255)
                                                                                                                            gACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatca                       <  7:893798/151‑1 (MQ=255)
                                                                                                                            gACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatc                        >  2:1542903/1‑150 (MQ=255)
                                                                                                                             aCATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatcat                      <  7:1518267/151‑1 (MQ=255)
                                                                                                                              cATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTCGCGATCAGCAGGACGGTAatcatc                     >  3:1416943/1‑151 (MQ=255)
                                                                                                                                  aaTTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACg                  >  1:1457589/1‑150 (MQ=255)
                                                                                                                                  aaTTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCAAACTGGGTGGCGATCAGCAGGACGGTAATCATCACg                  >  7:959306/1‑150 (MQ=255)
                                                                                                                                   aTTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGAt                >  4:548030/1‑151 (MQ=255)
                                                                                                                                     tCAGAAATCGACCAATCGAACTTATTGCCCACGGACATCGCTGAGATTTGCGCACTTTCACAGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCATTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg               <  5:1067427/150‑1 (MQ=255)
                                                                                                                                     tCAGAAATCGACCAATCGAACATATTCCCCACCGACAACGATGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGAGCAGCAGGACGGTAATCATCACGATg               <  6:1339156/150‑1 (MQ=255)
                                                                                                                                          aaTCGACCGATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCg         <  8:605365/151‑1 (MQ=255)
                                                                                                                                            tCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCt       >  1:241240/1‑151 (MQ=255)
                                                                                                                                            tCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCACGc        >  2:1333438/1‑150 (MQ=255)
                                                                                                                                              gACCAATCGAACATATTGCCCACGGACAACTCTGAGATTTGCGCACTTTGACCGGCAAAACTGAAACCTCCTTTTAGCTAAATTTCAGTGAAGTACGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTa     <  8:1236781/151‑1 (MQ=255)
                                                                                                                                                ccAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATc   >  6:328470/1‑151 (MQ=255)
                                                                                                                                                ccAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATc   >  4:995112/1‑151 (MQ=255)
                                                                                                                                                 cAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCa  <  4:1111305/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCA  >  NC_000913/4528891‑4529186

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