Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 179,327 A→T 27.6% I31F (ATC→TTC)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913179,3270AT27.6% 287.3 / 140.7 188I31F (ATC→TTC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (69/67);  new base T (24/28);  total (93/95)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 188 total aligned reads displayed.

GAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTG  >  NC_000913/179231‑179425
                                                                                                |                                                                                                  
gagGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCt                                                                                                 <  2:1695865/100‑1 (MQ=255)
  ggACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTc                                                                                               >  2:1461770/1‑100 (MQ=255)
     cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                            >  2:176279/1‑100 (MQ=255)
      aTGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCgaga                                                                                           <  2:2827584/100‑1 (MQ=255)
      aTGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCg                                                                                              <  2:142548/97‑1 (MQ=255)
       tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                          <  2:877423/100‑1 (MQ=255)
       tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                          >  1:391750/1‑100 (MQ=255)
       tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGagag                                                                                          <  1:2193605/100‑1 (MQ=255)
         gCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGagag                                                                                          <  2:1324109/98‑1 (MQ=255)
          cacaGATTGATTTGCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGa                                                                                       >  1:2790766/1‑100 (MQ=255)
          cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGa                                                                                       <  2:518892/100‑1 (MQ=255)
            caGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATc                                                                                     <  1:3746511/100‑1 (MQ=255)
                ttgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCgg                                                                                 <  2:3325991/100‑1 (MQ=255)
                ttgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCgg                                                                                 <  2:1305342/100‑1 (MQ=255)
                ttgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCgg                                                                                 <  2:2022632/100‑1 (MQ=255)
                 tgattTCCGAACAGAAATAAACTGGCGTCGTCGTTACCGTTCACCGCCGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGCGCGCgg                                                                                 >  2:3118527/1‑99 (MQ=255)
                 tgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCggg                                                                                <  1:3486389/100‑1 (MQ=255)
                  gattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGc                                                                               <  1:3832901/100‑1 (MQ=255)
                  gattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGc                                                                               <  1:1562478/100‑1 (MQ=255)
                   attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                              <  2:1515919/100‑1 (MQ=255)
                   attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                              <  1:2995272/100‑1 (MQ=255)
                   attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCg                                                                              <  2:1071966/100‑1 (MQ=255)
                     ttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTa                                                                            <  2:2485274/100‑1 (MQ=255)
                     ttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTa                                                                            <  1:443366/100‑1 (MQ=255)
                      tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTat                                                                           >  2:2213042/1‑100 (MQ=255)
                      tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTat                                                                           >  2:3160836/1‑100 (MQ=255)
                       ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatc                                                                          >  1:2265681/1‑100 (MQ=255)
                        cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatca                                                                         >  2:859898/1‑100 (MQ=255)
                        cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACAGCAGGGAGATAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGAACGCGGGCGt                                                                             >  2:3931554/1‑96 (MQ=255)
                         gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                        >  2:662624/1‑100 (MQ=255)
                         gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatcat                                                                        <  2:1976986/100‑1 (MQ=255)
                         gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatcat                                                                        >  1:114809/1‑100 (MQ=255)
                         gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatcat                                                                        >  1:3741087/1‑100 (MQ=255)
                          aaaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatcatc                                                                       >  1:850652/1‑100 (MQ=255)
                           tctaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                      <  2:1728860/97‑1 (MQ=255)
                             aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                    <  2:4119824/100‑1 (MQ=255)
                              aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCaaga                                                                   >  1:1728872/1‑98 (MQ=255)
                                aaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAActct                                                                 >  2:2142196/1‑100 (MQ=255)
                                    aaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                             >  2:3507393/1‑100 (MQ=255)
                                     aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                            <  2:619510/100‑1 (MQ=255)
                                     aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                            <  1:2140266/100‑1 (MQ=255)
                                         gCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAtt                                                        <  1:253339/100‑1 (MQ=255)
                                         gCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAtt                                                        >  2:2517400/1‑100 (MQ=255)
                                          cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc                                                       >  1:256360/1‑100 (MQ=255)
                                          cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc                                                       >  1:688234/1‑100 (MQ=255)
                                          cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGCCGTATCCTCAACTCTGCGGCAATtc                                                       >  1:3986846/1‑100 (MQ=255)
                                          cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGAGAGTATCATCAACTCTCCGGCAATtc                                                       >  2:2040412/1‑100 (MQ=255)
                                              gtcgtTACCGTTCACCGCAGGGCGTTAAGACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAGTCTCCGGCAATtcgtcg                                                   <  2:2872280/100‑1 (MQ=255)
                                              gtcgtTACCGTTCACCGCAGGGCGTTAAAGCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                   <  2:838531/100‑1 (MQ=255)
                                               tcgaTACCGTTCCCCGCAGGGCGTTAAAACCGCACATGAGATCCTTCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                   <  2:3630568/99‑1 (MQ=255)
                                                cgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCgtcgtc                                                 >  1:3431029/1‑100 (MQ=255)
                                                   tACCGTTCACCGCAGGGCGTTAAACCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                               <  2:3962638/99‑1 (MQ=255)
                                                   tACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc                                              <  1:132002/100‑1 (MQ=255)
                                                   tACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc                                              <  1:97834/100‑1 (MQ=255)
                                                    aCCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGca                                             >  2:256183/1‑100 (MQ=255)
                                                     ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                            >  1:3243137/1‑100 (MQ=255)
                                                       gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                          >  1:621622/1‑100 (MQ=255)
                                                       gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                          >  2:2951642/1‑100 (MQ=255)
                                                       gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                          >  1:87362/1‑100 (MQ=255)
                                                       gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaac                                           >  2:356547/1‑99 (MQ=255)
                                                          cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                       >  1:4022741/1‑100 (MQ=255)
                                                          cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                       <  1:62832/100‑1 (MQ=255)
                                                           aCCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAg                                      <  1:1191799/100‑1 (MQ=255)
                                                            ccGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAACAGa                                     >  2:1630104/1‑100 (MQ=255)
                                                             cGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                    <  1:858681/100‑1 (MQ=255)
                                                              gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc                                   <  2:4094963/100‑1 (MQ=255)
                                                                aGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCa                                 <  1:176280/100‑1 (MQ=255)
                                                                 gggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATAGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAg                                >  2:456414/1‑100 (MQ=255)
                                                                 gggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAg                                >  1:3585191/1‑100 (MQ=255)
                                                                  ggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                               <  2:1805056/100‑1 (MQ=255)
                                                                  ggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                               <  2:3993876/100‑1 (MQ=255)
                                                                    agTTACACCCGAACATAAGATCATGCGGATCTTCGAGAGCGATCGAGGGCGTATCATCAACTCTCCGGCATTGCGTCGTCTGCAACAGAAGACCCAGGtt                             >  2:2403307/2‑100 (MQ=255)
                                                                        aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                         <  2:1344081/100‑1 (MQ=255)
                                                                        aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                         >  2:2099875/1‑100 (MQ=255)
                                                                        aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                         <  1:2912506/100‑1 (MQ=255)
                                                                         ctaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                        <  1:207917/98‑1 (MQ=255)
                                                                         aaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                        <  2:1570852/100‑1 (MQ=255)
                                                                         aaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                         <  1:3390017/99‑1 (MQ=255)
                                                                          acccGAACATGAGCTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                        <  2:3384326/97‑1 (MQ=255)
                                                                            ccGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACt                     >  1:4110204/1‑100 (MQ=255)
                                                                             cGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                    <  2:3730216/100‑1 (MQ=255)
                                                                              gAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                   >  1:1433747/1‑100 (MQ=255)
                                                                              gAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                   <  1:766451/100‑1 (MQ=255)
                                                                              gAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                   <  2:3982966/100‑1 (MQ=255)
                                                                                 cATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgc                >  2:2855017/1‑100 (MQ=255)
                                                                                 cATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCACCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg                 <  1:1845079/99‑1 (MQ=255)
                                                                                 cATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgc                >  2:3249121/1‑100 (MQ=255)
                                                                                  aTGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcg               >  1:2632738/1‑100 (MQ=255)
                                                                                   tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc              <  2:3039952/100‑1 (MQ=255)
                                                                                    gagcTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCACCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa             <  2:3880591/100‑1 (MQ=255)
                                                                                      gaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt           >  2:3588093/1‑100 (MQ=255)
                                                                                      gaTCCTGCGCTTCTTCGAGAGCGATAGCGCGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGAAACAAAAGACCCAGCTTTTTCCACTGGAGCGCAAt           <  2:2528453/100‑1 (MQ=255)
                                                                                       aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATg          <  1:1803443/100‑1 (MQ=255)
                                                                                        tCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgc         <  2:2992226/100‑1 (MQ=255)
                                                                                          cgcaggATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccg       >  2:1436211/5‑100 (MQ=255)
                                                                                          cTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccg       >  1:3472250/1‑100 (MQ=255)
                                                                                             cGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc     >  2:3794370/1‑99 (MQ=255)
                                                                                               gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTg  <  1:1155426/100‑1 (MQ=255)
                                                                                               gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTg  >  1:1639813/1‑100 (MQ=255)
                                                                                               gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGATTTTCCACTGGAGCGCAATGCCGCCGTg  >  2:1905478/1‑100 (MQ=255)
                                                                                                |                                                                                                  
GAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTG  >  NC_000913/179231‑179425

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