Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,082,418 A→T 5.2% N59I (AAC→ATC)  efeO → inactive ferrous ion transporter EfeUOB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,082,4180AT 5.2% 238.9 / 2.9 77N59I (AAC→ATC) efeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base T (2/2);  total (34/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 = 8.15e-01

TTGCAGTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAC  >  NC_000913/1082270‑1082562
                                                                                                                                                    |                                                                                                                                                
ttGCAGTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGaa                                                                                                                                                  <  2:26533/149‑1 (MQ=255)
   cAGTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACca                                                                                                                                               <  5:245026/149‑1 (MQ=255)
   cAGTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACca                                                                                                                                               <  8:19762/149‑1 (MQ=255)
     gTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCATTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACcaca                                                                                                                                             <  7:95434/149‑1 (MQ=255)
      ttGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGaa                                                                                                                                                  <  6:113043/143‑1 (MQ=255)
      ttGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGaa                                                                                                                                                  >  5:113043/1‑143 (MQ=255)
        gAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGc                                                                                                                                           <  2:215126/148‑1 (MQ=255)
        gAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGc                                                                                                                                           <  5:182105/148‑1 (MQ=255)
                   cgcTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGgcgc                                                                                                                                 <  5:29394/147‑1 (MQ=255)
                        ttttCTTCTGGTTTTATGGCTAACGCCGCTGATGTGCAGCAGGTCAAAGTGACCGTGTCGGTTTAGCAGTGCGATGCGATGACCATTACGGTTATCGCCGGGAAAACACAGTTCATTATTCAGTACGACAGCCCGACGGCGCTGGAGt                                                                                                                           <  4:243919/148‑1 (MQ=255)
                        ttttCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAg                                                                                                                            <  3:120823/147‑1 (MQ=255)
                          ttcttcTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTgg                                                                                                                         <  2:236057/148‑1 (MQ=255)
                           tcttctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAAc                                                                                                                                                 <  1:233636/123‑1 (MQ=255)
                           tcttctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAAc                                                                                                                                                 >  2:233636/1‑123 (MQ=255)
                           tcttctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTggg                                                                                                                        <  5:307151/148‑1 (MQ=255)
                            cttctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTgg                                                                                                                         >  7:370520/1‑146 (MQ=255)
                             ttctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTgg                                                                                                                         >  2:349124/1‑145 (MQ=255)
                              tctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTgg                                                                                                                         >  1:233306/1‑144 (MQ=255)
                              tctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTgg                                                                                                                         >  4:232810/1‑144 (MQ=255)
                               ctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATc                                                                                                                   >  6:343180/1‑149 (MQ=255)
                                tGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATcc                                                                                                                  >  8:38372/1‑149 (MQ=255)
                                tGCTTTTATGGATAACTCCGCTGATGTGCCGCAGGTCAAATTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCAAAGCCAGAAGGCGCTGGAGTCGCAGATc                                                                                                                   >  4:36374/1‑148 (MQ=255)
                                  cTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATc                                                                                                                   <  1:244229/146‑1 (MQ=255)
                                   ttttATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTc                                                                                                                >  7:52262/1‑148 (MQ=255)
                                    tttATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCa                                                                                                               >  3:10167/1‑148 (MQ=255)
                                      tATGGCTAACTCCGCTGATGTGCAGCAGGTCAAAGTGAACGTGACGGATAAGCAGTGCGATCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACATCCAGAAGTCGCTGTAGTGGGAGATCCTCaaa                                                                                                             <  8:192156/148‑1 (MQ=255)
                                       aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAg                                                                                                            <  7:38372/148‑1 (MQ=255)
                                          gCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTggg                                                                                                                        >  7:79874/1‑133 (MQ=255)
                                          gCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTggg                                                                                                                        <  8:79874/133‑1 (MQ=255)
                                           cTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTg                                                                                                       <  4:214378/149‑1 (MQ=255)
                                            tAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGa                                                                                                      >  6:68415/1‑149 (MQ=255)
                                              aCGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtg                                                                                                    <  4:10167/149‑1 (MQ=255)
                                                 ccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtg                                                                                                 <  8:370520/149‑1 (MQ=255)
                                                   gcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAACCACAGTTCATTATTCAGCACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtggt                                                                                               >  5:214104/1‑149 (MQ=255)
                                                   gcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGGTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtggt                                                                                               >  6:159914/1‑149 (MQ=255)
                                                    cTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGAGCATTACGGTTAACGCCGGGATAACACAGTTCATTATTCAGAACCACAGCCTGAAGGCGCTGGAGTGGTAGATCCTCAAAGGCGTGAtggtggtg                                                                                              <  5:225363/149‑1 (MQ=255)
                                                           gCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTCCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGAGGTGGAAGAGc                                                                                       >  3:53117/1‑149 (MQ=255)
                                                             cGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCgg                                                                                     >  2:375544/1‑149 (MQ=255)
                                                                   tCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGaaaa                                                                                >  5:68741/1‑148 (MQ=255)
                                                                      aaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATc                                                                            >  7:50465/1‑149 (MQ=255)
                                                                         tGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcc                                                                         <  5:343180/149‑1 (MQ=255)
                                                                               tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACATTTCATTATTCAGAAACACAGCCAGAAGGCGCTGGAGTTGGAGATCCTCAAAGGCGGGATGGTGGTGTAAGAGCGGGAAAATATCGCCCCTGGc                                                                   <  1:83770/149‑1 (MQ=255)
                                                                               tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGATCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGc                                                                   >  6:209864/1‑149 (MQ=255)
                                                                               tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGATCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGc                                                                   >  8:102837/1‑149 (MQ=255)
                                                                                 aCGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCtt                                                                 <  3:232810/149‑1 (MQ=255)
                                                                                  cGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                <  8:50465/149‑1 (MQ=255)
                                                                                         gCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAg                                                          >  5:88344/1‑148 (MQ=255)
                                                                                         gCAGTGCGAACCGATGACCATTACGGTTAACGACGGGAAAACACAGTTCATTATTCAGAACCACAGCAAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTTGTGGAAGATCTGGAAAATATc                                                                            >  1:231944/1‑130 (MQ=255)
                                                                                          cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTAAAGCCAggg                                                        <  7:398240/149‑3 (MQ=255)
                                                                                            gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaaa                                                      <  7:139158/149‑1 (MQ=255)
                                                                                               cGACCCGCTGACCATTACGGGTAACGCCGGGACAACACACCTCATTATTCAGAACCACAGCCAGATGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGCGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                                     <  6:385514/147‑1 (MQ=255)
                                                                                               cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGa                                                                                            <  4:75940/108‑1 (MQ=255)
                                                                                               cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGa                                                                                            >  3:75940/1‑108 (MQ=255)
                                                                                               cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGa                                                   <  1:375544/149‑1 (MQ=255)
                                                                                                gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                  >  2:227642/1‑149 (MQ=255)
                                                                                                gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                  <  6:68741/149‑1 (MQ=255)
                                                                                                 aaCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                 <  2:233306/149‑1 (MQ=255)
                                                                                                 aaCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTGTAGCCAGAAAATGACg                                                 >  4:306312/1‑149 (MQ=255)
                                                                                                       tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAt                                           <  1:227642/149‑1 (MQ=255)
                                                                                                       tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAt                                           <  1:349124/149‑1 (MQ=255)
                                                                                                              tACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                                    >  6:339420/1‑149 (MQ=255)
                                                                                                                  gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCa                                                           <  7:190329/122‑1 (MQ=255)
                                                                                                                  gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCa                                                           >  8:190329/1‑122 (MQ=255)
                                                                                                                   ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACATCCAGAAGGAGCTGGAGTGGGAGATCCTCAAATGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCATAAAATGACGGCGAATTTACAGCCTGGc                               <  3:36374/149‑1 (MQ=255)
                                                                                                                   ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                               <  8:341869/149‑1 (MQ=255)
                                                                                                                   ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                               >  1:348545/1‑149 (MQ=255)
                                                                                                                   ttAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGc                               <  3:306312/149‑1 (MQ=255)
                                                                                                                            ggAAAACACAGTTCATTATTCAGATCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat                      <  5:209864/149‑1 (MQ=255)
                                                                                                                            ggAAAACACAGTTCATTATTCAGATCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat                      <  7:102837/149‑1 (MQ=255)
                                                                                                                            ggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat                      <  5:68415/149‑1 (MQ=255)
                                                                                                                            ggAAAACACAGTTCATTATTCAGAAAAACAGCAAGAAGGCGATGGATTGGTATCTCATCAAATGATAGATGGTGTAGTATGAGAGGGAAAATATAGCCCCTTGATTTAGCCAGAAAATGACGTCGAATTTACAGACTGGAGAATACGa                       >  7:284400/1‑148 (MQ=255)
                                                                                                                             gAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGata                     <  8:52262/149‑1 (MQ=255)
                                                                                                                             gAAAACACAGGTCATTAGTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGATTTTACAGCCTGGCGAATACGata                     <  1:26209/149‑1 (MQ=255)
                                                                                                                                 acacAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGAc                 >  6:182907/1‑149 (MQ=255)
                                                                                                                                   acaGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGTCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCt               >  2:303883/1‑149 (MQ=255)
                                                                                                                                      cttCATTATTCAGAACCACAGCCAGAAGCCGCTGGAGTGGGAGATCCTCACAGGCGTCATCGTGGTGGAAGACCGGGAAAATATCGCCCCTGGCGTTAGCCAGAAAATGACGGCGAATATACAGCCTTGCGAATACGATATGACCTGCg            <  5:159914/148‑1 (MQ=255)
                                                                                                                                       ttCATTATTCAGAACCCCAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCgg           >  3:345698/1‑149 (MQ=255)
                                                                                                                                           ttattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATa                          <  3:365784/130‑1 (MQ=255)
                                                                                                                                           ttattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATa                          >  4:365784/1‑130 (MQ=255)
                                                                                                                                                cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAc  >  4:31864/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTGCAGTTGAGCGTGGCTGCGCTGTTTTCTTCTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAC  >  NC_000913/1082270‑1082562

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: