Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,820,566 T→A 5.6% I217F (ATC→TTC)  chbC ← N,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,820,5660TA 5.6% 219.4 / 3.6 71I217F (ATC→TTC) chbCN,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base A (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

ACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCG  >  NC_000913/1820419‑1820713
                                                                                                                                                   |                                                                                                                                                   
aCGGCGTCATATTGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATc                                                                                                                                                    <  6:83700/149‑1 (MQ=255)
aCGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATc                                                                                                                                                    <  7:367304/149‑1 (MQ=255)
   gCGTCATAATGCCGTTGTCCAGTGCGGGCCGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGGGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGa                                                                                                                                                 <  6:282932/149‑1 (MQ=255)
         tAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGCGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCCCGCTAACCAACGATGCCAGTGGGGATCAGATGCTATCCATAATGATCTGATGGAAg                                                                                                                                           >  8:426476/1‑149 (MQ=255)
               cGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTg                                                                                                                                     >  7:95326/1‑149 (MQ=255)
                     ccAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAg                                                                                                                               >  2:304/1‑149 (MQ=255)
                      cAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAgg                                                                                                                              >  8:134169/1‑149 (MQ=255)
                      cAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAg                                                                                                                               <  2:235485/148‑1 (MQ=255)
                           cGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTAtt                                                                                                                          <  8:95326/148‑1 (MQ=255)
                               cagcgccagcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACg                                                                                                                     <  7:67745/149‑1 (MQ=255)
                                    ccagcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAg                                                                                                                >  6:363805/1‑149 (MQ=255)
                                    ccagcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAg                                                                                                                >  8:140260/1‑149 (MQ=255)
                                      agcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATgg                                                                                                                                              >  4:127789/1‑117 (MQ=255)
                                      agcgcGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATgg                                                                                                                                              <  3:127789/117‑1 (MQ=255)
                                           gcCATGTATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCCCAACGCTACCCAACGAAGACAGTGGGGTTGAGATGGTATCCATAATGATCTGAGGGATGTTGGTGCCCCAGGGATTCACCGCCCAGGCaatca                                                                                                         >  8:361395/1‑147 (MQ=255)
                                            ccATCAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCaataat                                                                                                        <  3:341592/149‑1 (MQ=255)
                                               tGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATccc                                                                                                     >  8:125814/1‑149 (MQ=255)
                                                       gaagaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                              >  2:265965/1‑148 (MQ=255)
                                                       gaagaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                              >  5:399373/1‑148 (MQ=255)
                                                          gaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATGATCCCCATCACGGaaa                                                                                          <  1:168440/149‑1 (MQ=255)
                                                          gaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGaaa                                                                                          <  2:386437/149‑1 (MQ=255)
                                                          gaaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGaaa                                                                                          <  6:399373/149‑1 (MQ=255)
                                                             ccAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCAGAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCCCGGAAAGaa                                                                                       >  5:261635/1‑149 (MQ=255)
                                                              cAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGaat                                                                                      <  7:134169/149‑1 (MQ=255)
                                                               agagCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGaata                                                                                     <  2:47770/149‑1 (MQ=255)
                                                                gagCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCCACGCCCAGGCAATAATCCCCATCACGGAAAGaataa                                                                                    >  2:139815/1‑149 (MQ=255)
                                                                 agCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGAACTGATGGaa                                                                                                                                            >  7:391839/1‑92 (MQ=255)
                                                                 agCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGAACTGATGGaa                                                                                                                                            >  7:386163/1‑92 (MQ=255)
                                                                 agCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGAACTGATGGaa                                                                                                                                            <  8:386163/92‑1 (MQ=255)
                                                                 agCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGAACTGATGGaa                                                                                                                                            <  8:391839/92‑1 (MQ=255)
                                                                        gAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCCGTGGGGGTGAGAGGGTATACATAATGATCTGATGGAAGTGGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCg                                                                            <  6:223379/149‑1 (MQ=255)
                                                                         aaCAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCgg                                                                           <  5:363805/149‑1 (MQ=255)
                                                                           cAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTg                                                                                                                                                  >  2:407491/1‑76 (MQ=255)
                                                                           cAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTg                                                                                                                                                  <  1:407491/76‑1 (MQ=255)
                                                                           cAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCggg                                                                          >  6:312875/1‑148 (MQ=255)
                                                                                  cacaTAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATACTCCCCATCACGGAAAGAATAATAACGCCGGGAATTAATg                                                                  >  7:332056/1‑149 (MQ=255)
                                                                                  cacaTAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATg                                                                  >  7:99495/1‑149 (MQ=255)
                                                                                       aGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGa                                                                         >  2:175483/1‑137 (MQ=255)
                                                                                       aGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGa                                                                         <  1:175483/137‑1 (MQ=255)
                                                                                          cccagcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGaa                                                           <  1:265965/148‑1 (MQ=255)
                                                                                          cccagcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAg                                                          <  8:99495/149‑1 (MQ=255)
                                                                                          cccagcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGACGGGAATTAATGCCGAGAAg                                                          >  3:151202/1‑149 (MQ=255)
                                                                                           ccagcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATAATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAg                                                          >  2:416819/1‑148 (MQ=255)
                                                                                           ccagcccaCCACGCTACCCAAAGATGCCAGAGGGGTTGAGATGGAATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGCGAATTAATGCCGAGAAg                                                          >  3:218414/1‑148 (MQ=255)
                                                                                              gcccaCCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAAc                                                      >  1:391112/1‑149 (MQ=255)
                                                                                                ccaccaCGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGCGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAAcgc                                                    <  6:119686/149‑1 (MQ=255)
                                                                                                 caccacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAAcgc                                                    >  7:383995/1‑148 (MQ=255)
                                                                                                    cacGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATa                                                <  7:125814/149‑1 (MQ=255)
                                                                                                        cTAACCAACGATGCCACTGGGGTTGAGATGGTATCCAAAAAGATCTGATGGAAGGTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAGGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGa                                            <  5:315339/149‑1 (MQ=255)
                                                                                                              aaCGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCGCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAgg                                      >  4:13366/1‑149 (MQ=255)
                                                                                                              aaCGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGaa                                           <  4:385852/144‑1 (MQ=255)
                                                                                                              aaCGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGaa                                           >  3:385852/1‑144 (MQ=255)
                                                                                                              aaCGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAgg                                      >  3:287888/1‑149 (MQ=255)
                                                                                                                cGATGCCAGTGGGGTTGAGATGGTATCCATAGTGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTa                                    >  6:301013/1‑149 (MQ=255)
                                                                                                                cGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAg                                          >  1:52834/1‑143 (MQ=255)
                                                                                                                cGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAg                                          <  2:52834/143‑1 (MQ=255)
                                                                                                                  aTGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACg                                  <  7:140260/149‑1 (MQ=255)
                                                                                                                       aGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTc                             <  8:383995/149‑1 (MQ=255)
                                                                                                                                  gaTGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTaa                  >  4:373642/1‑149 (MQ=255)
                                                                                                                                      gTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGa               <  3:373642/148‑1 (MQ=255)
                                                                                                                                      gTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGAc              >  8:110162/1‑149 (MQ=255)
                                                                                                                                      gTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGAc              <  2:130905/149‑1 (MQ=255)
                                                                                                                                           cATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAAt         >  2:255626/1‑149 (MQ=255)
                                                                                                                                             tAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTg       >  2:116035/1‑149 (MQ=255)
                                                                                                                                                tGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTgcgg    >  6:380695/1‑149 (MQ=255)
                                                                                                                                                 gATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCg                                                               <  8:180568/89‑1 (MQ=255)
                                                                                                                                                 gATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCg                                                               <  8:186353/89‑1 (MQ=255)
                                                                                                                                                 gATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCg                                                               >  7:180568/1‑89 (MQ=255)
                                                                                                                                                 gATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCg                                                               >  7:186353/1‑89 (MQ=255)
                                                                                                                                                 gATCTGATGGAAGTTGGTGCCACAGGTATTCAACGCCCAGGCAATAATCCCCATCACTGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTgcggc   >  5:85245/1‑149 (MQ=255)
                                                                                                                                                  aTCTGATGGAAGTTGGTGCCCCAGGTAATCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTgcggc   <  1:83996/148‑1 (MQ=255)
                                                                                                                                                   tCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTgcggcg  >  1:396176/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCG  >  NC_000913/1820419‑1820713

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: