Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,559,452 C→T 5.2% G153G (GGG→GGA)  ddpC ← D,D‑dipeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,559,4520CT 5.2% 245.8 / 2.2 76G153G (GGG→GGA) ddpCD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

TGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGC  >  NC_000913/1559304‑1559598
                                                                                                                                                    |                                                                                                                                                  
tGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGc                                                                                                                                                    <  2:1355243/149‑1 (MQ=255)
  aTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGccc                                                                                                                                                  >  1:471305/1‑149 (MQ=255)
  aTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGccc                                                                                                                                                  >  2:1523634/1‑149 (MQ=255)
  aTCAGATGCCAACGAGACGCGCAAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGccc                                                                                                                                                  >  1:721208/1‑149 (MQ=255)
   tCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCg                                                                                                                                                 >  2:599577/1‑149 (MQ=255)
   tCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACAAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCg                                                                                                                                                 >  1:990589/1‑149 (MQ=255)
      gATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGTCCgaga                                                                                                                                              >  1:105742/1‑149 (MQ=255)
      gATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGTCCgaga                                                                                                                                              >  1:880467/1‑149 (MQ=255)
          ccAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCgg                                                                                                                                          <  2:1329165/149‑1 (MQ=255)
          ccAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCgg                                                                                                                                          <  2:1000796/149‑1 (MQ=255)
            aaCGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCa                                                                                                                                        >  1:1539587/1‑149 (MQ=255)
                  ccgcgCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGt                                                                                                                                  <  2:1318948/148‑1 (MQ=255)
                       ccAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCAt                                                                                                                             >  1:128683/1‑149 (MQ=255)
                        cAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCAtt                                                                                                                            >  1:914898/1‑149 (MQ=255)
                          aaGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGt                                                                                                                          <  1:707447/149‑1 (MQ=255)
                          aaGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGt                                                                                                                          <  2:716222/149‑1 (MQ=255)
                          aaGGTGTTCGGCGCCTGAACATAGGTATACTGGCGTACAACTAATCCGTGCCCCCGCGCCAGGCGCACATAAAACGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGt                                                                                                                          <  1:385524/149‑1 (MQ=255)
                            ggTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCa                                                                                                                        <  1:787190/149‑1 (MQ=255)
                              ttttCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGTCCGAGAGCGGCAGCCAGTGCCATTGTCAGt                                                                                                                      <  2:105742/149‑1 (MQ=255)
                              ttttCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGTCCGAGAGCGGCAGCCAGTGCCATTGTCAGt                                                                                                                      <  2:880467/149‑1 (MQ=255)
                                   gccgccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCaa                                                                                                                 >  1:287327/1‑149 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  2:1248720/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  1:947766/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  2:287327/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  2:914898/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  1:1022645/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGaa                                                                                                             <  2:1539587/149‑1 (MQ=255)
                                       ccTGAACATAGGTATACTGGCGTACAACTAATGCTTGCACCCGCTCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGTCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAGCgca                                                                                                             >  2:998588/1‑147 (MQ=255)
                                               tAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGa                                                                                                     <  2:1104492/149‑1 (MQ=255)
                                                aGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGAc                                                                                                    <  1:391232/149‑1 (MQ=255)
                                                aGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGAc                                                                                                    <  1:1500327/149‑1 (MQ=255)
                                                 ggTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACa                                                                                                   >  1:518700/1‑149 (MQ=255)
                                                  gTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAg                                                                                                  >  1:555234/1‑149 (MQ=255)
                                                    ataCTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCa                                                                                                <  1:551564/149‑1 (MQ=255)
                                                    ataCTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCa                                                                                                <  1:1156953/149‑1 (MQ=255)
                                                    ataCTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCa                                                                                                <  2:868979/149‑1 (MQ=255)
                                                     taCTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCAt                                                                                               >  2:441457/1‑149 (MQ=255)
                                                       cTGGCGTACAACTAATGCTTGCACCCGCGCCAGGCTCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTGAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTCCCAACGAAGGAAGCGACAGCATaa                                                                                             >  2:769535/1‑149 (MQ=255)
                                                          gCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGt                                                                                          >  1:793329/1‑149 (MQ=255)
                                                             tACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCa                                                                                       >  1:1145913/1‑149 (MQ=255)
                                                              aCAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCat                                                                                      <  2:1185917/149‑1 (MQ=255)
                                                                aaCTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCatga                                                                                    >  2:760071/1‑149 (MQ=255)
                                                                aaCTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCatga                                                                                    >  2:1054096/1‑149 (MQ=255)
                                                                    aaTGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATgcg                                                                                >  2:1209077/1‑149 (MQ=255)
                                                                      tGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCTGCATGGCTTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCa                                                                              <  1:1209077/149‑1 (MQ=255)
                                                                      tGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCa                                                                              >  2:984809/1‑149 (MQ=255)
                                                                       gCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCaa                                                                                                                 <  2:656789/113‑1 (MQ=255)
                                                                       gCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCaa                                                                                                                 >  1:656789/1‑113 (MQ=255)
                                                                           gCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCatgatg                                                                         >  2:803981/1‑149 (MQ=255)
                                                                                        gcgcACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTgcgc                                                            >  1:440730/1‑149 (MQ=255)
                                                                                            acaTAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGc                                                                    <  1:1412406/137‑1 (MQ=255)
                                                                                            acaTAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGc                                                                    >  2:1412406/1‑137 (MQ=255)
                                                                                                 aaaGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCGGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCaa                                                   <  2:471305/149‑1 (MQ=255)
                                                                                                 aaaGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCaa                                                   <  2:555234/149‑1 (MQ=255)
                                                                                                  aaGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAg                                                  <  1:1523634/149‑1 (MQ=255)
                                                                                                  aaGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAg                                                  <  2:1291930/149‑1 (MQ=255)
                                                                                                  aaGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAg                                                  <  1:535457/149‑1 (MQ=255)
                                                                                                    gggAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGca                                                >  2:1370739/1‑149 (MQ=255)
                                                                                                       aaTTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCacac                                             <  1:1492413/149‑1 (MQ=255)
                                                                                                            gCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGAt                                        >  2:912586/1‑149 (MQ=255)
                                                                                                             cacaATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATa                                       >  2:937984/1‑149 (MQ=255)
                                                                                                             cacaATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATa                                       <  2:1145913/149‑1 (MQ=255)
                                                                                                                     caatagcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCg                               >  2:1034534/1‑149 (MQ=255)
                                                                                                                      aatagcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGa                              <  1:525994/149‑1 (MQ=255)
                                                                                                                       atagcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCCCCAAGCACACCGGATAGACATCCGAg                             >  1:1442446/1‑149 (MQ=255)
                                                                                                                       atagcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGGCCACCGTATAGAACTCCGAg                             >  2:143843/1‑149 (MQ=255)
                                                                                                                       atagcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAg                             >  1:594283/1‑149 (MQ=255)
                                                                                                                         agcCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTa                           <  2:518700/149‑1 (MQ=255)
                                                                                                                           ccAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGc                         >  2:211872/1‑149 (MQ=255)
                                                                                                                              gCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGaa                      >  1:403596/1‑149 (MQ=255)
                                                                                                                               cATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAAc                     >  1:904778/1‑149 (MQ=255)
                                                                                                                                  ggCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCaa                  <  2:793329/149‑1 (MQ=255)
                                                                                                                                     gTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCa               <  2:990589/149‑1 (MQ=255)
                                                                                                                                        aaacaaacTCTGCCCGATAGCGGCAGCCAGTGCCTTTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATAc            <  1:330009/149‑1 (MQ=255)
                                                                                                                                        aaacaaacTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATAc            >  1:921457/1‑149 (MQ=255)
                                                                                                                                         aacaaacTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATAcc           >  1:124231/1‑149 (MQ=255)
                                                                                                                                          acaaacTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATATTGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATAccc          >  1:971392/1‑149 (MQ=255)
                                                                                                                                              acTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCaa      <  1:441457/149‑1 (MQ=255)
                                                                                                                                                tCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATg    >  2:824882/1‑149 (MQ=255)
                                                                                                                                                  ggCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGAAATCGACAGCATAATGTCCATGATGAGAATAATGAAGGAGTCTTCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCAGACCCGCAAGGGc  >  2:1128002/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGCCAGTGCCATTGTCAGTACCAACGAAGGAATCGACAGCATAATGTCCATGATGCGCATGATGATGGCGTCTGCGCGTCCACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGC  >  NC_000913/1559304‑1559598

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