Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,974 A→T 6.7% 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.7% 276.8 / 7.8 99S20T (TCG→ACG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base T (2/5);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGT  >  NC_000913/4336831‑4337121
                                                                                                                                               |                                                                                                                                                   
gTCGAGGAACGACAGTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACTATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAATCACCTGCCCCCATAATATTCACAGACACCAt                                                                                                                                               >  2:325726/1‑150 (MQ=255)
  cGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                                                                                             >  1:1070771/1‑150 (MQ=255)
     ggtaCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                          <  1:1156340/147‑1 (MQ=255)
     ggAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                          <  1:491217/150‑1 (MQ=255)
      gAACGACATTTTGGCGTATACCCTCGAGAGCCCCAGTGCACCGATTATCGTCCCCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTAGCAGGTAACAGACCGGCCCCGGCCCCCATAATATTCCCCGCCCCCATCAgggg                                                                                                                                         >  2:305483/1‑149 (MQ=255)
        aCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCAAGGTTTGCATGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACCCCATCAGGGTGa                                                                                                                                       >  2:359147/1‑150 (MQ=255)
         cGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAc                                                                                                                                      >  1:361589/1‑150 (MQ=255)
           aCATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgg                                                                                                                                    <  1:237593/150‑1 (MQ=255)
           aCATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCACCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgg                                                                                                                                    <  2:377306/150‑1 (MQ=255)
            cATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCAATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGCCCCCATAATATTCCCCGACACCATCAGGGTGACggg                                                                                                                                   >  1:119210/1‑150 (MQ=255)
             aTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                  <  1:1344045/150‑1 (MQ=255)
             aTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                  <  2:821493/150‑1 (MQ=255)
              ttttGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCGGGGa                                                                                                                                 <  1:397816/150‑1 (MQ=255)
              ttttGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGa                                                                                                                                 <  1:540114/150‑1 (MQ=255)
               tttGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                <  1:1426253/150‑1 (MQ=255)
                   gCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCCCCAACCATCCATAAATGGCAATCCCCCCAGTACACGCCAGGTTCGCAGGTAGCAGAAAAACACCAGACCCCAAAATATTCCCCGACACCATCAGGGGGACGGGGATTAAg                                                                                                                            >  2:1412149/1‑150 (MQ=255)
                   gCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAg                                                                                                                            <  2:471799/150‑1 (MQ=255)
                    cGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTATCAGAAAAACACCTGACCCCAGAATATTCCCCGACACCATCAGGGTGACGGGGTTTAAGc                                                                                                                           >  1:1723039/1‑150 (MQ=255)
                         aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                      <  1:1386926/150‑1 (MQ=255)
                         aCCATCCAGAGCCCCAGTGCACCGAAAATCGTCACCAACCATCCATAAATCGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                      <  1:305483/150‑1 (MQ=255)
                          ccATCGAGAGCCCCAGTGCACCGATAATCGTCCCCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                                                                                                     <  1:1610863/150‑1 (MQ=255)
                          ccATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                                                                                                     >  1:1808867/1‑150 (MQ=255)
                          ccATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                                                                                                     >  2:452293/1‑150 (MQ=255)
                             tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgt                                                                                                                  <  2:1357967/150‑1 (MQ=255)
                              cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                 <  2:1786291/150‑1 (MQ=255)
                               gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCTCCAGTAGCGGCCAGGTTTGCAGGTAACAGAACAACACCGGACCCCATAATATTCCCCGCCAACATAAGGTTGACGGGGATTAAGCCCACGTtgtgc                                                                                                                >  2:121671/1‑149 (MQ=255)
                               gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGa                                                                                                                >  1:432413/1‑150 (MQ=255)
                               gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGa                                                                                                                <  2:822660/150‑1 (MQ=255)
                               gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGa                                                                                                                <  1:1641935/150‑1 (MQ=255)
                                agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                               <  2:1375679/150‑1 (MQ=255)
                                    ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                           <  1:1365093/150‑1 (MQ=255)
                                    ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                           >  2:1568372/1‑150 (MQ=255)
                                     cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatca                                                                                                          >  2:1376047/1‑150 (MQ=255)
                                      ccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCATGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACAACATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcag                                                                                                         <  2:944205/150‑1 (MQ=255)
                                      ccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcag                                                                                                         <  2:836411/150‑1 (MQ=255)
                                      ccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcag                                                                                                         <  1:1478211/150‑1 (MQ=255)
                                      ccAGTGCACCGATAATCGGCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcag                                                                                                         <  1:1634672/150‑1 (MQ=255)
                                        aGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagca                                                                                                       <  2:945410/150‑1 (MQ=255)
                                            cACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                   >  1:685002/1‑150 (MQ=255)
                                               cGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                <  2:1375400/150‑1 (MQ=255)
                                                 aTAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGAc                                                                                              <  2:642825/150‑1 (MQ=255)
                                                      cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCCCTGTGTGAGCATCAGCATCCGAAGACATAAt                                                                                         >  1:1216020/1‑150 (MQ=255)
                                                      cGTCACCAACCATCCATAAATGGCAATCCCGCAAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACTCCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                                                                                         >  2:375484/1‑150 (MQ=255)
                                                       gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAACACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACTGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                        >  2:1453187/1‑150 (MQ=255)
                                                       gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                        >  1:1522620/1‑150 (MQ=255)
                                                       gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                        >  1:1515517/1‑150 (MQ=255)
                                                        tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                       >  1:106246/1‑150 (MQ=255)
                                                        tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                       <  2:685851/150‑1 (MQ=255)
                                                           ccaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAACCACCTGCCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCAGCCGAAGACATAATAAAAc                                                                                    >  1:1165810/1‑150 (MQ=255)
                                                           ccaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAc                                                                                    <  1:1089902/150‑1 (MQ=255)
                                                             aaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                  <  2:1768904/150‑1 (MQ=255)
                                                                   ccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCTACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                            <  1:407745/150‑1 (MQ=255)
                                                                     atAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAg                                                                          <  2:1037660/150‑1 (MQ=255)
                                                                      tAAATTGCACTCCCGCCAGTAGAGGCCAGTTTTGCAGGTAACAGAAAAACACCTGACCCCATAATAATCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                         <  1:344525/150‑1 (MQ=255)
                                                                      tAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                         >  1:1166082/1‑150 (MQ=255)
                                                                        aaTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCg                                                                       <  1:617673/150‑1 (MQ=255)
                                                                        aaTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGAATCCGAAGACATAATTAAACTCCTGCGAAGTCg                                                                       <  2:931134/150‑1 (MQ=255)
                                                                         aTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                      <  1:585485/150‑1 (MQ=255)
                                                                         aTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                      >  1:1611813/1‑150 (MQ=255)
                                                                              aaTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                 >  1:411982/1‑150 (MQ=255)
                                                                               aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                >  1:1703514/1‑150 (MQ=255)
                                                                                tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGGATCCGAAGACATAATTAACCTCATGCGAAGGCGAGcggcgt                                                               >  2:297703/1‑145 (MQ=255)
                                                                                   cGCCAGTAGAGGCCAGGTTTGCAGGTAACACAAAATCCCCTGACCCCATAATATTCCCCGACACCATCATGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                            >  2:757428/1‑150 (MQ=255)
                                                                                    gCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                           <  1:790884/150‑1 (MQ=255)
                                                                                    gCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                           <  1:550157/150‑1 (MQ=255)
                                                                                       aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                        >  2:249220/1‑150 (MQ=255)
                                                                                           gagGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAACCTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                    >  1:1154685/1‑150 (MQ=255)
                                                                                           gagGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                    >  1:16199/1‑150 (MQ=255)
                                                                                             ggCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                  >  2:1152724/1‑150 (MQ=255)
                                                                                                cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                               <  1:1634985/150‑1 (MQ=255)
                                                                                                  ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                             <  1:547023/150‑1 (MQ=255)
                                                                                                      tGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                         <  2:1070771/150‑1 (MQ=255)
                                                                                                         aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                      >  2:1538057/1‑150 (MQ=255)
                                                                                                            tAACAGAAAAACACCTGCCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                   >  2:353575/1‑150 (MQ=255)
                                                                                                               cAGAAAATCACCTGACCCCATTATATTCACCGTCACAATCAGTTTGACGGGGATTAAGCCCACTTTTTTAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTga                                <  1:1674300/150‑2 (MQ=255)
                                                                                                                 gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                              >  2:189356/1‑150 (MQ=255)
                                                                                                                 gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                              >  1:1766710/1‑150 (MQ=255)
                                                                                                                   aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                            <  2:361589/150‑1 (MQ=255)
                                                                                                                    aaaCACCTGACCCCATAATATTCCCCTACACCATCTGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACACTAAATTACGTCATATCATACGCCTGCATTCg                           <  2:1110407/150‑1 (MQ=255)
                                                                                                                    aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCCTCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCCTGACAATAAATTACGTCATATCATACGCCTGCAttag                           >  1:889645/1‑148 (MQ=255)
                                                                                                                    aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                           >  2:1451803/1‑150 (MQ=255)
                                                                                                                      acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                         <  1:841006/150‑1 (MQ=255)
                                                                                                                      acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                         <  2:120121/150‑1 (MQ=255)
                                                                                                                       cacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                        >  2:612883/1‑150 (MQ=255)
                                                                                                                        acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                       >  1:5157/1‑150 (MQ=255)
                                                                                                                         ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                      >  1:1300401/1‑150 (MQ=255)
                                                                                                                          cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGtt                     >  1:1645123/1‑150 (MQ=255)
                                                                                                                           tGACCCCATAATATTCCCAGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACCTAATTAAACTCCGGCGAAGGCGAGCTTCGTGACAATAAATTACGGCATATCATACGCCTGCATTCGATAGTTg                    <  1:121671/150‑1 (MQ=255)
                                                                                                                               cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                >  1:1279013/1‑150 (MQ=255)
                                                                                                                               cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                <  2:106246/150‑1 (MQ=255)
                                                                                                                               cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                <  2:1808867/150‑1 (MQ=255)
                                                                                                                               cccAAAATATTCCCCGACACCATCAGGGTGACTGGGATTAATCCCACATTGTGAGCCTCAGAATCCTAAGACATAAATAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                <  2:294218/150‑1 (MQ=255)
                                                                                                                                ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatta               >  2:1668045/1‑150 (MQ=255)
                                                                                                                                ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCGTCGTGACAAGAAATTACGTCATATCATACGCCTGCATTCGAGAGTTGAagta               >  1:278388/1‑150 (MQ=255)
                                                                                                                                 cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGTCATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat              <  2:390785/150‑1 (MQ=255)
                                                                                                                                  ggcttaatCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAGAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt             >  1:291835/8‑150 (MQ=255)
                                                                                                                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta            >  2:804709/1‑150 (MQ=255)
                                                                                                                                      tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAGTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACATTAAATTACGTCATAACATACTCATGCATTCGATAGTTGAATTATTATca         >  1:1203918/1‑149 (MQ=255)
                                                                                                                                          ccccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg     >  1:624058/1‑150 (MQ=255)
                                                                                                                                             ctacacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCTACGACATAATTAAACTCCTGCGAAGGCGAGCTTCGGGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTTAATTATAATCGACCGGGt  <  1:1603966/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGT  >  NC_000913/4336831‑4337121

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: ‑