Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,921,515 C→T T179T (ACG→ACA atpB ← F0 sector of membrane‑bound ATP synthase, subunit a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,921,5150CT98.7% 237.9 / ‑4.8 76T179T (ACG→ACAatpBF0 sector of membrane‑bound ATP synthase, subunit a
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (37/38);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AATGAAAATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921383‑3921647
                                                                                                                                    |                                                                                                                                    
aaTGAAAATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 >  2:1232943/1‑138 (MQ=255)
aaTGAAAATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 <  1:1232943/138‑1 (MQ=255)
            cTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAg                                                                                                                    <  1:994539/139‑1 (MQ=255)
                   gCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGa                                                                                                             <  1:789911/139‑1 (MQ=255)
                           gTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTc                                                                                                     >  2:826159/1‑139 (MQ=255)
                           gTTACCGAACAGTCACAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGCATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGACGCCGATGCCTTTc                                                                                                     >  2:483699/1‑139 (MQ=255)
                             tACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      >  2:1264281/1‑136 (MQ=255)
                             tACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      <  1:1264281/136‑1 (MQ=255)
                               ccGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACGGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCAttt                                                                                                 <  1:826159/139‑1 (MQ=255)
                                  aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  1:342142/139‑1 (MQ=255)
                                  aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  1:1141852/139‑1 (MQ=255)
                                  aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAAAGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  2:733994/139‑1 (MQ=255)
                                       tCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTCGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         >  1:119775/1‑139 (MQ=255)
                                                 aGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         >  1:654884/1‑129 (MQ=255)
                                                 aGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         <  2:654884/129‑1 (MQ=255)
                                                 aGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               >  1:734152/1‑139 (MQ=255)
                                                  gTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGTTGTAGAACAGAAt                                                                              >  1:387208/1‑139 (MQ=255)
                                                  gTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAAt                                                                              <  2:207758/139‑1 (MQ=255)
                                                     aaaCTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               <  2:669112/135‑1 (MQ=255)
                                                     aaaCTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               >  1:669112/1‑135 (MQ=255)
                                                           gTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATaaa                                                                     >  1:1091980/1‑139 (MQ=255)
                                                           gTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATaaa                                                                     >  1:579221/1‑139 (MQ=255)
                                                                 aCAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               <  2:43313/139‑1 (MQ=255)
                                                                   agcagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:1091980/139‑1 (MQ=255)
                                                                   agcagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:242892/139‑1 (MQ=255)
                                                                   agcagGCTTACCCCTTCAAGGATTAAGTCGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:119775/139‑1 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGcc                                                                                                         <  2:1163953/93‑1 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGcc                                                                                                         >  1:1163953/1‑93 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCa                                                                                                    >  1:897103/1‑98 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               <  1:432106/119‑1 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  1:502640/1‑137 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:502640/137‑1 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGGCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               >  2:432106/1‑119 (MQ=255)
                                                                     cagGCTTACCCCTTCAAGGATTAAGTTGACAGCAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCa                                                                                                    <  2:897103/98‑1 (MQ=255)
                                                                             ccccTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 >  1:1008192/1‑61 (MQ=255)
                                                                             ccccTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 <  2:1008192/61‑1 (MQ=255)
                                                                                       gATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACg                                         >  1:364421/1‑139 (MQ=255)
                                                                                        aTTAAGTTGACAGGAATGAACGCTCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:327313/118‑1 (MQ=255)
                                                                                        aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:541070/118‑1 (MQ=255)
                                                                                        aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  1:327313/1‑118 (MQ=255)
                                                                                        aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  2:541070/1‑118 (MQ=255)
                                                                                           aaGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCt                                                                                                        <  2:601013/72‑1 (MQ=255)
                                                                                           aaGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCt                                                                                                        >  1:601013/1‑72 (MQ=255)
                                                                                            aGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCAc                                    >  1:523133/1‑139 (MQ=255)
                                                                                                 aCAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCg                               <  2:523133/139‑1 (MQ=255)
                                                                                                  cAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTGAACTCTTTCGTGAAGCCGCCGATGTCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGACATAGACAGCGTTACGTTCACGTCCGc                              >  2:735301/1‑139 (MQ=255)
                                                                                                        tGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGaa                                                                                    >  1:216747/1‑79 (MQ=255)
                                                                                                        tGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGaa                                                                                    <  2:216747/79‑1 (MQ=255)
                                                                                                          aaCGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGaa                      >  1:213622/1‑139 (MQ=255)
                                                                                                           aCGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGcc                                                        <  1:1009521/104‑1 (MQ=255)
                                                                                                           aCGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGcc                                                        >  2:1009521/1‑104 (MQ=255)
                                                                                                               ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACcaca                 >  1:463154/1‑139 (MQ=255)
                                                                                                                 aGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAAt                                                                              <  2:854437/76‑1 (MQ=255)
                                                                                                                 aGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAAt                                                                              >  1:854437/1‑76 (MQ=255)
                                                                                                                 aGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACg               <  2:213622/139‑1 (MQ=255)
                                                                                                                    gATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCg                               <  2:1113345/120‑1 (MQ=255)
                                                                                                                    gATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCg                               >  1:1113345/1‑120 (MQ=255)
                                                                                                                      ttGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGTTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTg          <  2:387208/139‑1 (MQ=255)
                                                                                                                          aCGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcagg      <  2:579221/139‑1 (MQ=255)
                                                                                                                          aCGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcagg      <  2:364421/139‑1 (MQ=255)
                                                                                                                           cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  1:684897/1‑83 (MQ=255)
                                                                                                                           cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:684897/83‑1 (MQ=255)
                                                                                                                           cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggc     <  2:734152/139‑1 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGa                                                                                <  2:1095682/62‑1 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGa                                                                                >  1:1095682/1‑62 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGTCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggcag   >  2:23426/1‑139 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAAcc                    >  2:560101/1‑122 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAAcc                    <  1:560101/122‑1 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggcag   >  1:1030406/1‑139 (MQ=255)
                                                                                                                             gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggcag   >  2:286207/1‑139 (MQ=255)
                                                                                                                              cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  <  2:122155/139‑1 (MQ=255)
                                                                                                                              cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  >  1:1252630/1‑139 (MQ=255)
                                                                                                                              cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  >  1:1106876/1‑139 (MQ=255)
                                                                                                                              cTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATATTGATGCTGTAGAATAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  <  2:213274/139‑1 (MQ=255)
                                                                                                                                gCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            <  1:95408/47‑1 (MQ=255)
                                                                                                                                gCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            >  2:95408/1‑47 (MQ=255)
                                                                                                                                    |                                                                                                                                    
AATGAAAATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921383‑3921647

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