Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,427,622 T→G 6.9% N57T (AAC→ACC)  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,6220TG 6.9% 249.6 / 10.9 73N57T (AAC→ACC) argTlysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base G (2/3);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GGAGAAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGA  >  NC_000913/2427474‑2427767
                                                                                                                                                    |                                                                                                                                                 
ggAGAAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGt                                                                                                                                                   >  1:104770/1‑149 (MQ=255)
   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                                                                                                                                          >  2:146790/1‑149 (MQ=255)
         aaTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGAt                                                                                                                                          >  4:205000/1‑149 (MQ=255)
          aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                         <  2:398595/149‑1 (MQ=255)
          aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                         <  1:451074/149‑1 (MQ=255)
          aTCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATc                                                                                                                                         <  4:286211/149‑1 (MQ=255)
           tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                        <  2:41308/149‑3 (MQ=255)
           tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                        <  2:38453/149‑3 (MQ=255)
           tctcCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGGTAACGAGAtgt                                                                                                                                        <  1:226826/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                                                                                                 <  1:268533/149‑1 (MQ=255)
                                                  taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                 <  2:335370/149‑1 (MQ=255)
                                                  taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                 <  4:201605/149‑1 (MQ=255)
                                                  taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                 <  4:490560/149‑1 (MQ=255)
                                                  taataGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGAt                                                                                                 <  4:475219/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                                    >  3:167339/1‑149 (MQ=255)
                                                                                                               ccAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGc                                    >  1:490074/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:420013/149‑1 (MQ=255)
                                                                                                                                               tctcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACaaa    <  2:407826/149‑1 (MQ=255)
                                                                                                                                                 tcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGa  >  2:44197/1‑149 (MQ=255)
                                                                                                                                                 tcGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGa  <  1:451570/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGAGAAGGCAATCTCCTGCTGACGTTTATCGGTAATGGAAAGCGACGAAATAATAGCGTCGATTTTTTTCGCTTTCAGTGAGGGGATCAGCGCGTCAAAGTCACTGGCAACCCAGGTACATTTCACCTGCATCCGTTTGCACATCTCGTTACCGAGATCGATATCAAAGCCAACAAAATCACCTTTAGCATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGA  >  NC_000913/2427474‑2427767

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