Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,427,625 C→A 6.9% G56V (GGT→GTT)  argT ← lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,427,6250CA 6.9% 248.1 / 8.1 72G56V (GGT→GTT) argTlysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (2/3);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

GAAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAG  >  NC_000913/2427477‑2427772
                                                                                                                                                    |                                                                                                                                                   
gaAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTAc                                                                                                                                                     >  2:517691/1‑149 (MQ=255)
gaAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTAc                                                                                                                                                     <  1:54148/149‑1 (MQ=255)
   ggCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACga                                                                                                                                                  >  1:41308/1‑149 (MQ=255)
    gCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCgag                                                                                                                                                 >  2:451074/1‑149 (MQ=255)
     cAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCgaga                                                                                                                                                <  2:140097/149‑1 (MQ=255)
      aaTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGAt                                                                                                                                               >  4:205000/1‑149 (MQ=255)
      aaTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGAt                                                                                                                                               >  2:146790/1‑149 (MQ=255)
       aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                              <  4:286211/149‑1 (MQ=255)
       aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                              <  1:451074/149‑1 (MQ=255)
       aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                              <  2:398595/149‑1 (MQ=255)
        tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                             <  1:226826/149‑3 (MQ=255)
        tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                             <  2:41308/149‑3 (MQ=255)
        tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                             <  2:38453/149‑3 (MQ=255)
               ctgACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCaa                                                                                                                                      >  1:180948/1‑149 (MQ=255)
                   cGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGcc                                                                                                                                  <  3:548332/149‑1 (MQ=255)
                    gTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCca                                                                                                                                 >  2:451570/1‑149 (MQ=255)
                        aTCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCcaacaa                                                                                                                             <  4:523075/149‑1 (MQ=255)
                               aTGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCAcc                                                                                                                      <  2:523084/149‑1 (MQ=255)
                                  gAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCttt                                                                                                                   <  4:509805/149‑1 (MQ=255)
                                   aaaGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTa                                                                                                                  <  2:68399/149‑1 (MQ=255)
                                   aaaGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTa                                                                                                                  <  2:139896/149‑1 (MQ=255)
                                   aaaGCGACGAAATAACAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACCTTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTa                                                                                                                  <  2:458843/149‑1 (MQ=255)
                                    aaGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAg                                                                                                                 <  3:205000/149‑1 (MQ=255)
                                        gacgaAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCAACTGCATCCGTTTGCACATCACGTTACCGAGATCGATATCAAAGCCAACAAAATCTCCTTTaccttc                                                                                                             >  2:291647/1‑144 (MQ=255)
                                              ataataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGa                                                                                                       <  3:577577/149‑1 (MQ=255)
                                               taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                      <  4:475219/149‑1 (MQ=255)
                                               taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                      <  1:268533/149‑1 (MQ=255)
                                               taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                      <  4:201605/149‑1 (MQ=255)
                                               taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                      <  2:335370/149‑1 (MQ=255)
                                               taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                      <  4:490560/149‑1 (MQ=255)
                                                 ataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATga                                                                                                    <  1:120448/149‑1 (MQ=255)
                                                    gCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGaa                                                                                                 >  3:98921/1‑149 (MQ=255)
                                                        cGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGt                                                                                             <  4:98921/149‑1 (MQ=255)
                                                            ttttttCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGt                                                                                         <  2:125992/149‑1 (MQ=255)
                                                                ttCGCTTTCAGTGAGGGGATCGGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAggt                                                                                     <  3:65295/149‑1 (MQ=255)
                                                                           tGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTATAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTc                                                                          <  4:273088/149‑1 (MQ=255)
                                                                             aGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCg                                                                        >  3:136400/1‑149 (MQ=255)
                                                                              ggggATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGa                                                                       >  1:81991/1‑149 (MQ=255)
                                                                                 gATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATAc                                                                    <  4:294156/149‑1 (MQ=255)
                                                                                 gATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATAc                                                                    >  3:243145/1‑149 (MQ=255)
                                                                                  aTCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACg                                                                   >  3:291090/1‑149 (MQ=255)
                                                                                  aTCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGAACCGATAcc                                                                   <  1:480617/149‑2 (MQ=255)
                                                                                   tCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGt                                                                  >  4:211906/1‑149 (MQ=255)
                                                                                   tCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGAACCGATAcca                                                                  <  3:288841/149‑3 (MQ=255)
                                                                                      gcgcGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTAcc                                                               >  2:246342/1‑149 (MQ=255)
                                                                                           tCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGtctc                                                          >  2:150856/1‑149 (MQ=255)
                                                                                             aaaGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCg                                                        >  2:565899/1‑149 (MQ=255)
                                                                                              aaGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCgg                                                       >  1:473608/1‑149 (MQ=255)
                                                                                                gTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTa                                                     >  3:62726/1‑149 (MQ=255)
                                                                                                          aCCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATa                                           >  3:181196/1‑149 (MQ=255)
                                                                                                            ccAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGc                                         >  1:490074/1‑149 (MQ=255)
                                                                                                            ccAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGc                                         >  3:167339/1‑149 (MQ=255)
                                                                                                              aGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTg                                       >  1:585744/1‑149 (MQ=255)
                                                                                                              aGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTg                                       >  1:260903/1‑149 (MQ=255)
                                                                                                              aGGGACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTg                                       >  1:511503/1‑149 (MQ=255)
                                                                                                                      ttCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGc                               <  1:371106/149‑1 (MQ=255)
                                                                                                                         aCCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCtgt                            <  4:516099/149‑1 (MQ=255)
                                                                                                                          ccTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCtgtg                           >  4:135264/1‑149 (MQ=255)
                                                                                                                            tGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGga                         <  2:180948/149‑1 (MQ=255)
                                                                                                                             gCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCtgt                            <  3:539238/145‑1 (MQ=255)
                                                                                                                             gCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCtgt                            >  4:539238/1‑145 (MQ=255)
                                                                                                                             gCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGgag                        <  2:104770/149‑1 (MQ=255)
                                                                                                                               aTCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGgagag                      <  3:269569/149‑1 (MQ=255)
                                                                                                                                tCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGgagaga                     >  1:266848/1‑149 (MQ=255)
                                                                                                                                   gTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCg                  >  2:347236/1‑149 (MQ=255)
                                                                                                                                    tttGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGa                 <  4:243145/149‑1 (MQ=255)
                                                                                                                                       atcgaTCTCGGTAACGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTa              >  3:385792/5‑149 (MQ=255)
                                                                                                                                        cacaTCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTaa             <  2:528025/149‑1 (MQ=255)
                                                                                                                                           aTCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaa          >  4:69914/1‑149 (MQ=255)
                                                                                                                                            tctcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa         <  2:407826/149‑1 (MQ=255)
                                                                                                                                            tctcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa         <  2:420013/149‑1 (MQ=255)
                                                                                                                                              tcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGa       <  1:451570/149‑1 (MQ=255)
                                                                                                                                              tcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGa       >  2:44197/1‑149 (MQ=255)
                                                                                                                                                   aCCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACagag  >  4:463828/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAG  >  NC_000913/2427477‑2427772

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