Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,113,849 G→A 5.9% S308N (AGC→AAC)  clsA → cardiolipin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,113,8490GA 5.9% 209.2 / 2.3 69S308N (AGC→AAC) clsAcardiolipin synthase
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base A (2/2);  total (37/32)
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.71e-01

ACACCCACGTACAGATCAGCGGGCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCGGTGCTACGAGGAGT  >  NC_002947/6113702‑6113996
                                                                                                                                                   |                                                                                                                                                   
acacCCACGTACAGATCAGCGGGCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCg                                                                                                                                                   <  1:349834/150‑1 (MQ=255)
    ccACGTACAGATCAGCGGGCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGccc                                                                                                                                               <  1:216044/150‑1 (MQ=255)
              aTCAGCGGGCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCca                                                                                                                                     >  1:999924/1‑150 (MQ=255)
                      gCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAAc                                                                                                                             >  1:697190/1‑150 (MQ=255)
                       ccAGTACTGGCCTACCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAAcc                                                                                                                            >  2:355109/1‑150 (MQ=255)
                        cAGTACTGCCCTGCCTGCAAGATTCGTTCGCCGAGGACTGTTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCt                                                                                                                           <  2:560616/150‑1 (MQ=255)
                              tGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGc                                                                                                                     >  1:1123314/1‑150 (MQ=255)
                              tGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAACGGCCCGGCAGACCCACAAGAAACCTGCTCGc                                                                                                                     >  2:1067533/1‑150 (MQ=255)
                              tGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAACGGCCCGGCAGACCCACAAGAAACCTGCTCGc                                                                                                                     >  2:1067566/1‑150 (MQ=255)
                              tGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGACCTGTGTCAGGCCTTGGCCAGCg                                                                                                                                                   >  1:963276/1‑120 (MQ=255)
                              tGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGACCTGTGTCAGGCCTTGGCCAGCg                                                                                                                                                   <  2:963276/120‑1 (MQ=255)
                                gcctgcctGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCca                                                                                                                                     >  1:55748/1‑132 (MQ=255)
                                gcctgcctGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCca                                                                                                                                     <  2:55748/132‑1 (MQ=255)
                                gcctgcctGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCAAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTg                                                                                                                   >  2:735393/1‑150 (MQ=255)
                                  ctgcctgcAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGtt                                                                                                                 <  1:949561/150‑1 (MQ=255)
                                   tgcctgcAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACGTACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGttc                                                                                                                >  2:232715/1‑150 (MQ=255)
                                   tgcctgcAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGttc                                                                                                                <  2:697190/150‑1 (MQ=255)
                                        gcAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTACTCGCTGTTCTTCCt                                                                                                           >  2:975828/1‑150 (MQ=255)
                                          aaGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTg                                                                                                         >  2:761492/1‑150 (MQ=255)
                                            gagTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGaa                                                                                                       <  1:355109/150‑1 (MQ=255)
                                             agTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAg                                                                                                      <  1:538808/150‑1 (MQ=255)
                                                      ccGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGattcatt                                                                                             >  2:817937/1‑150 (MQ=255)
                                                         aGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCAAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCgg                                                                                          <  1:735393/150‑1 (MQ=255)
                                                              tGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCt                                                                                                                    >  1:1006952/1‑119 (MQ=255)
                                                              tGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCt                                                                                                                    <  2:1006952/119‑1 (MQ=255)
                                                              tGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCAcc                                                                                     >  2:33287/1‑150 (MQ=255)
                                                                    tGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAACGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAgg                                                                               <  1:1067566/150‑1 (MQ=255)
                                                                    tGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAACGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTCAAGCGATTCATTCGGCCACCCGCAgg                                                                               <  1:1067533/150‑1 (MQ=255)
                                                                      ggCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGt                                                                                                                  <  2:1195455/113‑1 (MQ=255)
                                                                      ggCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGt                                                                                                                  >  1:1195455/1‑113 (MQ=255)
                                                                      ggCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCt                                                                                                           >  2:417632/1‑120 (MQ=255)
                                                                      ggCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCt                                                                                                           <  1:417632/120‑1 (MQ=255)
                                                                         cACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGtgtg                                                                          >  1:300848/1‑150 (MQ=255)
                                                                                 aGCTACCGCCGCTGATCCTGCCCGATACGTACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAAc                                                                                                                             >  1:515479/1‑91 (MQ=255)
                                                                                 aGCTACCGCCGCTGATCCTGCCCGATACGTACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAAc                                                                                                                             <  2:515479/91‑1 (MQ=255)
                                                                                         ccgcTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGtgtg                                                                          >  1:1208260/1‑134 (MQ=255)
                                                                                         ccgcTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGtgtg                                                                          <  2:1208260/134‑1 (MQ=255)
                                                                                               aTCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATc                                                    <  1:33287/150‑1 (MQ=255)
                                                                                                 ccTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATccc                                                  >  2:217693/1‑150 (MQ=255)
                                                                                                    gCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACTCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGa                                               >  1:885455/1‑150 (MQ=255)
                                                                                                     cccGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGAc                                              >  2:226713/1‑150 (MQ=255)
                                                                                                         atacatacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGttcttc                                                                                                             <  2:237490/83‑1 (MQ=255)
                                                                                                         atacatacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGttcttc                                                                                                             >  1:237490/1‑83 (MQ=255)
                                                                                                          tacatacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCAc                                                                    >  1:387928/1‑123 (MQ=255)
                                                                                                          tacatacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCAc                                                                    <  2:387928/123‑1 (MQ=255)
                                                                                                           acatacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCACCGGCCCGGCAGACCCACAAGAAACCTGCTCGCGTTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGcc                                        <  1:681579/150‑1 (MQ=255)
                                                                                                             atacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGgtgt                                                                           >  1:685393/1‑113 (MQ=255)
                                                                                                             atacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGgtgt                                                                           <  2:685393/113‑1 (MQ=255)
                                                                                                              tacCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCgtg                                     >  2:1120354/1‑150 (MQ=255)
                                                                                                               acCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCgtgt                                    >  2:1113725/1‑150 (MQ=255)
                                                                                                               acCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCgtgt                                    <  1:1120354/150‑1 (MQ=255)
                                                                                                                 ccGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTc                                  <  1:366282/150‑1 (MQ=255)
                                                                                                                  cGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTcg                                 >  1:422326/1‑150 (MQ=255)
                                                                                                                   ggATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTcgc                                <  1:217693/150‑1 (MQ=255)
                                                                                                                    gATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTACTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTcgcc                               <  1:975828/150‑1 (MQ=255)
                                                                                                                                cTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTcc                                                                                                            >  1:831439/1‑61 (MQ=255)
                                                                                                                                cTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTcc                                                                                                            <  2:831439/61‑1 (MQ=255)
                                                                                                                                cTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTg                   >  1:199344/1‑150 (MQ=255)
                                                                                                                                 tgtgTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTgg                  >  1:335522/1‑150 (MQ=255)
                                                                                                                                    gtCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCgg               <  2:422326/150‑1 (MQ=255)
                                                                                                                                     tCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCgt                                      <  2:679234/126‑1 (MQ=255)
                                                                                                                                     tCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCgt                                      >  1:679234/1‑126 (MQ=255)
                                                                                                                                       aGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTgg                  <  1:622651/144‑1 (MQ=255)
                                                                                                                                       aGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTgg                  >  2:622651/1‑144 (MQ=255)
                                                                                                                                         gCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGc                        >  2:119323/1‑136 (MQ=255)
                                                                                                                                         gCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGc                        <  1:119323/136‑1 (MQ=255)
                                                                                                                                          ccTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCGGTGCTAc         >  1:768736/1‑150 (MQ=255)
                                                                                                                                           cTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCGGTGCTACg        >  1:816954/1‑150 (MQ=255)
                                                                                                                                                 cAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCGGTGCTACGAGGAGt  <  2:999924/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACACCCACGTACAGATCAGCGGGCCAGTACTGGCCTGCCTGCAAGAGTCGTTCGCCGAGGACTGGTACTGGGCCACGCGCCAGCTACCGCCGCTGATCCTGCCCGATACATACCCGGATAACGGCGTCCTGTGTCAGGCCTTGGCCAGCGGCCCGGCAGACCCACAAGAAACCTGCTCGCTGTTCTTCCTTGAAGCGATTCATTCGGCCACCCGCAGGGTGTGGATCACCAGCCCGTACTTCATCCCTGACGAGGCCGTGTTCGCCGCCTTGCGGCTGGCGGTGCTACGAGGAGT  >  NC_002947/6113702‑6113996

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