Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,307,634 T→A 100% I338F (ATC→TTC)  clsA ← cardiolipin synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,307,6340TA98.9% 292.2 / ‑4.3 88I338F (ATC→TTC) clsAcardiolipin synthase 1
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (41/46);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGTG  >  NC_000913/1307539‑1307730
                                                                                               |                                                                                                
cgGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACaaa                                                                                              >  2:378413/1‑100 (MQ=255)
   aaGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAatat                                                                                           <  1:2486998/100‑1 (MQ=255)
       ataataCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCAc                                                                                       >  2:1571938/1‑100 (MQ=255)
          ataCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCAcgcg                                                                                    <  1:2973069/100‑1 (MQ=255)
           taCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACgcgc                                                                                   >  1:2568340/1‑100 (MQ=255)
             cTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                 >  1:767072/1‑100 (MQ=255)
             cTGACATCCACCCCACGCAACAACGCCGTCCAAATAACATGAAATAAAACATCGCTCGTCACAAAACAGAGCTAGGACATGAACAAAAATTCAcgc                                                                                     >  2:804450/1‑96 (MQ=255)
              ttaCATCCTCCCCTCGCTGTGCCGCCTTGCAAATCGCATGTAGTAAATCATCCCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGa                                                                                <  1:2731684/98‑1 (MQ=255)
               gACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGaa                                                                               >  1:1272787/1‑100 (MQ=255)
               gACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGaa                                                                               >  1:2776137/1‑100 (MQ=255)
               gACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGaa                                                                               >  2:3017753/1‑100 (MQ=255)
                aCATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAAt                                                                              <  2:30200/100‑1 (MQ=255)
                 cATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTATGGCGTGGTCATGAACAAATATTCACGCGCCGAATa                                                                             <  1:1836833/100‑1 (MQ=255)
                 cATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATa                                                                             >  2:1688437/1‑100 (MQ=255)
                  aTCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATaa                                                                            >  2:2026414/1‑100 (MQ=255)
                   tCCACCCCGCGCTGCGCCGCCGTACAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATAAACAAATATTCACGCGCCGAATAAg                                                                           >  1:991763/1‑100 (MQ=255)
                    ccACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                          <  2:2124210/100‑1 (MQ=255)
                    ccACCCCGCGCTGCGCCGCCGTGCAAATCGCAGGAAGTACATTATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                          <  2:2129301/100‑1 (MQ=255)
                     cACCCCGCGCTGCGCCGCCGTGCAAATCGCATGCAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGGGCCGAATAAGc                                                                          <  2:1591686/99‑1 (MQ=255)
                     cACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                          <  1:2237285/99‑1 (MQ=255)
                     cACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGAGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGcc                                                                         >  1:569893/1‑100 (MQ=255)
                       ccccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGcc                                                                         <  2:1201742/98‑1 (MQ=255)
                         ccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTATGGCGTGGTCATGAACAGCTATTCACGCGCCGAATAAGCCgt                                                                       <  2:2168770/98‑2 (MQ=255)
                          cgcgCTGCTCCGCCGTGCAAATCGTATGAAGTAAATCATCGCTTTGCACAAAGTAGTGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                 <  2:1262858/87‑1 (MQ=255)
                            cgcTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCa                                                                  <  2:2498964/100‑1 (MQ=255)
                              cTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaat                                                                <  2:716796/100‑1 (MQ=255)
                              cTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaat                                                                <  1:1842182/100‑1 (MQ=255)
                                gcgcCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaataa                                                              >  2:948409/1‑100 (MQ=255)
                                gcgcCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGACATGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaataa                                                              <  2:444673/100‑1 (MQ=255)
                                  gccgccGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATg                                                            >  2:1721094/1‑100 (MQ=255)
                                     gccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCt                                                         <  1:172376/100‑1 (MQ=255)
                                     gccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCt                                                         <  1:923119/100‑1 (MQ=255)
                                       cgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTgg                                                       <  1:967340/100‑1 (MQ=255)
                                         tGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTg                                                     <  2:1896869/100‑1 (MQ=255)
                                             aaTCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATg                                                 <  1:1297098/100‑1 (MQ=255)
                                              aTCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATga                                                >  2:513091/1‑100 (MQ=255)
                                              aTCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATga                                                <  2:771080/100‑1 (MQ=255)
                                                 gCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGAt                                             <  2:1074895/100‑1 (MQ=255)
                                                 gCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGAt                                             >  1:3028492/1‑100 (MQ=255)
                                                 gCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGAt                                             >  2:2633991/1‑100 (MQ=255)
                                                   aTGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCt                                           >  1:530869/1‑100 (MQ=255)
                                                     gAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTc                                         >  1:1873178/1‑100 (MQ=255)
                                                     gAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTc                                         >  2:2284537/1‑100 (MQ=255)
                                                        gTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCgg                                      <  2:2509638/100‑1 (MQ=255)
                                                        gTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCgg                                      <  2:2863411/100‑1 (MQ=255)
                                                        gTAAATCATCACTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACTCGCCGAATAAGACGCAGTCAATAATGCCTGGTGAATGAGATCTTCCgg                                      >  1:407543/1‑100 (MQ=255)
                                                          aaaTCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGaa                                    <  1:874274/100‑1 (MQ=255)
                                                          aaaTCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGaa                                    >  2:981087/1‑100 (MQ=255)
                                                          aaaTCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGaa                                    >  1:2799024/1‑100 (MQ=255)
                                                              catcGCTTGGCACAAAGTAGGGCGTGGTCATTAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGcc                                <  1:1120952/100‑1 (MQ=255)
                                                               atcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCg                               >  1:939601/1‑100 (MQ=255)
                                                               atcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCg                               >  2:2235306/1‑100 (MQ=255)
                                                               atcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGATAAATGAGATCTTCCGGAAAGCCg                               >  1:1413534/1‑100 (MQ=255)
                                                               atcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCAAATAAACCGCAGACAATAATGACTGGTGAATGAGATCTTCCGGAAAGCCg                               >  1:121130/1‑100 (MQ=255)
                                                                 cGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCggg                             <  1:2266886/100‑1 (MQ=255)
                                                                  gCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATTAGATCTTCCGGAAAGCCggg                             >  1:723703/1‑99 (MQ=255)
                                                                  gCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCgggg                            >  2:520845/1‑100 (MQ=255)
                                                                  gCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCgggg                            <  1:659000/100‑1 (MQ=255)
                                                                  gCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCggcg                            >  2:1840607/1‑98 (MQ=255)
                                                                   cTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGc                           <  2:1883286/100‑1 (MQ=255)
                                                                    ttGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGcc                          <  1:2378480/100‑1 (MQ=255)
                                                                     tGGCACAAAGTAGGGCGTGGTCATGAATAAATATTCACGCGCTGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCa                         <  1:1761050/100‑1 (MQ=255)
                                                                     tGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCa                         <  2:324494/100‑1 (MQ=255)
                                                                       ggtgaaaGGAGGGCGGGGTCATGAACAAAGAGTCAGGCGCGGAATAAGTGGCAGGTCATAATGCCTGGTGAATGAGATCGTCCGGAAAGCCGGGGCCAg                        <  2:2313147/95‑1 (MQ=255)
                                                                       gCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGa                       >  2:1364715/1‑100 (MQ=255)
                                                                       gCAAAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAg                        <  2:1332895/99‑1 (MQ=255)
                                                                          caAAGTAGGGCGTTGTCATGATCAAATATTCACGTGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGc                    <  1:393534/100‑1 (MQ=255)
                                                                            aaGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTTGTGAATGAGATCTTCCGGAAAGCCTGGGCCAGAAGCaa                  <  2:1873156/100‑1 (MQ=255)
                                                                             aGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAAt                 <  2:2635940/100‑1 (MQ=255)
                                                                              gTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAAAGAAATCTTCCGGAAAGCCGGGGCCAGAAGCAAtt                >  1:992819/1‑100 (MQ=255)
                                                                               tAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtg               <  1:2942031/100‑1 (MQ=255)
                                                                                aGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgt              >  1:1384725/1‑100 (MQ=255)
                                                                                aGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgt              <  1:1401434/100‑1 (MQ=255)
                                                                                 ttccgatctCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtg             <  1:2213299/92‑1 (MQ=255)
                                                                                 gggCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtg             >  1:2842605/1‑100 (MQ=255)
                                                                                 gggCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtg             >  1:2388953/1‑100 (MQ=255)
                                                                                  ggCGTGGTCATGAACAAATCTTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtg               >  1:1065207/1‑97 (MQ=255)
                                                                                  ggCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgt            >  2:2987057/1‑100 (MQ=255)
                                                                                   gCGTGGTCATGAACATATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgtg           <  1:3000305/100‑1 (MQ=255)
                                                                                   gCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgtg           <  1:2026437/100‑1 (MQ=255)
                                                                                   gCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgtg           <  1:1688457/100‑1 (MQ=255)
                                                                                     gTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATTCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGaa         <  1:316841/100‑1 (MQ=255)
                                                                                          cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgagatcggaa    >  2:2213276/1‑91 (MQ=255)
                                                                                          cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg    >  1:491620/1‑100 (MQ=255)
                                                                                          cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg    <  2:1679403/100‑1 (MQ=255)
                                                                                           aTGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtgt   <  1:1133341/100‑1 (MQ=255)
                                                                                            tGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGtgtg  >  1:1977234/1‑100 (MQ=255)
                                                                                            tGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGtgtg  >  1:2713308/1‑100 (MQ=255)
                                                                                               |                                                                                                
CGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGTG  >  NC_000913/1307539‑1307730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: