Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,415 A→C 8.2% L272R (CTG→CGG) ‡ gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,4150AC 8.2% 315.5 / 19.6 97L272R (CTG→CGG) ‡gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base C (4/4);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTG  >  NC_000913/3963268‑3963559
                                                                                                                                                    |                                                                                                                                                
aTAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                 >  2:580939/1‑149 (MQ=255)
aTAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                                                                                                                                 >  3:511448/1‑149 (MQ=255)
 tAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAgg                                                                                                                                                >  2:213134/1‑149 (MQ=255)
    aGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATc                                                                                                                                             <  1:873745/149‑1 (MQ=255)
    aGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATc                                                                                                                                             >  3:555442/1‑149 (MQ=255)
       ggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGcc                                                                                                                                          <  4:541322/149‑1 (MQ=255)
        gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                                                                                                                                         <  1:567537/149‑1 (MQ=255)
         cgTTGAATCTTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCATCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCTGTAAAAATGGCGATCCTGATCGCCat                                                                                                                                        <  3:775211/149‑2 (MQ=255)
          gTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCatt                                                                                                                                       <  1:46947/149‑3 (MQ=255)
           ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                                                                                                                                       <  1:851535/148‑3 (MQ=255)
           ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                                                                                                                                       >  2:851535/1‑146 (MQ=255)
           ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                                                                                                                                      >  4:499824/1‑149 (MQ=255)
           ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                                                                                                                                      <  4:511448/149‑1 (MQ=255)
                  ttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                               >  4:47480/1‑149 (MQ=255)
                  ttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg                                                                                                                               >  3:273405/1‑149 (MQ=255)
                   tACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCggg                                                                                                                              <  1:580939/149‑1 (MQ=255)
                    tcGCAGCCGTACGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGa                                                                                                                             >  3:263205/2‑149 (MQ=255)
                     cGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAATTGGCGATCAGGATCGCCAGACCACTCGGGaa                                                                                                                            <  2:4300/149‑1 (MQ=255)
                     cGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa                                                                                                                            >  1:887152/1‑149 (MQ=255)
                       cAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAaca                                                                                                                          <  1:420575/149‑1 (MQ=255)
                                ggCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCa                                                                                                                  >  4:633532/1‑149 (MQ=255)
                                   tGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGt                                                                                                               >  3:362969/1‑149 (MQ=255)
                                     cGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTc                                                                                                             <  3:474944/149‑1 (MQ=255)
                                      gAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTcc                                                                                                            >  2:376201/1‑149 (MQ=255)
                                       aaTATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCTCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCa                                                                                                           <  3:762932/149‑1 (MQ=255)
                                       aaTATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCa                                                                                                           >  2:93678/1‑149 (MQ=255)
                                            ccTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTc                                                                                                      <  1:542627/149‑1 (MQ=255)
                                            ccTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTc                                                                                                      <  1:619276/149‑1 (MQ=255)
                                                 tCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGccc                                                                                                 >  3:879550/1‑149 (MQ=255)
                                                   gACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGt                                                                                               <  3:848648/149‑1 (MQ=255)
                                                       gCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAt                                                                                           <  3:26333/149‑1 (MQ=255)
                                                        ccAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                          <  1:143295/149‑1 (MQ=255)
                                                         cAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtct                                                                                         >  1:895234/1‑149 (MQ=255)
                                                          aGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtctc                                                                                        <  4:113694/149‑1 (MQ=255)
                                                          aGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtctc                                                                                        <  3:664669/149‑1 (MQ=255)
                                                           gATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTcc                                                                                       >  3:848860/1‑149 (MQ=255)
                                                           gATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTcc                                                                                       >  1:770903/1‑149 (MQ=255)
                                                             tGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAg                                                                                     <  1:213134/149‑1 (MQ=255)
                                                              gCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGt                                                                                    >  2:774071/1‑149 (MQ=255)
                                                                aaCATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtc                                                                                  >  1:56387/1‑149 (MQ=255)
                                                                   aTGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTtcctc                                                                               >  1:835825/1‑149 (MQ=255)
                                                                              cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATCTTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGc                                                                    >  2:871454/1‑149 (MQ=255)
                                                                              cAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGc                                                                    >  3:884744/1‑149 (MQ=255)
                                                                               aGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCa                                                                   >  2:393007/1‑149 (MQ=255)
                                                                               aGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaa                                                                  >  2:595277/1‑149 (MQ=255)
                                                                                ggCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCTGTCAATCTCCAGTTCCTCCAGCCGACCGCaa                                                                  <  1:287091/149‑1 (MQ=255)
                                                                                    ttCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaa                                                              <  2:887152/149‑1 (MQ=255)
                                                                                     tCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaat                                                             <  3:523684/149‑1 (MQ=255)
                                                                                     tCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaat                                                             <  1:880991/149‑1 (MQ=255)
                                                                                      cACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaatg                                                            <  4:20560/149‑1 (MQ=255)
                                                                                      cACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatgaatg                                                            <  3:320460/149‑1 (MQ=255)
                                                                                       aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATgg                                                           <  3:818068/149‑1 (MQ=255)
                                                                                       aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATgg                                                           <  1:592405/149‑1 (MQ=255)
                                                                                       aCGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATgg                                                           <  3:47480/149‑1 (MQ=255)
                                                                                        cGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGc                                                          >  2:936066/1‑149 (MQ=255)
                                                                                        cGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGc                                                          >  3:565999/1‑149 (MQ=255)
                                                                                         gCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCt                                                         >  1:234731/1‑149 (MQ=255)
                                                                                          cAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTc                                                        <  1:393007/149‑1 (MQ=255)
                                                                                             cgcACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCt                                                     <  3:653069/149‑1 (MQ=255)
                                                                                               cACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGt                                                   >  3:719296/1‑149 (MQ=255)
                                                                                                    cccgccAGGGTCATACCCTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTGCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                              <  3:612496/149‑1 (MQ=255)
                                                                                                    cccgccAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTGCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCaa                                              >  2:65989/1‑149 (MQ=255)
                                                                                                         cAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgct                                         >  4:334741/1‑149 (MQ=255)
                                                                                                         cAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATtgct                                         <  2:770903/149‑1 (MQ=255)
                                                                                                              tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGc                                                                                                         <  3:874884/80‑1 (MQ=255)
                                                                                                              tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGc                                                                                                         >  4:874884/1‑80 (MQ=255)
                                                                                                              tCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAAc                                    >  4:452166/1‑149 (MQ=255)
                                                                                                                aTACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACtt                                  >  4:795800/1‑149 (MQ=255)
                                                                                                                 tACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACttt                                 >  2:254950/1‑149 (MQ=255)
                                                                                                                  acacTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACtttt                                >  4:907995/1‑149 (MQ=255)
                                                                                                                    acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                          <  3:810253/89‑1 (MQ=255)
                                                                                                                    acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                          >  4:887819/1‑89 (MQ=255)
                                                                                                                    acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                          <  3:887819/89‑1 (MQ=255)
                                                                                                                    acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAAtc                                                                                          >  4:810253/1‑89 (MQ=255)
                                                                                                                    acTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACACGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTcc                              <  3:761812/149‑1 (MQ=255)
                                                                                                                     cTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCa                             >  1:648245/1‑149 (MQ=255)
                                                                                                                      tGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAg                            >  1:927313/1‑149 (MQ=255)
                                                                                                                       gAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTACCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAgg                           <  3:775093/149‑1 (MQ=255)
                                                                                                                           atTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTaa                       <  3:101463/149‑1 (MQ=255)
                                                                                                                           atTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTaa                       <  4:884744/149‑1 (MQ=255)
                                                                                                                            ttCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAAt                      >  2:304568/1‑149 (MQ=255)
                                                                                                                             tCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATg                     >  1:489936/1‑149 (MQ=255)
                                                                                                                               aGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCg                   >  2:142698/1‑149 (MQ=255)
                                                                                                                                 ttCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGtt                 >  2:348319/1‑149 (MQ=255)
                                                                                                                                 ttCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGtt                 >  2:419105/1‑149 (MQ=255)
                                                                                                                                      tAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc            >  2:234772/1‑149 (MQ=255)
                                                                                                                                      tAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc            >  4:295372/1‑149 (MQ=255)
                                                                                                                                       aaaaaTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcca           >  1:559589/1‑149 (MQ=255)
                                                                                                                                           ctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCaa                                     <  1:65343/118‑1 (MQ=255)
                                                                                                                                           ctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCaa                                     >  2:65343/2‑119 (MQ=255)
                                                                                                                                           ctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc       >  2:626954/2‑149 (MQ=255)
                                                                                                                                           ctGGCGATCCTGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc       >  3:642717/2‑149 (MQ=255)
                                                                                                                                           aTGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc       <  4:432319/149‑1 (MQ=255)
                                                                                                                                               cGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgt   <  2:387597/149‑1 (MQ=255)
                                                                                                                                                gATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg  <  4:565999/149‑1 (MQ=255)
                                                                                                                                                 aTCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCatccatcc        >  3:906299/1‑142 (MQ=255)
                                                                                                                                                 aTCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCatccatcc        <  4:906299/142‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTG  >  NC_000913/3963268‑3963559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: