Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,696 C→T A221T (GCC→ACC)  yaaW ← UPF0174 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,6960CT97.4% 1394.6 / 15.6 428A221T (GCC→ACC) yaaWUPF0174 family protein
Reads supporting (aligned to +/- strand):  ref base C (3/8);  new base T (231/183);  total (235/193)
Fisher's exact test for biased strand distribution p-value = 7.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 428 total aligned reads displayed.

TGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTG  >  NC_000913/10548‑10843
                                                                                                                                                    |                                                                                                                                                   
tGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCGAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTg                                                                                                                                                   >  6:1706017/1‑151 (MQ=255)
tGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTg                                                                                                                                                   >  4:1383377/1‑151 (MQ=255)
 gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTgg                                                                                                                                                  >  5:1497938/1‑151 (MQ=255)
   aGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGa                                                                                                                                                 <  2:591927/150‑1 (MQ=255)
     tgtTTTCAGGGCGGGGCATAATAATCATCCAGTGTGTCAGTTTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTCCAGTACGGTTGGAAt                                                                                                                                               <  1:1382076/150‑1 (MQ=255)
     tgtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATc                                                                                                                                              >  4:106831/1‑151 (MQ=255)
     ggtTTTCAGGGCGCGGCATAATCATCAGCCAGTGCTGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATc                                                                                                                                              <  4:1184608/150‑1 (MQ=255)
      gtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCg                                                                                                                                             >  1:568099/1‑151 (MQ=255)
        tttCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGt                                                                                                                                            <  8:436552/150‑1 (MQ=255)
         ttCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTc                                                                                                                                           >  8:297083/1‑150 (MQ=255)
          tCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCa                                                                                                                                          <  2:229836/150‑1 (MQ=255)
           cAGGGCGCGGCATAATAATGAGCCAGTGGCGCTGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCAcg                                                                                                                                        <  4:952906/151‑1 (MQ=255)
           cAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCAcg                                                                                                                                        <  2:1299607/151‑1 (MQ=255)
             gggCGCTGCATAATAATCTGGCAGTGGGGCAGTGTCTTCGATCTTTTGAGGGGAACATGAATATTTTCCCCGGTTTGCGGTTTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCAcgcg                                                                                                                                      <  1:416287/151‑1 (MQ=255)
              ggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCAcgcg                                                                                                                                      >  6:569932/1‑150 (MQ=255)
                cgcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGAt                                                                                                                                    >  5:843935/1‑150 (MQ=255)
                cgcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATa                                                                                                                                   <  4:735431/151‑1 (MQ=255)
                cgcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATa                                                                                                                                   <  8:819730/151‑1 (MQ=255)
                  cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAgg                                                                                                                                 <  8:1208384/151‑1 (MQ=255)
                     cATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGgcgc                                                                                                                              >  5:26025/1‑151 (MQ=255)
                      ataataATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCt                                                                                                                             >  8:132683/1‑151 (MQ=255)
                      ataataATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCt                                                                                                                             >  8:923302/1‑151 (MQ=255)
                         ataatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTACCCATCCGGCGCAGGCAGGCGAGTTGCAGTACGGTTGGAATCGTCACGCGCTAGGCGCTGcc                                                                                                                          >  2:1175787/1‑151 (MQ=255)
                         ataatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGcc                                                                                                                          <  2:157765/151‑1 (MQ=255)
                           aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGc                                                                                                                        >  6:1257427/1‑151 (MQ=255)
                            atCAGCCAGTGGGGCAGTCTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCt                                                                                                                       <  1:1003016/151‑1 (MQ=255)
                               aGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGAc                                                                                                                    >  3:307616/1‑151 (MQ=255)
                                gCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGAcc                                                                                                                   <  4:1293297/151‑1 (MQ=255)
                                 ccAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCg                                                                                                                  >  7:944059/1‑151 (MQ=255)
                                       gggCAGTGTCTACGATCTTTTGA‑GGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTatc                                                                                                           <  5:558618/151‑3 (MQ=255)
                                           aGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACGATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcccc                                                                                                        <  6:1097996/151‑1 (MQ=255)
                                              gtCTACGATGTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAGCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTAGGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAtt                                                                                                     <  1:1215025/151‑1 (MQ=255)
                                              gtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAtt                                                                                                     <  4:173702/151‑1 (MQ=255)
                                                cTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                    >  8:1131883/1‑150 (MQ=255)
                                                 tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                    >  3:282729/1‑149 (MQ=255)
                                                 tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                    >  1:859520/1‑149 (MQ=255)
                                                  aCGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                    <  7:540266/148‑1 (MQ=255)
                                                  aCGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                    <  8:1697125/148‑1 (MQ=255)
                                                   cGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTg                                                                                                >  6:1040351/1‑151 (MQ=255)
                                                       cTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAATCCCATTTAGTGCCGc                                                                                            >  8:1020243/1‑151 (MQ=255)
                                                       cTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                                                                                            >  2:956225/1‑151 (MQ=255)
                                                        ttttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                                                                                            >  2:1120697/1‑150 (MQ=255)
                                                         tttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGc                                                                                               <  1:263280/146‑1 (MQ=255)
                                                         tttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGc                                                                                            <  6:289978/149‑1 (MQ=255)
                                                           tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCt                                                                                        >  5:1148174/1‑151 (MQ=255)
                                                           tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCt                                                                                        >  2:1462202/1‑151 (MQ=255)
                                                             aGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGGCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                      >  3:139506/1‑151 (MQ=255)
                                                             aGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                      >  2:866237/1‑151 (MQ=255)
                                                                  aaaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGc                                                                                 <  5:1087811/151‑1 (MQ=255)
                                                                     atgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTc                                                                              <  7:454762/151‑1 (MQ=255)
                                                                     atgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTc                                                                              >  3:783004/1‑151 (MQ=255)
                                                                     atgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTc                                                                              <  2:457088/151‑1 (MQ=255)
                                                                        aaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCa                                                                           <  1:750642/151‑1 (MQ=255)
                                                                          aatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGc                                                                         >  8:1470666/1‑151 (MQ=255)
                                                                           atTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGCCCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGcc                                                                        >  4:304880/1‑151 (MQ=255)
                                                                           atTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGcc                                                                        >  2:727505/1‑151 (MQ=255)
                                                                             tttCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGcccc                                                                      >  7:192479/1‑151 (MQ=255)
                                                                              ttCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTACCCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccg                                                                     >  3:1630011/1‑151 (MQ=255)
                                                                               tCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                                                                    >  4:604924/1‑151 (MQ=255)
                                                                               tCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccg                                                                     <  6:1472771/150‑1 (MQ=255)
                                                                                ccccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcg                                                                   <  1:1476937/151‑1 (MQ=255)
                                                                                ccccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                                                                    <  1:682159/150‑1 (MQ=255)
                                                                                    ggTTTCCGCTATCAGCCCTGAGTGGCGCTAACCATCCGGCGCAGCCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAATTAGTGCCGCACCTACAGCGCCTCCCACCCccgcgccgc                                                               >  4:425468/1‑151 (MQ=255)
                                                                                     gTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCGGCCGCAGGCAGGCGATTTGGAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACGGCTTTACGCCCATTTAGTGCCGCTCCTACAGGGCCTCCCAGCCccgcgccgcc                                                              >  4:233510/1‑150 (MQ=255)
                                                                                     gTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgccgcg                                                              >  6:728170/1‑151 (MQ=255)
                                                                                       ttCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCcgcgccgcgc                                                             >  3:899436/1‑150 (MQ=255)
                                                                                          cGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCgcagc                                                          >  5:1627731/1‑150 (MQ=255)
                                                                                           ggTATCAGACCTGAGTGGCGCTAACCATCGGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCgcagca                                                         >  6:826633/1‑150 (MQ=255)
                                                                                           ggTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCgcagca                                                         >  2:1006475/1‑150 (MQ=255)
                                                                                            gTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCaaa                                                       <  8:203958/151‑1 (MQ=255)
                                                                                               tCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAAcc                                                     >  5:1430417/1‑150 (MQ=255)
                                                                                                 aGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATg                                                  >  3:1507720/1‑151 (MQ=255)
                                                                                                  gACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGc                                                 >  7:107311/1‑151 (MQ=255)
                                                                                                  gACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCACTACGGTTGCAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACACGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGc                                                 >  5:1266195/1‑151 (MQ=255)
                                                                                                       gAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCACACCATGCCCa                                              >  8:759875/1‑149 (MQ=255)
                                                                                                       gAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCaa                                             <  7:1234180/150‑1 (MQ=255)
                                                                                                        cgTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGt                                           <  1:533376/150‑1 (MQ=255)
                                                                                                             cgcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCt                                      >  1:1662263/1‑151 (MQ=255)
                                                                                                              gcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTc                                     >  6:1591227/1‑151 (MQ=255)
                                                                                                              gcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTc                                     >  4:15987/1‑151 (MQ=255)
                                                                                                                 aaCCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCAtt                                  <  7:1358047/151‑1 (MQ=255)
                                                                                                                 aaCCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCAtt                                  >  7:1455293/1‑151 (MQ=255)
                                                                                                                   ccATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTACCCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATtgc                                >  6:661006/1‑151 (MQ=255)
                                                                                                                         ggCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTg                          >  1:1615208/1‑151 (MQ=255)
                                                                                                                         ggCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTg                          >  2:324784/1‑151 (MQ=255)
                                                                                                                          gcgcAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTgg                         <  6:1148877/151‑1 (MQ=255)
                                                                                                                           cgcAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTggg                        <  4:921166/151‑1 (MQ=255)
                                                                                                                            gcaggcagGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGt                       <  1:1110503/151‑1 (MQ=255)
                                                                                                                              aggcaggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGc                     >  3:19446/1‑151 (MQ=255)
                                                                                                                                  aggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTa                  >  7:1157116/1‑150 (MQ=255)
                                                                                                                                     cGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAAt              >  7:1531792/1‑151 (MQ=255)
                                                                                                                                      gATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATg             >  7:335475/1‑151 (MQ=255)
                                                                                                                                       aTTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCTGCAAACTATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGc            >  7:656608/1‑151 (MQ=255)
                                                                                                                                       aTTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGc            >  3:1027953/1‑151 (MQ=255)
                                                                                                                                        tttGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCg           >  5:95861/1‑151 (MQ=255)
                                                                                                                                           gCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCggg         <  4:1441452/150‑1 (MQ=255)
                                                                                                                                           gCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGt        <  3:1009676/151‑1 (MQ=255)
                                                                                                                                            cAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTc       <  3:1261483/151‑1 (MQ=255)
                                                                                                                                             aGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCa      <  5:1338420/151‑1 (MQ=255)
                                                                                                                                               tACGGTCGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGt    >  6:972038/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTG  >  NC_000913/10548‑10843

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