Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,679 G→T 7.5% intergenic (+68/‑779) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,6790GT 7.5% 261.1 / 9.8 80intergenic (+68/‑779)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (30/43);  new base T (3/3);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTT  >  NC_000913/770537‑770823
                                                                                                                                              |                                                                                                                                                
gAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                                                                             <  1:830166/148‑1 (MQ=255)
 aaaaTCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACgg                                                                                                                                            <  2:553782/148‑1 (MQ=255)
     tCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGa                                                                                                                                       >  3:355675/1‑149 (MQ=255)
      ctaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGaa                                                                                                                                      <  3:236414/147‑1 (MQ=255)
      cAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGaa                                                                                                                                      <  4:149532/149‑1 (MQ=255)
       aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                                                          >  1:433138/1‑144 (MQ=255)
       aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                                                          <  2:433138/144‑1 (MQ=255)
       aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAc                                                                                                                                     <  4:370555/149‑1 (MQ=255)
       aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGaa                                                                                                                                      <  2:913113/148‑1 (MQ=255)
         ggTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAAGGCATGAACGc                                                                                                                                   <  4:89739/149‑1 (MQ=255)
         ctctgcaaTCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACg                                                                                                                                    >  2:918109/7‑148 (MQ=255)
          gTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCt                                                                                                                                  <  4:260294/149‑1 (MQ=255)
            tCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCttt                                                                                                                                <  4:292345/149‑1 (MQ=255)
              aaaTCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCGCATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAt                                                                                                                              >  4:815670/1‑149 (MQ=255)
                aTCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAtt                                                                                                                             <  2:542329/148‑1 (MQ=255)
                  cTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                                                                           <  1:563333/148‑1 (MQ=255)
                  cTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCg                                                                                                                           <  2:28262/148‑1 (MQ=255)
                   tATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCCTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                                                          <  1:588136/148‑1 (MQ=255)
                    aTCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATATCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                                                                         >  1:791143/1‑148 (MQ=255)
                    aTCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                                                                         <  1:795628/148‑1 (MQ=255)
                     tCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTTAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                        <  1:184998/148‑1 (MQ=255)
                           ccTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaaa                                                                                                                 <  4:408375/149‑1 (MQ=255)
                           ccTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaa                                                                                                                  >  2:313392/1‑148 (MQ=255)
                             ttttttACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                                                               >  3:422622/1‑149 (MQ=255)
                              attttACCCGCCCACTCCCGGATGTTCAGTGACCGGCTTGTCTTAATAGCAATCAAACCGAAGCTAAATAGGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAATGCTTTATTCGGTCTACAAAAgga                                                                                                              <  4:752160/148‑3 (MQ=255)
                                     ccGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAAGGCATGAACGCTTTACTCGGTCTACAAAAGCAGGCaaa                                                                                                        <  2:481623/148‑1 (MQ=255)
                                       ggCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                                                                      <  1:106769/148‑1 (MQ=255)
                                        gCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTc                                                                                                     >  2:742748/1‑148 (MQ=255)
                                         ccAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                   <  4:355675/149‑1 (MQ=255)
                                         ccAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCa                                                                                                    <  2:711468/148‑1 (MQ=255)
                                         ccAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCAGTCTACAAAAGCAGGCAAATTCaa                                                                                                   <  3:545227/149‑1 (MQ=255)
                                              cACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaa                                                                                                                  <  1:437826/129‑1 (MQ=255)
                                              cACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaa                                                                                                                  >  2:437826/1‑129 (MQ=255)
                                              cACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCAACATATGCGCGGCCAGATTGTTGACCAAGGGCGCTTTGTTCATGCCGGATACGGCTTGAACGCATTATTCGGTCTACAAAAGCAGGCAAATTCAatata                                                                                               >  1:234816/1‑148 (MQ=255)
                                              cACGGACGATCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAatata                                                                                               <  2:791143/148‑1 (MQ=255)
                                                      ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                                                                       <  1:918109/148‑1 (MQ=255)
                                                        cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                     >  2:271046/1‑148 (MQ=255)
                                                         aGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATTCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgatga                                                                                   <  3:169149/149‑1 (MQ=255)
                                                               cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                     <  3:92174/141‑1 (MQ=255)
                                                               cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                     >  4:92174/1‑141 (MQ=255)
                                                                     gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTg                                                                        <  1:137702/148‑1 (MQ=255)
                                                                        tGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCa                                                                     >  1:785696/1‑148 (MQ=255)
                                                                            tAGCAATCAAACCGAAGCCACATATGCGCTGCCAGATTGTGGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTATACACTAGCAAGCGt                                                                >  4:484409/1‑149 (MQ=255)
                                                                            tAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCg                                                                 >  2:11121/1‑148 (MQ=255)
                                                                             aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                                                                >  1:289530/1‑148 (MQ=255)
                                                                             aGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGt                                                                >  1:112035/1‑148 (MQ=255)
                                                                                 aTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                                                                  <  1:770630/62‑1 (MQ=255)
                                                                                 aTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                                                                  >  2:770630/1‑62 (MQ=255)
                                                                                  tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGAATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCa                                                          >  3:620680/1‑149 (MQ=255)
                                                                                    aaaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATc                                                        >  3:662468/1‑149 (MQ=255)
                                                                                     aaCCGAAGCCACATATGCTCGGCCAGATTGTTCACAAAGGGCGCTCTGTTCATGCTGGATACGGCATGAACGTTTTATGCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATc                                                        <  2:234816/148‑1 (MQ=37)
                                                                                      aCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTCGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg                                                      >  4:205867/1‑149 (MQ=255)
                                                                                          aaGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                                                             >  3:615526/1‑138 (MQ=255)
                                                                                          aaGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                                                             <  4:615526/138‑1 (MQ=255)
                                                                                              cacaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAttt                                               >  2:113571/1‑148 (MQ=255)
                                                                                               acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAtttt                                              <  1:363874/148‑1 (MQ=255)
                                                                                               acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTg                                             >  3:305755/1‑149 (MQ=255)
                                                                                                caTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGGGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTg                                             >  1:516000/1‑148 (MQ=255)
                                                                                                 atatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGc                                            <  1:909443/148‑1 (MQ=255)
                                                                                                 atatGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCg                                           <  3:767641/149‑1 (MQ=255)
                                                                                                               aTTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCActct                             >  3:895915/1‑149 (MQ=255)
                                                                                                                      aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGcc                       >  1:445893/1‑148 (MQ=255)
                                                                                                                      aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGcc                       >  1:84080/1‑148 (MQ=255)
                                                                                                                       cAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCa                      >  1:868482/1‑148 (MQ=255)
                                                                                                                         aaGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACg                    >  1:588264/1‑148 (MQ=255)
                                                                                                                          aGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                  >  4:932709/1‑149 (MQ=255)
                                                                                                                            ggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGtata                >  3:236187/1‑149 (MQ=255)
                                                                                                                              cgcTTTGTGCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat               <  1:737665/148‑1 (MQ=255)
                                                                                                                                cTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt             <  1:664264/148‑1 (MQ=255)
                                                                                                                                   tGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGGGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGc          <  2:516000/148‑1 (MQ=255)
                                                                                                                                   tGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGc          <  1:857635/148‑1 (MQ=255)
                                                                                                                                     ttCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCa                                 >  3:930696/1‑123 (MQ=255)
                                                                                                                                     ttCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCa                                 <  4:930696/123‑1 (MQ=255)
                                                                                                                                      tCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCttt       >  1:175596/1‑148 (MQ=255)
                                                                                                                                       cATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTa      >  1:240651/1‑148 (MQ=255)
                                                                                                                                       cATGCCGGATACGTCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCTAGCCACGTATATGTGGCTTTAt     <  4:88632/149‑1 (MQ=255)
                                                                                                                                       cATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTa      <  2:588264/148‑1 (MQ=255)
                                                                                                                                        aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAGGCAGGCAAATTCAATATATTGCAGAGGTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACGCTCAAGCCACGTATGCGTGGCttgag     <  2:265237/148‑4 (MQ=37)
                                                                                                                                         tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt    <  1:84825/148‑1 (MQ=255)
                                                                                                                                          gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt  <  3:576333/149‑1 (MQ=255)
                                                                                                                                           ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt  <  2:289530/148‑1 (MQ=255)
                                                                                                                                           ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt  <  2:112184/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTT  >  NC_000913/770537‑770823

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