Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,134,556 G→T 6.1% L142L (CTG→CTT yghT → putative ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,134,5560GT 6.1% 289.9 / 7.3 82L142L (CTG→CTTyghTputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base T (2/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GCTTACTGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATAACGCCA  >  NC_000913/3134408‑3134703
                                                                                                                                                    |                                                                                                                                                   
gCTTACTGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCAtggctg                                                                                                                                                     <  1:453575/149‑1 (MQ=255)
    aCTGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCaa                                                                                                                                                 >  2:469881/1‑149 (MQ=255)
    aCTGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCaa                                                                                                                                                 >  1:388860/1‑149 (MQ=255)
      tGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGc                                                                                                                                               >  4:247277/1‑149 (MQ=255)
           tttCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTTGCCAGCCATGa                                                                                                                                          <  2:211334/149‑1 (MQ=255)
            ttCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGaa                                                                                                                                         <  4:562218/149‑1 (MQ=255)
               gTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGt                                                                                                                                      >  2:503522/1‑149 (MQ=255)
                tCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                                                                                                     >  1:144162/1‑149 (MQ=255)
                tCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                                                                                                     >  2:530546/1‑149 (MQ=255)
                tCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                                                                                                     <  3:465664/149‑1 (MQ=255)
                   tatGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAg                                                                                                                                  >  4:396717/1‑149 (MQ=255)
                   tatGCTTGCCTGTCGTCGTCGCGGCTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAg                                                                                                                                  <  2:388860/149‑1 (MQ=255)
                      gCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                                                                                               <  3:521765/149‑1 (MQ=255)
                       cTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGtt                                                                                                                              >  2:99955/1‑149 (MQ=255)
                        ttGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGttt                                                                                                                             >  3:574811/1‑149 (MQ=255)
                            cTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACt                                                                                                                         >  2:512727/1‑149 (MQ=255)
                             tGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTg                                                                                                                        >  4:282664/1‑149 (MQ=255)
                                 gtcgtcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATg                                                                                                                    >  3:76832/1‑149 (MQ=255)
                                        cgGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                                                                                                             <  1:469881/149‑1 (MQ=255)
                                              aaTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGc                                                                                                       <  2:5529/149‑1 (MQ=255)
                                               aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                                                                                                      >  2:526557/1‑149 (MQ=255)
                                               aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                                                                                                      >  4:451441/1‑149 (MQ=255)
                                                   ttCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGAt                                                                                                  >  3:486677/1‑149 (MQ=255)
                                                    tCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATc                                                                                                 <  1:424391/149‑1 (MQ=255)
                                                       aaCCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTgg                                                                                              <  4:222470/149‑1 (MQ=255)
                                                        accgaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGt                                                                                             >  2:285586/1‑149 (MQ=255)
                                                        accgaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGt                                                                                             >  3:352604/1‑149 (MQ=255)
                                                         ccgaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtc                                                                                            <  3:396717/149‑1 (MQ=255)
                                                           gaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                                                                                                             >  3:210160/1‑130 (MQ=255)
                                                           gaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                                                                                                             <  4:210160/130‑1 (MQ=255)
                                                           gaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcat                                                                                          >  4:5810/1‑149 (MQ=255)
                                                            accgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatc                                                                                         >  4:129832/1‑149 (MQ=255)
                                                              cgTTATCCTCAGGACCAAATTCCTGTCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTTGCCAGCCATGa                                                                                                                                          <  4:479264/98‑1 (MQ=255)
                                                              cgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTTGCCAGCCATGa                                                                                                                                          >  3:479264/1‑98 (MQ=255)
                                                               gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATCAGCCTGATCTGGTCATCaaa                                                                                      <  1:474685/149‑1 (MQ=255)
                                                               gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCaaa                                                                                      <  2:107955/149‑1 (MQ=255)
                                                                ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc                                                                                     <  2:201409/149‑1 (MQ=255)
                                                                ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc                                                                                     >  2:363252/1‑149 (MQ=255)
                                                                ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc                                                                                     >  1:588751/1‑149 (MQ=255)
                                                                   tCCTCAGGCCCAAATTCCTGGCGCTTACGATGGTACGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCaa                                                                                 <  3:212019/149‑1 (MQ=255)
                                                                     cTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAAt                                                                                <  3:264628/149‑1 (MQ=255)
                                                                     cTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAAt                                                                                <  3:268940/149‑1 (MQ=255)
                                                                     cTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAAt                                                                                <  4:28861/149‑1 (MQ=255)
                                                                      tCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATg                                                                               >  4:23904/1‑149 (MQ=255)
                                                                          gACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGa                                                                           <  1:526557/149‑1 (MQ=255)
                                                                           aCCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGAc                                                                          <  3:349256/149‑1 (MQ=255)
                                                                            ccAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGAcc                                                                         >  2:461370/1‑149 (MQ=255)
                                                                            ccAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGAcc                                                                         >  4:205066/1‑149 (MQ=255)
                                                                               aaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTg                                                                      >  4:71208/1‑149 (MQ=255)
                                                                                 ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGaa                                                                    <  3:8496/149‑1 (MQ=255)
                                                                                   ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGt                                                                  <  1:363252/149‑1 (MQ=255)
                                                                                    cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTTgccagccatga                                                                                                                                          >  3:491115/1‑73 (MQ=255)
                                                                                    cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTTGCCAGCCATGa                                                                                                                                          <  4:491115/76‑1 (MQ=255)
                                                                                        cgcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTTGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCt                                                             <  4:375304/149‑1 (MQ=255)
                                                                                        cgcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCt                                                             <  1:521568/149‑1 (MQ=255)
                                                                                        cgcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCt                                                             <  4:199440/149‑1 (MQ=255)
                                                                                          cTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCtgt                                                           <  1:461370/149‑1 (MQ=255)
                                                                                            tACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGc                                                         <  1:503522/149‑1 (MQ=255)
                                                                                              cGATTGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTTGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCAc                                                       <  4:339143/149‑1 (MQ=255)
                                                                                              cGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCAc                                                       <  3:282664/149‑1 (MQ=255)
                                                                                              cGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCAc                                                       <  4:76832/149‑1 (MQ=255)
                                                                                                    tACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAAc                                                 >  2:73972/1‑149 (MQ=255)
                                                                                                     aCGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAAcc                                                >  4:99655/1‑149 (MQ=255)
                                                                                                       ggTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCa                                                                                  <  1:353833/113‑1 (MQ=255)
                                                                                                       ggTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCa                                                                                  >  2:353833/1‑113 (MQ=255)
                                                                                                       ggTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCCGTCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCg                                              >  4:359510/1‑149 (MQ=255)
                                                                                                              ccACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATaa                                       >  1:539897/1‑149 (MQ=255)
                                                                                                                aCCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGc                                                               <  2:45494/123‑1 (MQ=255)
                                                                                                                aCCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGc                                                               >  1:45494/1‑123 (MQ=255)
                                                                                                                aCCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAAc                                     <  2:100186/149‑1 (MQ=255)
                                                                                                                  cTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACgg                                   <  3:5810/149‑1 (MQ=255)
                                                                                                                   tAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGAcc                                           >  3:329492/1‑140 (MQ=255)
                                                                                                                   tAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGAcc                                           <  4:329492/140‑1 (MQ=255)
                                                                                                                   tAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAACGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACggg                                  >  1:175070/1‑149 (MQ=255)
                                                                                                                     aCGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGaa                                <  3:451441/149‑1 (MQ=255)
                                                                                                                          gAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCt                           <  2:464759/149‑1 (MQ=255)
                                                                                                                                      tttGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGAt               <  3:71208/149‑1 (MQ=255)
                                                                                                                                         gTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGc            <  1:99955/149‑1 (MQ=255)
                                                                                                                                          tctcATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGcc           <  3:99655/149‑1 (MQ=255)
                                                                                                                                            tcATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCAt         >  1:350226/1‑149 (MQ=255)
                                                                                                                                            tcATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCAt         >  1:154341/1‑149 (MQ=255)
                                                                                                                                              aTGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATaa       >  2:376435/1‑149 (MQ=255)
                                                                                                                                              aTGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATaa       <  1:512727/149‑1 (MQ=255)
                                                                                                                                                   tggcAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATAACGCca  >  4:504950/1‑149 (MQ=255)
                                                                                                                                                    ggcAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCATGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATAAc      >  1:103019/1‑144 (MQ=255)
                                                                                                                                                    ggcAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATAAc      <  2:103019/144‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTTACTGCGCTTTCGTCATATGCTTGCCTGTCGTCGTCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAATGTTGACCTTGAAGTTGCCTGTGCACGTAAACCCGACCATAAACGGGAATCGCTGGCGAGGAAGATTGCCATAACGCCA  >  NC_000913/3134408‑3134703

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