Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 179,373 G→A 5.9% R46H (CGT→CAT)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913179,3730GA 5.9% 319.8 / 4.8 101R46H (CGT→CAT) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (43/51);  new base A (3/3);  total (46/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CGCGGGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAG  >  NC_000913/179224‑179522
                                                                                                                                                     |                                                                                                                                                     
cgcgGGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCATGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg                                                                                                                                                       <  2:1224042/150‑1 (MQ=255)
  cgGGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc                                                                                                                                                     <  2:332378/150‑1 (MQ=255)
    ggAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt                                                                                                                                                   >  2:1553640/1‑150 (MQ=255)
     gagaGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCgtcgtc                                                                                                                                                  <  2:226970/150‑1 (MQ=255)
     gagaGGACATGGCACAGATTGATTTCAGAAAAAAAATAAACTGGCATCGGCGTTACCGGTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCgtcgtc                                                                                                                                                  <  2:74674/150‑1 (MQ=255)
       gagGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                                                                                                                                <  1:1182808/150‑1 (MQ=255)
       cagGACATTGCACAGATTGATTTACGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                                                                                                                                <  2:1693368/149‑1 (MQ=255)
           aCATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGAGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaac                                                                                                                                            >  2:1198592/1‑150 (MQ=255)
            cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                                                                                                                           <  2:731537/150‑1 (MQ=255)
            cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                                                                                                                           >  1:1680894/1‑150 (MQ=255)
            cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                                                                                                                           <  1:1317748/150‑1 (MQ=255)
              tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGTCAATTCGTCGTATGCAACaaa                                                                                                                                         >  2:474919/1‑150 (MQ=255)
              tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaa                                                                                                                                         >  1:559206/1‑150 (MQ=255)
              tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaa                                                                                                                                         >  1:650662/1‑150 (MQ=255)
                gCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAg                                                                                                                                       <  1:1158991/150‑1 (MQ=255)
                 cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCATCGTCTGCAACAAAAgc                                                                                                                                      >  1:1031678/1‑149 (MQ=255)
                 cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGAGGATCTTCGAGAGCGATCGCGGGCGTATCATCACCTCTCCGGCAATTCATCGTCTGCAACAAAAGa                                                                                                                                      >  1:1128834/1‑150 (MQ=255)
                 cacaGATTGATTTCCGAAAAAAAATAAACTGCCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGACATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATAATCAACTCTCCGGCAATTCATCGTCTGCAACTAAAGa                                                                                                                                      >  1:1128829/1‑150 (MQ=255)
                  acaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                                                                                                                     <  1:1195505/150‑1 (MQ=255)
                   caGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc                                                                                                                                    >  2:492195/1‑150 (MQ=255)
                   caGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc                                                                                                                                    >  2:336251/1‑150 (MQ=255)
                     gattgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCa                                                                                                                                  >  1:1499638/1‑150 (MQ=255)
                          attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                                                                                                                             <  2:650662/150‑1 (MQ=255)
                          attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                                                                                                                             <  2:559206/150‑1 (MQ=255)
                          attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                                                                                                                             <  1:434578/150‑1 (MQ=255)
                            ttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttttt                                                                                                                           >  1:896028/1‑150 (MQ=255)
                             tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCCGATCTACGAGAGCGATCGCGGGCGTATCATCAACTCTCAGGCAATTCGACGTCGGCAACAAAAGACCCAGGAGTTTc                                                                                                                          >  1:595552/1‑150 (MQ=255)
                              ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                                                                                                                         <  1:1287879/150‑1 (MQ=255)
                              ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                                                                                                                         <  2:1598584/150‑1 (MQ=255)
                              ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                                                                                                                         >  2:407776/1‑150 (MQ=255)
                               cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCa                                                                                                                        >  2:235418/1‑150 (MQ=255)
                                   aaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGGTAAAAGCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTAGCATCAACTCTCCGGCAATTCGTCGGCTGCAACAAAAGACCCAGGTTTTTCCACTgg                                                                                                                    <  1:282616/150‑1 (MQ=255)
                                     aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAATCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg                                                                                                                  <  1:73130/150‑1 (MQ=255)
                                     aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg                                                                                                                  >  1:681975/1‑150 (MQ=255)
                                     aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg                                                                                                                  <  2:125560/150‑1 (MQ=255)
                                        aTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc                                                                                                               >  2:409888/1‑150 (MQ=255)
                                                 cATCTTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc                                                                                                      <  2:1222475/150‑1 (MQ=255)
                                                 cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCACCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGCGCGCAAAgccgcc                                                                                                      >  1:1251048/1‑150 (MQ=255)
                                                 cATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc                                                                                                      >  1:1378707/1‑150 (MQ=255)
                                                  aTCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGccgccg                                                                                                     <  1:336251/150‑1 (MQ=255)
                                                    cgtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCGACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTg                                                                                                   <  1:287526/150‑1 (MQ=255)
                                                     gtcgtTACCGTTCACCGCAGGGCGTTAAATCCGCACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCGTCAGCGCTCCGGCAATTCGTCGTGGGCACCGAAAGGCCCAGGTTTTTCCCCGGGCGCGCAATGccgccgggg                                                                                                  >  2:1483952/1‑147 (MQ=255)
                                                     gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGCCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgc                                                                                                  >  2:1057367/1‑150 (MQ=255)
                                                       cgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgcgc                                                                                                >  2:1712906/1‑150 (MQ=255)
                                                        gtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCTGCCGTGCGCa                                                                                               <  1:490076/150‑1 (MQ=255)
                                                        gtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCa                                                                                               <  1:559680/150‑1 (MQ=255)
                                                            ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCAcgcg                                                                                           >  2:591659/1‑150 (MQ=255)
                                                              gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTc                                                                                         <  1:409888/150‑1 (MQ=255)
                                                               ttCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCt                                                                                        <  2:450849/150‑1 (MQ=255)
                                                                tCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCtt                                                                                       <  1:918331/150‑1 (MQ=255)
                                                                 cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCATGCGCACGCGTCTTa                                                                                      >  2:933440/1‑150 (MQ=255)
                                                                      cAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCAc                                                                                 <  1:223539/150‑1 (MQ=255)
                                                                            gtaaaaaCAGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGGCTGCAACATAAGACCCAGGTTTTTCCACTGGAGAGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATg                                                                           <  2:1041989/147‑1 (MQ=255)
                                                                            gTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATg                                                                           <  2:820148/150‑1 (MQ=255)
                                                                            gTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGGTTTTCCACTGGAGCGCACTGCCGCCGTGCGCACGCGTCTTACCCACTCGATg                                                                           <  2:95564/150‑1 (MQ=255)
                                                                              tAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGa                                                                         >  2:212537/1‑150 (MQ=255)
                                                                               aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGaa                                                                        >  1:951773/1‑150 (MQ=255)
                                                                                aaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCCACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGACGTGCGCACGCGTCTTACCCACTCGATGGAAg                                                                       <  2:811998/150‑1 (MQ=255)
                                                                                aaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAg                                                                       <  2:1378707/150‑1 (MQ=255)
                                                                                 aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGt                                                                      >  2:772310/1‑150 (MQ=255)
                                                                                      aaCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCcagc                                                                 >  1:1695019/1‑150 (MQ=255)
                                                                                      aaCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCcagc                                                                 <  1:1591399/150‑1 (MQ=255)
                                                                                          tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGt                                                             <  1:407776/150‑1 (MQ=255)
                                                                                          tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGt                                                             >  1:720302/1‑150 (MQ=255)
                                                                                           gagaTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACCCCTCTTACCCACACGATGGAAGTCCAGCAGGTg                                                            >  2:228352/1‑150 (MQ=255)
                                                                                                  tGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGCCCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACGCACTCGATGGAAGTCCAGCAGGTGGGGCGCt                                                     >  2:1140082/1‑150 (MQ=255)
                                                                                                  tGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCt                                                     >  1:947817/1‑150 (MQ=255)
                                                                                                       aTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGCCCCAGGTTTTTCCACTGGAGCGCCATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATc                                                >  1:714327/1‑150 (MQ=255)
                                                                                                       aTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCGTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATc                                                >  1:708059/1‑150 (MQ=255)
                                                                                                         cTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCATGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGc                                              <  1:933440/150‑1 (MQ=255)
                                                                                                          ttCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGcc                                             <  2:947817/150‑1 (MQ=255)
                                                                                                          ttCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGGGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGcc                                             <  2:371325/150‑1 (MQ=255)
                                                                                                           tCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCa                                            <  1:235418/150‑1 (MQ=255)
                                                                                                             gagaGCGATCGCGGGTGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGAGGGAACTCCAGCCGGTGGGGCGCTACATCGCCaaa                                          >  2:848486/1‑150 (MQ=255)
                                                                                                             gagaGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCaaa                                          >  1:1032142/1‑150 (MQ=255)
                                                                                                              cgagCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                         <  2:617761/149‑1 (MQ=255)
                                                                                                              agagCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                         <  1:634296/150‑1 (MQ=255)
                                                                                                                 gCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGaaa                                      >  1:1487565/1‑150 (MQ=255)
                                                                                                                 ccGATCGCGGGCGTATCATCAACTCTCCGGCAAATAGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCGGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGGGGGGCGCTACATCGCCAAAGaaa                                      <  1:474919/149‑1 (MQ=255)
                                                                                                                      cgcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGGCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTa                                 <  1:1553640/150‑1 (MQ=255)
                                                                                                                                 catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaa                      >  2:1277118/1‑150 (MQ=255)
                                                                                                                                  atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaaa                     <  1:980033/150‑1 (MQ=255)
                                                                                                                                  atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaaa                     <  1:591659/150‑1 (MQ=255)
                                                                                                                                  atcaACTCTCCGCCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGACGGTGGGGCGCTACATCGCCAAAGAAATGTTAAGCCGTCTGaaa                     <  2:96758/150‑1 (MQ=255)
                                                                                                                                  atcaAATCTCCGGCAATTCTTCTTCTTCAACAAAAGAGCCAGGTATTTCCACTGGAGCGCAATGCCGACGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaaa                     <  2:88324/150‑1 (MQ=255)
                                                                                                                                   tcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAag                    <  1:1712906/150‑1 (MQ=255)
                                                                                                                                     aaCTCTCCGGCACTTCGTCGGCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                  <  1:848486/150‑1 (MQ=255)
                                                                                                                                     aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                  <  2:1032142/150‑1 (MQ=255)
                                                                                                                                     aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                  >  1:1469188/1‑150 (MQ=255)
                                                                                                                                      aCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGc                 <  1:1198592/150‑1 (MQ=255)
                                                                                                                                           ccGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTCCCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACTTCGCCAAAGAAATTTTAAGCAGTCTGAAAGAGCGTaaa            >  2:1491159/1‑150 (MQ=255)
                                                                                                                                              gCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTa         >  1:963866/1‑150 (MQ=255)
                                                                                                                                              gCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTa         >  1:669862/1‑150 (MQ=255)
                                                                                                                                               cAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTAc        <  2:681975/150‑1 (MQ=255)
                                                                                                                                                aaTTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCCACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACt       <  2:243141/150‑1 (MQ=255)
                                                                                                                                                  ttCATCGTCTGCATCAAAATACCCAGTTTTTTCGACTGGAGCGCAATGCAGCCGTGCGCACGCGTCTTACGCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCATAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTgg     <  2:1128834/150‑1 (MQ=255)
                                                                                                                                                  ttCATCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTgg     <  2:1031678/150‑1 (MQ=255)
                                                                                                                                                  ttCATCGGCTGCAACAAAATACCCAGGGATTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTAGATGTAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTATAAGCCGTCTGAAAGAGCTTAAATTACTgg     <  2:1128829/150‑1 (MQ=255)
                                                                                                                                                    cgtcgtCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGaa   <  2:896028/150‑1 (MQ=255)
                                                                                                                                                     gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAg  >  2:1432653/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGGGAGAGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAG  >  NC_000913/179224‑179522

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: ‑