Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,935,650 C→T 6.5% R222H (CGC→CAC)  zwf ← glucose‑6‑phosphate 1‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,935,6500CT 6.5% 191.1 / 3.2 62R222H (CGC→CAC) zwfglucose‑6‑phosphate 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTGGTCGATA  >  NC_000913/1935504‑1935793
                                                                                                                                                  |                                                                                                                                               
gagaCTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAgcgc                                                                                                                                               >  4:208093/1‑149 (MQ=255)
gagaCTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAgcgc                                                                                                                                               <  2:141098/149‑1 (MQ=255)
     ttCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTc                                                                                                                                          >  7:232494/1‑149 (MQ=255)
     ttCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTc                                                                                                                                          >  2:173164/1‑149 (MQ=255)
       cAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTc                                                                                                                                          <  3:281631/147‑1 (MQ=255)
        aGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGAt                                                                                                                                       >  8:363033/1‑149 (MQ=255)
            cTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATcccc                                                                                                                                   <  2:167194/149‑1 (MQ=255)
            cTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATccc                                                                                                                                    <  2:129855/148‑1 (MQ=255)
            cTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATccc                                                                                                                                    <  7:251255/148‑1 (MQ=255)
             tttCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGTCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGCCCGTCTTTATCAAAATAGCCCCAGCGCCCTTCGATCcaca                                                                                                                                  >  1:91595/1‑149 (MQ=255)
              ttCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCacac                                                                                                                                 >  7:147281/1‑149 (MQ=255)
                      ttCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTtcttct                                                                                                                          <  8:13463/148‑1 (MQ=255)
                      ttCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTtcttct                                                                                                                          <  1:173164/148‑1 (MQ=255)
                       tCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTg                                                                                                                         <  3:102626/148‑1 (MQ=255)
                              ggATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGt                                                                                                                  >  6:160681/1‑148 (MQ=255)
                                   cTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAAtctc                                                                                                            <  1:264708/149‑1 (MQ=255)
                                   cTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAAtctc                                                                                                            <  5:160681/149‑1 (MQ=255)
                                    tGTATGCGCTCAGGTCAGACGGCGGATACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAAtctc                                                                                                            >  1:211040/1‑148 (MQ=255)
                                      tCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAAc                                                                                                         <  6:33076/149‑1 (MQ=255)
                                               aGGTCAGACGACGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGGGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAACTAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAAt                                                                                                <  6:162181/149‑1 (MQ=255)
                                                gttCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTTCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCATAATAGCCCCATCGCCCTTCGATCCCCACTTCTTCTGCCACGGTACTCTCAACATGATCAATg                                                                                               <  1:296863/147‑1 (MQ=255)
                                                ggTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGTGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATg                                                                                               >  8:166575/1‑149 (MQ=255)
                                                 gTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGTGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATgg                                                                                              >  8:238518/1‑149 (MQ=255)
                                                          cggAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGt                                                                                     >  3:133736/1‑149 (MQ=255)
                                                              gaCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCa                                                                                 <  3:208093/149‑1 (MQ=255)
                                                               aCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAg                                                                                >  1:330135/1‑149 (MQ=255)
                                                                cATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAg                                                                                <  5:132750/148‑1 (MQ=255)
                                                                 aTCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATc                                                                                                                                      >  8:286103/1‑93 (MQ=255)
                                                                 aTCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATc                                                                                                                                      <  7:286103/93‑1 (MQ=255)
                                                                 aTCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGt                                                                               >  5:297928/1‑148 (MQ=255)
                                                                  tCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTa                                                                             >  2:361536/1‑149 (MQ=255)
                                                                  tCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATCATCAATGGTGCGATTTTCCCAGtga                                                                             >  7:369328/1‑147 (MQ=255)
                                                                    gCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGTGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTAtt                                                                           <  7:166575/149‑1 (MQ=255)
                                                                               aaGAATTTTCAGGATGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCATAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAggg                                                                <  8:369328/149‑1 (MQ=255)
                                                                               aaGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAgg                                                                 >  2:165849/1‑148 (MQ=255)
                                                                                         agcagGTGGGTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGGCCCAGTTATTCACAAACAGGGAGTTAGCaaa                                                      <  2:91595/149‑1 (MQ=255)
                                                                                          gcagGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGCTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCaaaa                                                     <  7:363033/149‑1 (MQ=255)
                                                                                                      ggATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTca                                                                         >  5:64680/1‑117 (MQ=255)
                                                                                                      ggATAATGTCGCGCATCTGACCGGCTTTATCATAATAGACCCAGCGCCCTTAGATCCCCACTTCTTCTGCAACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTca                                                                         <  6:64680/117‑1 (MQ=255)
                                                                                                           aTGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGtt                                    <  8:232494/149‑1 (MQ=255)
                                                                                                           aTGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGtt                                    <  2:330135/149‑1 (MQ=255)
                                                                                                            tGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGtt                                    >  6:379180/1‑148 (MQ=255)
                                                                                                                gcgcATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCa                               >  1:12509/1‑149 (MQ=255)
                                                                                                                 cgcATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACg                                                   <  8:22531/128‑1 (MQ=255)
                                                                                                                 cgcATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACg                                                   >  7:22531/1‑128 (MQ=255)
                                                                                                                   cATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCg                            <  2:12509/149‑1 (MQ=255)
                                                                                                                       tGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTc                        >  1:101798/1‑149 (MQ=255)
                                                                                                                         aCCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGGAGCGCCAACAGGTTCAGCACCGttt                         <  1:165849/146‑1 (MQ=255)
                                                                                                                         aCCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCtt                      >  5:359151/1‑149 (MQ=255)
                                                                                                                         aCCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCtt                      >  8:327907/1‑149 (MQ=255)
                                                                                                                          ccGGCTTTATCAAAATAGCCCCCTCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCtt                      <  8:147281/148‑1 (MQ=255)
                                                                                                                           cGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGTTTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTa                    >  2:181371/1‑149 (MQ=255)
                                                                                                                            ggCTTTATCAAAATAGCCCCAGTGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTAc                   <  7:238518/149‑1 (MQ=255)
                                                                                                                            ggCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTAc                   <  6:297928/149‑1 (MQ=255)
                                                                                                                                  aTCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGa              >  3:240701/1‑148 (MQ=255)
                                                                                                                                   tCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCTATGGTGCGATTGTCCCAGATATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATa            <  2:211040/149‑1 (MQ=255)
                                                                                                                                   tCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTATACCAAGa              >  3:2646/1‑147 (MQ=255)
                                                                                                                                      aaaTAGCCCCAGAGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTg         >  4:225476/1‑149 (MQ=255)
                                                                                                                                      aaaTACCCCCAGCCCCCTTCGAACCCCACTACTGCTGCCCCGGTAAACTCAACATCATCAAAGGTCCGATAGTCCCATTTAATCACAAACAGGGAGTTAGCAAAACGAAGCGCCAACAGGTTCAGCACACTTTCTTTACCAAGATAgag         >  2:129000/1‑147 (MQ=255)
                                                                                                                                            ccccAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGTTTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTGGTCGAt   <  1:181371/149‑1 (MQ=255)
                                                                                                                                            ccccAGCGCCCTTCGATACCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTGCAGCACCGTTTCTTTACCAAGAGAGGGGTCGa    >  4:327583/1‑148 (MQ=255)
                                                                                                                                             ccccGCGCCCTTCGAGCCCCCCTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGtt                          <  6:373130/125‑1 (MQ=255)
                                                                                                                                             cccAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTGGTCGATa  <  1:129000/149‑1 (MQ=255)
                                                                                                                                             cccAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTGGTCGATa  >  6:123371/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAAGATAGTGGTCGATA  >  NC_000913/1935504‑1935793

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