Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,166,155 G→C 5.1% Q522E (CAG→GAG) ‡ ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,1550GC 5.1% 368.8 / 6.7 112Q522E (CAG→GAG) ‡ygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (60/46);  new base C (4/2);  total (64/48)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCTGCG  >  NC_000913/3166006‑3166304
                                                                                                                                                     |                                                                                                                                                      
tCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTg                                                                                                                                                        >  2:64424/1‑150 (MQ=255)
 cAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCTAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGCCGTGCTCCGACGGCGATAAAGGCATCAATCAGCCTTGAAATACAACTCATTGCTGTAGACATAATCCTgg                                                                                                                                                       >  2:637551/1‑150 (MQ=255)
 cAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                       >  2:65435/1‑150 (MQ=255)
  aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTggg                                                                                                                                                      <  1:1671072/150‑1 (MQ=255)
   aaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                                                                                                                                                     <  2:17419/150‑1 (MQ=255)
   aaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCAAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                                                                                                                                                     <  2:1578688/150‑1 (MQ=255)
          aaTATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCATTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGc                                                                                                                                              <  2:1303778/150‑1 (MQ=255)
          aaTATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGc                                                                                                                                              <  1:706278/150‑1 (MQ=255)
            tatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAg                                                                                                                                            <  2:886758/150‑1 (MQ=255)
            tatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCCCAGGATTATgtct                                                                                                                                            <  1:1703602/150‑4 (MQ=255)
             atTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGc                                                                                                                                           <  1:457593/150‑1 (MQ=255)
                gCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGa                                                                                                                                        <  1:1261801/150‑1 (MQ=255)
                 ccGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa                                                                                                                                       <  1:220016/150‑1 (MQ=255)
                              aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                          >  1:426991/1‑150 (MQ=255)
                              aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                          <  2:752474/150‑1 (MQ=255)
                              aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                          >  1:616693/1‑150 (MQ=255)
                               ttATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGCCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAg                                                                                                                        >  2:62343/1‑150 (MQ=255)
                                 aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAggg                                                                                                                      <  1:554648/150‑1 (MQ=255)
                                  tCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACACCTCTTTGCTGTAGACATAATCCTGGGGATTATGAAGCGGAAAACCGCCAACGTAGGGt                                                                                                                     >  1:1106921/1‑150 (MQ=255)
                                   cTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGtt                                                                                                                    >  1:1483931/1‑150 (MQ=255)
                                   cTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGtt                                                                                                                    >  1:204125/1‑150 (MQ=255)
                                    tGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGttt                                                                                                                   <  2:89156/150‑1 (MQ=255)
                                     tccGCTTCGAGTCACAGGATATCAACGTGAATAGCGCCGCCGGACTTGCTCCGACGGCGATAAAGGCATAAATTTGCCTTGATATACAACTATTTGCTGTAGAAATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTa                                                                                                                  <  1:1047353/149‑1 (MQ=255)
                                        gcgTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGc                                                                                                               <  1:342501/150‑1 (MQ=255)
                                            cGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCCGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                           >  2:1085788/1‑150 (MQ=255)
                                            cGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGtt                                                                                                           <  2:383928/150‑1 (MQ=255)
                                             gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGttt                                                                                                          <  2:1062838/150‑1 (MQ=255)
                                             gAGGCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGttt                                                                                                          <  2:1489489/150‑1 (MQ=255)
                                              aGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTg                                                                                                         <  1:1613612/150‑1 (MQ=255)
                                               gTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGa                                                                                                        >  1:1165833/1‑150 (MQ=255)
                                                tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGAt                                                                                                       <  1:1170119/150‑1 (MQ=255)
                                                 cacaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGAtt                                                                                                      >  2:362456/1‑150 (MQ=255)
                                                 cacaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGAtt                                                                                                      >  1:511560/1‑150 (MQ=255)
                                                        gaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCt                                                                                               <  1:655653/150‑1 (MQ=255)
                                                         ctCAACGTGTATTCCGCCTCCGGAAGTTCTCCTACGTAGATAAAGGCATAAATGTGCCTTGATATCCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTg                                                                                              <  2:1140021/149‑1 (MQ=255)
                                                           cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAt                                                                                            >  2:843376/1‑150 (MQ=255)
                                                            aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCCGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAta                                                                                           >  2:375624/1‑150 (MQ=255)
                                                            aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAta                                                                                           >  1:955488/1‑150 (MQ=255)
                                                             aCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAtag                                                                                          >  1:1143838/1‑150 (MQ=255)
                                                             aCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAtag                                                                                          <  2:1724841/150‑1 (MQ=255)
                                                                tGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGAtagtag                                                                                       <  2:1725274/150‑1 (MQ=255)
                                                                tGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGGAGGGTTTAAGCAGTTTGATTAACGGCTGAtagtag                                                                                       >  1:671470/1‑150 (MQ=255)
                                                                 gAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGa                                                                                      >  2:1115026/1‑150 (MQ=255)
                                                                       cgccgcCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACAAAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGaa                                                                                <  1:445251/150‑1 (MQ=255)
                                                                        gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAAt                                                                               >  2:126634/1‑150 (MQ=255)
                                                                           gccgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAg                                                                            >  1:818365/1‑150 (MQ=255)
                                                                             cgGACGCGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGAGAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                          >  2:871429/1‑150 (MQ=255)
                                                                                 cGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGc                                                                      >  2:195742/1‑150 (MQ=255)
                                                                                  gTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                     >  2:922140/1‑150 (MQ=255)
                                                                                  gTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                     <  1:843376/150‑1 (MQ=255)
                                                                                   tGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCt                                                                    >  1:602944/1‑150 (MQ=255)
                                                                                    gCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTg                                                                   >  1:1535242/1‑150 (MQ=255)
                                                                                      tCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgtt                                                                 >  2:164859/1‑150 (MQ=255)
                                                                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgtt                                                                 >  1:1494290/1‑150 (MQ=255)
                                                                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttg                                                                <  1:1641375/150‑1 (MQ=255)
                                                                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttg                                                                <  2:25011/150‑1 (MQ=255)
                                                                                        cGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttgg                                                               <  1:1635499/150‑1 (MQ=255)
                                                                                         gACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggt                                                              >  2:1638483/1‑150 (MQ=255)
                                                                                           cGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggttg                                                            >  1:1209378/1‑150 (MQ=255)
                                                                                             gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                          <  2:798851/150‑1 (MQ=255)
                                                                                             gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                          <  1:86667/150‑1 (MQ=255)
                                                                                             gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGAAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                          <  2:75662/150‑1 (MQ=255)
                                                                                              cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCCAGTTAGGGTTTAAGCAGTTTGATTAACGTCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                         >  1:127160/1‑150 (MQ=255)
                                                                                              cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                         <  2:489251/150‑1 (MQ=255)
                                                                                              cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                         >  1:874461/1‑150 (MQ=255)
                                                                                                  aaaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCa                                                     <  1:661510/150‑1 (MQ=255)
                                                                                                  aaaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCa                                                     <  1:1365587/150‑1 (MQ=255)
                                                                                                       cATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAATCCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTTATAGTAGATAGGAATCAGCGGTTCCTGTTGGTTGATGTTCACTTct                                                >  2:78531/1‑150 (MQ=255)
                                                                                                          cAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcc                                             >  1:560595/1‑150 (MQ=255)
                                                                                                          cAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcc                                             <  2:813743/150‑1 (MQ=255)
                                                                                                            aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAGTTTGCAGCGGAAAACCGCCAACGTAGGGTGTAAGGAGTTTGAGTACCGGCTGATTGTAGATAGGAATCAGAGGTGGATGTTGGGTGATGATCCCTTctgcctg                                           >  2:272860/1‑150 (MQ=255)
                                                                                                                gCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGa                                       <  1:1479119/150‑1 (MQ=255)
                                                                                                                 ccTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                      >  1:738577/1‑150 (MQ=255)
                                                                                                                 ccTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                      >  2:1506872/1‑150 (MQ=255)
                                                                                                                  cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATa                                     >  1:1441503/1‑150 (MQ=255)
                                                                                                                     gATATACAACTCTTTGCTGTAGGCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACaa                                  <  2:1005541/150‑1 (MQ=255)
                                                                                                                     gATATACAACTCTTTGCTGTAGACATAATCCTGGGGAGTATGCAGCGGAAAACCGCAAACGTAGGGTTTAAGCAGTTTGATTAAGGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACaa                                  >  2:669680/1‑150 (MQ=255)
                                                                                                                       tataCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACg                                >  1:338389/1‑150 (MQ=255)
                                                                                                                       tataCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACg                                >  2:654031/1‑150 (MQ=255)
                                                                                                                           cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTACCGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCAGATACAACGCAtt                            >  2:299667/1‑150 (MQ=255)
                                                                                                                           cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCAtt                            >  2:738384/1‑150 (MQ=255)
                                                                                                                               tctTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGGTTAAGCCGTTTGAGTAACGGCTGAGTGTAGAGAGGACTCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGc                        >  2:172056/1‑150 (MQ=255)
                                                                                                                                ctTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                       >  1:823566/1‑150 (MQ=255)
                                                                                                                                ctTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                       >  2:635099/1‑150 (MQ=255)
                                                                                                                                 tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCtt                      <  2:204125/150‑1 (MQ=255)
                                                                                                                                 tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCtt                      <  2:1449604/150‑1 (MQ=255)
                                                                                                                                  ttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCttt                     <  2:426991/150‑1 (MQ=255)
                                                                                                                                     cTGTAGACATAATCCTGGTGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTc                  <  1:458476/150‑1 (MQ=255)
                                                                                                                                       gTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGc                >  2:1180671/1‑150 (MQ=255)
                                                                                                                                       gTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGACAACGTAGGCTTTAAGCAGTTTGAATAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAAAGCATTACGCTTTGTCGc                >  2:1214736/1‑150 (MQ=255)
                                                                                                                                       cgctgcATAATCCCCAGGATTATGTAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGc                >  1:1187529/6‑150 (MQ=255)
                                                                                                                                        gctgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCa               >  1:1521216/5‑150 (MQ=255)
                                                                                                                                        gctgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCa               >  1:1128885/5‑150 (MQ=255)
                                                                                                                                         ctgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCAt              >  1:425303/4‑150 (MQ=255)
                                                                                                                                           aCATAATCCTGGGGATTATGCAGCGGAAAACCTCCAACGTAGGGTTTAATCAGTTTGATTAACGGCTGATAGTAGATAGGTATCAGCGGTGCCTGTTTGTTGATGATCACTTCTGCCGGCTGATTCAACGCATTACGCTTTGTCGCATCa            >  2:210005/1‑150 (MQ=255)
                                                                                                                                           aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCCCTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCa            >  1:1355284/1‑150 (MQ=255)
                                                                                                                                           aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCa            >  1:202770/1‑150 (MQ=255)
                                                                                                                                              tAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTg         >  2:57532/1‑150 (MQ=255)
                                                                                                                                              tAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTg         >  1:1157342/1‑150 (MQ=255)
                                                                                                                                                aTCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTTCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGAt       >  1:113191/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCTGCG  >  NC_000913/3166006‑3166304

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: