Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,974 A→T 6.3% S20T (TCG→ACG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9740AT 6.3% 306.8 / 6.7 118S20T (TCG→ACG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (66/43);  new base T (3/5);  total (70/48)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 122 total aligned reads displayed.

TGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTT  >  NC_000913/4336828‑4337122
                                                                                                                                                  |                                                                                                                                                    
tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac                                                                                                                                                   <  1:1462188/150‑1 (MQ=255)
  ggTCGAGGCACGACATTTTGGCGTATACCATCGAGAGCCCCAGGGCACCGATAATCGTCCCCAACGCTCCAGAAATGGCAAACCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGTAAAACACCGGACCCCATAATATTCCCCGacacaa                                                                                                                                                 >  1:854790/1‑148 (MQ=255)
  ggTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacaaca                                                                                                                                                 >  2:882717/1‑150 (MQ=255)
   gTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGCCCCCAAAATATTCCCCGCCACCAt                                                                                                                                                >  1:902444/1‑150 (MQ=255)
    tCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATc                                                                                                                                               <  1:1383813/150‑1 (MQ=255)
        ggAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                           >  1:512223/1‑150 (MQ=255)
          aaCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTg                                                                                                                                         >  2:1624340/1‑150 (MQ=255)
            cGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAc                                                                                                                                       >  1:1543576/1‑150 (MQ=255)
              aCATTTTGGCGTATACCATCGAGAGCCCCAGTGAACCGATAATCGTCACCAACCCTCCATAAATGGCAATCCCGCCACTAGCGGCCAGGTTTGCAGGTAACAGAAAAACACCAGACCCCATAATATTCCCCGACACCATCAGGGAGAAgg                                                                                                                                     >  1:1264174/1‑150 (MQ=255)
               cATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGAATATCGTCACCAACCATACATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCggg                                                                                                                                    <  1:879794/150‑1 (MQ=255)
                aTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGGTTGCAGGTAACAGAAAAACACCTGACCCCATAATATACCCCGTCACCCTGATGGTGTCgggg                                                                                                                                   <  1:133010/150‑1 (MQ=255)
                aTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACAATCCATAAATGGCAATCCCGCCAGTAGAGGACAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACAGCATCAGGGTGACgggg                                                                                                                                   >  1:1229326/1‑150 (MQ=255)
                aTTTTGGCGTATACCATCGAGAGCCCCAGTACACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                   >  2:964798/1‑150 (MQ=255)
                     ggCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCTACCATCCATATATGGCAATCCCGCCAGTAGAGTCCAGGTTTTCAGGATACAGAAAAACACCTGACCCCATAATATTCCCCGACACCTTCAGGGTGACGGGGATTaa                                                                                                                              <  1:1329502/150‑1 (MQ=255)
                           taCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGGCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                                                                                        >  1:528124/1‑150 (MQ=255)
                            aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                       >  1:1529351/1‑150 (MQ=255)
                            aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                       >  2:1454820/1‑150 (MQ=255)
                            aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                       <  1:1332693/150‑1 (MQ=255)
                             ccATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAAAGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGGAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCaggt                                                                                                                      >  1:1471178/1‑147 (MQ=255)
                              cATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACttt                                                                                                                     >  2:1035195/1‑150 (MQ=255)
                              cATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGGGTGCAGGTAACAGAAAACCCCCTGACCCCATAACCTTCCCCGACCCCATCAGGGTGCCGGGGATTACGcccccgtt                                                                                                                     >  1:254731/1‑145 (MQ=255)
                                 cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                  <  2:1646532/150‑1 (MQ=255)
                                  gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCCAACCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATAAAGCCCACTTtgtta                                                                                                                 >  1:1053349/1‑148 (MQ=255)
                                  gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGa                                                                                                                 >  1:1160858/1‑150 (MQ=255)
                                      gCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcat                                                                                                             >  1:1728418/1‑150 (MQ=255)
                                       ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGATAAACCCCTGACCCCATAATATTCCCCGACAACATCAGGGTGACGGGGATTAAGCCCACTATGTGagcaac                                                                                                            >  1:278165/1‑148 (MQ=255)
                                       ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                            >  1:1573075/1‑150 (MQ=255)
                                        cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatca                                                                                                           <  2:763113/150‑1 (MQ=255)
                                         ccAGTGCACCGATAATCGGCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACGTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcag                                                                                                          <  1:1060507/150‑1 (MQ=255)
                                               cACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                    >  2:1049985/1‑150 (MQ=255)
                                               cACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                    >  2:687579/1‑150 (MQ=255)
                                               cACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACGTTGTGAGcatccgcatccg                                                                                                    >  2:247485/1‑150 (MQ=255)
                                                  cGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                 <  2:131247/150‑1 (MQ=255)
                                                   gATAATCGTCACCACCCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGATAAACCCCTGCCCCCCTAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCAGCCGCATCCGAAgc                                                                                                >  2:309482/1‑149 (MQ=255)
                                                          gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                         <  1:1659266/150‑1 (MQ=255)
                                                           tCACCAACCATCCATAAATGGCAATCCAGCCAGTAGAGGCCAGGTTTGCAGGTCACTGATAAACAACTGACCCCTTAATATTCCCCGACACCATCTGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                        <  1:346172/150‑1 (MQ=255)
                                                            cACCAACAATCCATAAATGGCAATTCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                       >  1:502387/1‑150 (MQ=255)
                                                              ccaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAc                                                                                     <  1:1520252/150‑1 (MQ=255)
                                                               caaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt                                                                                    >  1:395854/1‑150 (MQ=255)
                                                               caaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt                                                                                    <  2:1358341/150‑1 (MQ=255)
                                                                    atccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCGGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                               >  2:1725376/1‑150 (MQ=255)
                                                                    atccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATAATCCCCGACACCATCAGGGTGACGGGGATTCAGCCCACTTTGAGAGCATCAGCATAAGAAGACATAATTAAACTCCTGc                                                                               >  2:1068189/1‑150 (MQ=255)
                                                                    atccatAAATGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCAGACCCCATAATATTCCCCGACACCATCAGCGTGACTGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                               <  1:256607/150‑1 (MQ=255)
                                                                      ccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                             >  2:755623/1‑150 (MQ=255)
                                                                       catAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGGCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGaa                                                                            >  2:555709/1‑150 (MQ=255)
                                                                        atAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAg                                                                           <  1:895890/150‑1 (MQ=255)
                                                                          aaaTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGc                                                                         <  2:338876/150‑1 (MQ=255)
                                                                          aaaTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGc                                                                         <  1:1430186/150‑1 (MQ=255)
                                                                             tGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                      <  1:1049985/150‑1 (MQ=255)
                                                                                cAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                   >  2:1482984/1‑150 (MQ=255)
                                                                                cAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                   <  2:512223/150‑1 (MQ=255)
                                                                                   tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                <  1:1035195/150‑1 (MQ=255)
                                                                                    cccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                               >  2:1248708/1‑150 (MQ=255)
                                                                                     ccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCAGCCGAAGACATAATTAAACTCCTGCGAAGGCGCGCTTCGTGa                                                              >  2:1076239/1‑150 (MQ=255)
                                                                                      cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                             <  1:22023/150‑1 (MQ=255)
                                                                                      cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                             <  2:1165661/150‑1 (MQ=255)
                                                                                       gCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                            <  2:421294/150‑1 (MQ=255)
                                                                                           gTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                        >  1:670434/1‑150 (MQ=255)
                                                                                              gagGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                     >  2:347228/1‑150 (MQ=255)
                                                                                                  ccAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                 <  1:572369/150‑1 (MQ=255)
                                                                                                  ccAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                 >  1:643575/1‑150 (MQ=255)
                                                                                                   cAGGTTTGCAGGTAACAGAAAAACACCTGCCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                                >  2:1330510/1‑150 (MQ=255)
                                                                                                   cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                                <  2:1248930/150‑1 (MQ=255)
                                                                                                        tggCAGGTAACAGAAAACCACCTGACCCCATAACATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatc                                           <  2:254731/148‑1 (MQ=255)
                                                                                                           cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                        >  2:613074/1‑150 (MQ=255)
                                                                                                             ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACg                                      >  2:482349/1‑150 (MQ=255)
                                                                                                              gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                     >  2:944428/1‑150 (MQ=255)
                                                                                                               tAACAGAAAAACACCTGACCCCATTATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                    >  1:254846/1‑150 (MQ=255)
                                                                                                                aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                   >  2:520100/1‑150 (MQ=255)
                                                                                                                aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAAGGGGATTAAGCCCACTTTGTGAGCATCTGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                   >  2:1097760/1‑150 (MQ=255)
                                                                                                                 aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGTGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                  <  2:1160858/150‑1 (MQ=255)
                                                                                                                 aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                  <  1:1194477/150‑1 (MQ=255)
                                                                                                                  cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                                 >  1:1070983/1‑150 (MQ=255)
                                                                                                                  cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                                 >  2:200965/1‑150 (MQ=255)
                                                                                                                  cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACTCCTGc                                 >  1:1615489/1‑150 (MQ=255)
                                                                                                                    gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                               >  2:1539497/1‑150 (MQ=255)
                                                                                                                     aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGCCAATAAATTACGTCATATCATACGCCGGCAtt                              >  2:266722/1‑150 (MQ=255)
                                                                                                                     aaaaaCACCTGACCCCATAATATTCCCAGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAATGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                              <  2:1229326/150‑1 (MQ=255)
                                                                                                                         acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                          >  1:204011/1‑150 (MQ=255)
                                                                                                                           acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                        >  1:148122/1‑150 (MQ=255)
                                                                                                                           aacTGACCCCATAATATTCCACGACACCATCAGGGTGACGGGGATTAAGCCGACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                        <  1:132480/148‑1 (MQ=255)
                                                                                                                              tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                     <  1:347228/150‑1 (MQ=255)
                                                                                                                              tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                     <  1:1666751/150‑1 (MQ=255)
                                                                                                                                aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                   <  2:1728418/150‑1 (MQ=255)
                                                                                                                                 ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                  >  2:104642/1‑150 (MQ=255)
                                                                                                                                 ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                  >  1:519699/1‑150 (MQ=255)
                                                                                                                                  cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                 <  1:555709/150‑1 (MQ=255)
                                                                                                                                  cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                 <  2:902444/150‑1 (MQ=255)
                                                                                                                                   ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGCGCATCCGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAAGCATGCATTCGATAGTTGAattt                >  2:1304554/1‑149 (MQ=255)
                                                                                                                                   ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatta                >  2:1205571/1‑150 (MQ=255)
                                                                                                                                     ataataTTCCCCGACCCCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCATCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt              <  2:1530198/150‑1 (MQ=255)
                                                                                                                                        atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattata           >  1:1178564/1‑149 (MQ=255)
                                                                                                                                        atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc           >  2:1304935/1‑150 (MQ=255)
                                                                                                                                         tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATcc          <  2:502387/150‑2 (MQ=255)
                                                                                                                                           atCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGACATATCAGACGCCTGCCTTCGATACTATAATTATTAACGaa        >  1:285310/2‑149 (MQ=255)
                                                                                                                                             ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg      >  2:1664915/1‑150 (MQ=255)
                                                                                                                                              cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATAGACCgg     >  2:334576/1‑150 (MQ=255)
                                                                                                                                              cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAGAGTTGAATTATTATCGCCCgg     >  2:1139392/1‑150 (MQ=255)
                                                                                                                                               ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg    >  1:967590/1‑150 (MQ=255)
                                                                                                                                                 gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGtt  >  1:454121/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTT  >  NC_000913/4336828‑4337122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: