Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,020 G→T 7.2% A306S (GCG→TCG)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,0200GT 7.2% 220.7 / 5.8 69A306S (GCG→TCG) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base T (2/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAC  >  NC_002947/5961878‑5962157
                                                                                                                                                 |                                                                                                                                          
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                                                                                                                                     >  2:165706/1‑131 (MQ=255)
                   gAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGaa                                                                                                                                     <  1:165706/131‑1 (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                                                                                        <  4:264553/98‑1 (MQ=255)
                                                                                                aTACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCa                                                                                        >  3:264553/1‑98 (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                                         >  6:272119/1‑140 (MQ=255)
                                                                                                     tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTc                                         <  5:272119/140‑1 (MQ=255)
                                                                                                     tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  8:31673/1‑149 (MQ=255)
                                                                                                     tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  2:78186/1‑149 (MQ=255)
                                                                                                     tCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATc                                >  5:218749/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                               >  5:53759/1‑139 (MQ=255)
                                                                                                                gAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCt                               <  6:53759/139‑1 (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                                                                                    <  6:328758/70‑1 (MQ=255)
                                                                                                                                ccATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAgcgc                                                                                    >  5:328758/1‑70 (MQ=255)
                                                                                                                                   tcatc‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAc  <  7:305862/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CATCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATC‑‑GATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGGTCAAGCCACGTGCAGGCGAAGGCCGCGCCTACGACCTGGTCGGCCCTATCTGCGAGACTGGCGACTTCCTCGGCAAGGAC  >  NC_002947/5961878‑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: