Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 179,327 A→T 72.9% 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,3270AT72.9% 237.0 / 101.1 137I31F (ATC→TTC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base T (56/43);  total (75/62)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 142 total aligned reads displayed.

GGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGC  >  NC_000913/179228‑179421
                                                                                                   |                                                                                              
   gagGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGa                                                                                                <  1:1789118/97‑1 (MQ=255)
   gagGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCt                                                                                             >  2:2511613/1‑100 (MQ=255)
   gagGACATGGCACAGATAGATTACCTAAAAAAAATAAACTGTCATCGTCGTGACCGTTCACCGCAGGGCGTTAAAGCCGAACATGAGATCCTGCGGTTc                                                                                              <  1:1158338/99‑1 (MQ=255)
     ggACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCAACGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGttcttc                                                                                           >  2:603432/1‑100 (MQ=255)
      gACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCg                                                                                          <  1:1698480/100‑1 (MQ=255)
      gACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCg                                                                                          <  2:1963615/100‑1 (MQ=255)
       aCATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCga                                                                                         >  1:935485/1‑100 (MQ=255)
        cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                        <  2:2494644/100‑1 (MQ=255)
        cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCgag                                                                                        >  1:776210/1‑100 (MQ=255)
         aTGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgaga                                                                                       >  1:473473/1‑100 (MQ=255)
          tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGCACATGAGATCCTGCGGTTCTTCGagag                                                                                      >  2:25945/1‑100 (MQ=255)
          tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                      <  2:1796289/100‑1 (MQ=255)
          tGGCACAGATTGATTTCAGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                      <  2:2665999/100‑1 (MQ=255)
             cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGa                                                                                   >  1:947522/1‑100 (MQ=255)
             cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGa                                                                                   >  2:2670181/1‑100 (MQ=255)
             cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGa                                                                                   >  2:1485381/1‑100 (MQ=255)
               caGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATc                                                                                 >  2:1383968/1‑100 (MQ=255)
                aGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcg                                                                                >  2:3290433/1‑100 (MQ=255)
                aGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcg                                                                                <  2:473475/100‑1 (MQ=255)
                   ttgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCgg                                                                             >  1:448519/1‑100 (MQ=255)
                   ttgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCgg                                                                             >  2:1725395/1‑100 (MQ=255)
                      attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCg                                                                          <  1:1455350/100‑1 (MQ=255)
                      attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCg                                                                          <  2:1689083/100‑1 (MQ=255)
                       tttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcgc                                                                               >  1:2120422/1‑94 (MQ=255)
                       tttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGt                                                                         >  1:3158516/1‑100 (MQ=255)
                         tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTat                                                                       >  1:1664371/1‑100 (MQ=255)
                          ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatc                                                                      <  2:1567578/100‑1 (MQ=255)
                          ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatc                                                                      <  2:2087739/100‑1 (MQ=255)
                          ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatc                                                                      <  1:2406128/100‑1 (MQ=255)
                            gAAAAAAAATAAACTGGCATTGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                    <  1:953093/100‑1 (MQ=255)
                            gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                    >  2:2584525/1‑100 (MQ=255)
                             aaaaaaaaTAAACTGGCATTGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                    <  1:951596/99‑1 (MQ=255)
                              aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                  <  1:2936241/100‑1 (MQ=255)
                              aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                  <  2:110649/100‑1 (MQ=255)
                               aaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCaa                                                                 <  2:2840022/100‑1 (MQ=255)
                               aaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCaa                                                                 <  2:407619/100‑1 (MQ=255)
                               aaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCaa                                                                 <  2:2131397/100‑1 (MQ=255)
                               aaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGACGTTACAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCaa                                                                 >  1:435282/1‑100 (MQ=255)
                                aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                <  2:837424/100‑1 (MQ=255)
                                aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                <  1:688973/100‑1 (MQ=255)
                                aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                <  1:347929/100‑1 (MQ=255)
                                 aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                               >  2:1758147/1‑100 (MQ=255)
                                 aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                               >  2:1162940/1‑100 (MQ=255)
                                  aaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAActc                                                              >  2:1873539/1‑100 (MQ=255)
                                  aaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAActc                                                              >  2:2173767/1‑100 (MQ=255)
                                   aaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAActct                                                             >  1:1193636/1‑100 (MQ=255)
                                      aaaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCg                                                          <  2:1538008/100‑1 (MQ=255)
                                       aaCTGGCATCGTCGTTACCGTTCACGGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                         <  1:824587/100‑1 (MQ=255)
                                        aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                        <  2:2242222/100‑1 (MQ=255)
                                         cTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa                                                       >  2:3124300/1‑100 (MQ=255)
                                         cTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa                                                       >  2:1342774/1‑100 (MQ=255)
                                         cTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGGTCTTCGAGAGCGATCACGGGCGTATCATCAACTCTACGGCa                                                       >  1:2705729/1‑100 (MQ=255)
                                            gCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAtt                                                    >  2:835360/1‑100 (MQ=255)
                                              aTCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg                                                  >  1:2797502/1‑100 (MQ=255)
                                               tcgtcgTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt                                                 >  2:156661/1‑100 (MQ=255)
                                                cgtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc                                                >  2:3157692/1‑100 (MQ=255)
                                                cgtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc                                                <  1:156661/100‑1 (MQ=255)
                                                 gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATttgtc                                                >  2:2896323/1‑99 (MQ=255)
                                                 gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcn                                               <  1:2066335/100‑2 (MQ=255)
                                                 gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                               <  1:235309/100‑1 (MQ=255)
                                                  tcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt                                              >  2:908697/1‑100 (MQ=255)
                                                   cgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCgtcgtc                                             >  1:344664/1‑100 (MQ=255)
                                                      tACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc                                          >  2:2597138/1‑100 (MQ=255)
                                                         cGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaac                                       >  1:2503738/1‑100 (MQ=255)
                                                          gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                      >  2:1300523/1‑100 (MQ=255)
                                                          gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                      >  1:878759/1‑100 (MQ=255)
                                                           ttCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaacaa                                     >  2:503423/1‑100 (MQ=255)
                                                            tCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACgaa                                    <  2:83494/100‑1 (MQ=255)
                                                             cACCGCAGGGCGGTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGAGCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                   >  1:565121/1‑100 (MQ=255)
                                                              cttccgatctgTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAg                                  <  1:3333268/90‑1 (MQ=255)
                                                              aCCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAg                                  >  2:3274552/1‑100 (MQ=255)
                                                                    gggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAg                            >  2:737742/1‑100 (MQ=255)
                                                                     ggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                           <  2:2067408/100‑1 (MQ=255)
                                                                      gCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgn                           <  2:1593284/99‑2 (MQ=255)
                                                                       cGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtt                         >  1:73743/1‑100 (MQ=255)
                                                                        gTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAATAAAAGagatcggaag                        >  2:3333296/1‑91 (MQ=255)
                                                                        gTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                        >  1:201638/1‑100 (MQ=255)
                                                                           aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                     <  2:2775901/100‑1 (MQ=255)
                                                                           aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                     >  2:530663/1‑100 (MQ=255)
                                                                           aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                     >  2:1240576/1‑100 (MQ=255)
                                                                            aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                    >  2:2093023/1‑100 (MQ=255)
                                                                            aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                    <  2:1100253/100‑1 (MQ=255)
                                                                            aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                    <  2:1812674/100‑1 (MQ=255)
                                                                             aaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                   >  2:1618303/1‑100 (MQ=255)
                                                                              aCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                  >  1:1230974/1‑100 (MQ=255)
                                                                              aCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                  <  1:2314436/100‑1 (MQ=255)
                                                                              aCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                  <  1:115353/100‑1 (MQ=255)
                                                                              aCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                  >  2:2925976/1‑100 (MQ=255)
                                                                                 gAACATGAGATCCTGCGGTTCTTCGAGAGCGATCTCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGATTTTCCACTgg               <  1:1782436/100‑1 (MQ=255)
                                                                                   aCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg             <  1:2043979/100‑1 (MQ=255)
                                                                                   aCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg             >  1:3367387/1‑100 (MQ=255)
                                                                                    cATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgc            >  1:210192/1‑100 (MQ=255)
                                                                                    cATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg             >  1:2665416/1‑99 (MQ=255)
                                                                                      tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc          >  2:70319/1‑100 (MQ=255)
                                                                                      tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACCCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc          >  1:2463355/1‑100 (MQ=255)
                                                                                      tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc          >  1:2854434/1‑100 (MQ=255)
                                                                                         gaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt       >  2:2861382/1‑100 (MQ=255)
                                                                                           tCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgc     >  1:2604828/1‑100 (MQ=255)
                                                                                           tCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgc     >  1:2955448/1‑100 (MQ=255)
                                                                                              tGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTTGAGCGCAATgccgc  <  1:3165460/100‑1 (MQ=255)
                                                                                                   |                                                                                              
GGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGC  >  NC_000913/179228‑179421

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