Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,015 A→C 7.0% D304A (GAT→GCT)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,0150AC 7.0% 213.9 / 6.4 71D304A (GAT→GCT) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base C (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CATCAAGGCCATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAC  >  NC_002947/5961869‑5962157
                                                                                                                                                    |                                                                                                                                                
cATCAAGGCCATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATg                                                                                                                                               >  5:148283/1‑149 (MQ=255)
 aTCAAGGCCATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATg                                                                                                                                               >  2:232088/1‑148 (MQ=255)
         cATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGa                                                                                                                                      <  8:22781/149‑1 (MQ=255)
          aTCCGCGAGCGCGAGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGACGCTACATCCTCGCCAACGCCGCCGACCTGCTGACCCGCCTGGAATACCTCAAGCACACCGAACACAAGGACTTCTCCATCATC‑‑GATCAAGCGATGaa                                                                                                                                     >  4:188066/1‑149 (MQ=255)
             cgcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGCCCCGCGTGGAATCCCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGa                                                                                                                                  >  7:338461/1‑149 (MQ=255)
             cgcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GCTGCATCGATGatgg                                                                                                                                   <  7:215714/148‑4 (MQ=255)
              gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGAc                                                                                                                                 >  2:165203/1‑149 (MQ=255)
                  gcgcGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGa                                                                                                                             >  1:239711/1‑149 (MQ=255)
                            gAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGaa                                                                                                                                     <  1:165706/131‑1 (MQ=255)
                            gAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGaa                                                                                                                                     >  2:165706/1‑131 (MQ=255)
                             aGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACCCAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGcc                                                                                                                 <  8:338461/149‑1 (MQ=255)
                                  gACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTa                                                                                                            >  6:137547/1‑149 (MQ=255)
                                          ccTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCt                                                                                                    <  6:148283/149‑1 (MQ=255)
                                          ccTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCt                                                                                                    <  8:152356/149‑1 (MQ=255)
                                          ccTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCCGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑ACTGCATCGATGatg                                                                                                                                    >  2:180947/1‑116 (MQ=255)
                                          ccTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCCGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑ACTGCATCGATGatg                                                                                                                                    <  1:180947/118‑3 (MQ=255)
                                            tGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTgg                                                                                                  >  5:212925/1‑149 (MQ=255)
                                                 ttCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATgg                                                                                              <  6:322717/148‑1 (MQ=255)
                                                   cGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGgt                                                                                            <  7:89974/148‑1 (MQ=255)
                                                    gAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGTCATGGGCATCATCCGGAGTGCAGCGATCAACGCCATGATCCGCCCGGCCCTGAACCAGGCCTGGATggt                                                                                             >  1:122477/1‑148 (MQ=255)
                                                      aCCAGGCTGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGCTCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        <  1:254868/149‑1 (MQ=255)
                                                       ccAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCACCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        <  1:232088/149‑1 (MQ=255)
                                                        cAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCa                                                                                                         >  3:245942/1‑130 (MQ=255)
                                                        cAGGCAGCTACATCGTCGCCAACGACGGCGTGCTTCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCTCCCGGCCCTGTACCa                                                                                                         <  4:245942/130‑1 (MQ=255)
                                                                  cATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACACGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGc                                                                            <  2:196272/149‑1 (MQ=255)
                                                                   aTCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGcc                                                                           >  8:305862/1‑149 (MQ=255)
                                                                         gCCAACGCCGGCGTGCTGCTGTCCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        <  2:299788/130‑1 (MQ=255)
                                                                         gCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        >  1:299788/1‑130 (MQ=255)
                                                                          ccAACGCCGGCGTGCCGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTg                                                                      >  6:328396/1‑147 (MQ=255)
                                                                             aCGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGc                                                                 >  8:284529/1‑149 (MQ=255)
                                                                                   gCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTc                                                                                         <  6:271579/119‑1 (MQ=255)
                                                                                   gCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTc                                                                                         >  5:271579/1‑119 (MQ=255)
                                                                                      tgctgcTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCcgc                                                        <  1:165203/149‑1 (MQ=255)
                                                                                         tgctgACCAGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACCTGCAGGCGAAGGCCGCGcc                                                     >  5:256815/1‑149 (MQ=255)
                                                                                          gctgACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGcc                                                     >  6:58671/1‑148 (MQ=255)
                                                                                                    gTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGt                                          >  6:200537/1‑149 (MQ=255)
                                                                                                         aTACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        >  3:264553/1‑98 (MQ=255)
                                                                                                         aTACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        <  4:264553/98‑1 (MQ=255)
                                                                                                            ccTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTa                                  >  4:77861/1‑149 (MQ=255)
                                                                                                             cTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTAt                                 <  2:239711/149‑1 (MQ=255)
                                                                                                             cTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTAt                                 >  3:163698/1‑149 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTc                                         <  5:272119/140‑1 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTc                                         >  6:272119/1‑140 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  5:218749/1‑149 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  8:31673/1‑149 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  2:78186/1‑149 (MQ=255)
                                                                                                              tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACATGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTAAGCGCGGTCAAGACACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTagc                                >  8:85066/1‑147 (MQ=255)
                                                                                                               cAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCt                               >  6:191887/1‑149 (MQ=255)
                                                                                                                aaGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCt                               <  5:58671/148‑1 (MQ=255)
                                                                                                                aaGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTg                              <  5:200537/149‑1 (MQ=255)
                                                                                                                 aGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGc                             <  2:122477/149‑1 (MQ=255)
                                                                                                                  gCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGc                             >  5:199256/1‑148 (MQ=255)
                                                                                                                     cacCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCAGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCgaga                         >  5:361058/1‑149 (MQ=255)
                                                                                                                       ccGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACATGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACt                       >  2:41645/1‑149 (MQ=255)
                                                                                                                         gAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCt                               <  6:53759/139‑1 (MQ=255)
                                                                                                                         gAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCt                               >  5:53759/1‑139 (MQ=255)
                                                                                                                          aaCACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGc                    <  3:188066/149‑1 (MQ=255)
                                                                                                                          aaCACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCCAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGc                    <  1:78186/149‑1 (MQ=255)
                                                                                                                           acacAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCg                   >  8:57947/1‑149 (MQ=255)
                                                                                                                             acaAGGACTTCGCCATCATC‑‑GATGCGGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTgg                                           <  5:89854/123‑1 (MQ=255)
                                                                                                                             acaAGGACTTCGCCATCATC‑‑GATGCGGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTgg                                           >  6:89854/1‑123 (MQ=255)
                                                                                                                             acaAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGAc                 <  5:137547/149‑1 (MQ=255)
                                                                                                                               aaGGACTTCGCCATCATC‑‑GATGCAGCGATGACCGCCCTGAGCCGCCCGGCCCTGTCCCAGGCCGGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCGATCTGCGAGACTGGCGagag               >  5:305058/1‑146 (MQ=255)
                                                                                                                               aaGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGc                    <  6:305058/144‑1 (MQ=255)
                                                                                                                                     ttCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCgg         >  4:190528/1‑149 (MQ=255)
                                                                                                                                       cGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCgg         <  5:328396/147‑1 (MQ=255)
                                                                                                                                         cgttcatc‑‑GCTGCATCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGc        >  8:148833/4‑146 (MQ=255)
                                                                                                                                         cgttcatc‑‑GCTGCATCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGc        <  7:148833/143‑1 (MQ=255)
                                                                                                                                         ccATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAgcgc                                                                                    >  5:328758/1‑70 (MQ=255)
                                                                                                                                         ccATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAgcgc                                                                                    <  6:328758/70‑1 (MQ=255)
                                                                                                                                            tcatc‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAc  <  7:305862/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CATCAAGGCCATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAC  >  NC_002947/5961869‑5962157

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