Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 95,070 G→A 15.1% V141I (GTA→ATA)  murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091395,0700GA15.1% 328.7 / 39.7 132V141I (GTA→ATA) murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (55/56);  new base A (11/9);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

Only 100 of 137 total aligned reads displayed.

TCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTC  >  NC_000913/94971‑95166
                                                                                                   |                                                                                                
tcctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGgtg                                                                                                  >  1:1935558/1‑100 (MQ=255)
tcctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGgtg                                                                                                  <  2:75128/100‑1 (MQ=255)
 cctccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGgtgt                                                                                                 >  2:1065141/1‑100 (MQ=255)
  gtccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGAAATCTCAACACCGACATCGgtgt                                                                                                 >  1:338110/2‑99 (MQ=255)
   tccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAc                                                                                               >  1:438558/1‑100 (MQ=255)
    ccGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTAcc                                                                                              <  2:1910058/100‑1 (MQ=255)
      ggCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCg                                                                                             >  1:1495638/1‑99 (MQ=255)
      ggCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGa                                                                                            <  2:23499/100‑1 (MQ=255)
       gCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGAt                                                                                           >  2:1672849/1‑100 (MQ=255)
            aCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTACACCGATGTCgt                                                                                      <  2:1886432/100‑2 (MQ=255)
            aCCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAATGCGGCAACACGCTTTATAGGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGc                                                                                      <  1:1410251/100‑1 (MQ=255)
             ccTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCt                                                                                     >  2:356133/1‑100 (MQ=255)
             ccTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCt                                                                                     >  2:1232374/1‑100 (MQ=255)
                ccGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGtt                                                                                  <  1:2270276/100‑1 (MQ=255)
                  gTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttga                                                                                <  2:2272698/100‑2 (MQ=25)
                    tAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgaga                                                                              <  1:2265481/100‑4 (MQ=21)
                    tAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgaga                                                                              <  1:282647/100‑4 (MQ=25)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCTACACCGACATCGGTGTACCGATGACGCTGTTGCGCt                                                                             <  1:347604/100‑1 (MQ=255)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgagat                                                                             <  2:564772/100‑5 (MQ=21)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgagat                                                                             <  2:1557019/100‑5 (MQ=21)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCt                                                                             >  1:2252520/1‑100 (MQ=255)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCt                                                                             >  1:1217972/1‑100 (MQ=255)
                     aaaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCt                                                                             <  2:843901/100‑1 (MQ=255)
                      aaGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCtt                                                                            >  2:495298/1‑100 (MQ=255)
                       agagATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTa                                                                           <  2:2363237/100‑1 (MQ=255)
                        gagaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTaa                                                                          >  1:1624604/1‑100 (MQ=255)
                          gaTGTCGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAGCGTCATCGGTACACCGATGTCGttgttgagattgc                                                                          <  1:1147743/98‑8 (MQ=18)
                          gaTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACg                                                                        >  2:1621741/1‑100 (MQ=255)
                           aTGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGc                                                                       >  2:2111785/1‑100 (MQ=255)
                            tGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGcc                                                                      >  1:201081/1‑100 (MQ=255)
                            tGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGcc                                                                      >  2:530405/1‑100 (MQ=255)
                              aCGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCgg                                                                    <  2:1409497/100‑1 (MQ=255)
                              aCGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCgg                                                                    <  1:2310811/100‑1 (MQ=255)
                               cggcggCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCgg                                                                    >  2:2100163/1‑99 (MQ=255)
                               cggcggCGAGTTTAAGCCAGTGTGGCAACACGCTTTATGCGGCAGGCACTCTCAGCAGCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCgn                                                                    <  1:313459/99‑2 (MQ=255)
                                ggcggcGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGa                                                                   <  1:592072/99‑1 (MQ=255)
                                  cggcgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATa                                                                >  1:149415/1‑100 (MQ=255)
                                     cgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                             >  2:761243/1‑100 (MQ=255)
                                     cgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                             >  1:1221097/1‑100 (MQ=255)
                                     cgATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGa                                                             <  2:696019/100‑1 (MQ=255)
                                      gATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAt                                                            >  1:653895/1‑100 (MQ=255)
                                       aTTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGCCATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGAtt                                                           >  1:611059/1‑100 (MQ=255)
                                        ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTa                                                          >  1:495372/1‑100 (MQ=255)
                                        ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTa                                                          <  1:718922/100‑1 (MQ=255)
                                        ttttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTa                                                          >  2:1037686/1‑100 (MQ=255)
                                        ttttAAGCCACTGCGGCAACACGCCTTATACGGCAGGCAATCTAAACAACGACATTGGTTTACCGATGACGCTGATGCGCTCCACGCCGCAATCCGAtt                                                           >  1:2375311/1‑99 (MQ=255)
                                          ttAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACg                                                        <  2:1881029/100‑1 (MQ=255)
                                            aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGc                                                       <  1:2077307/99‑1 (MQ=255)
                                            aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCa                                                      <  2:1327074/100‑1 (MQ=255)
                                            aaGCCAGTGCGGCAACACGCTTTATACGGCAGGCAAACTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACACCGGAATACGATTACg                                                        >  2:2026610/1‑98 (MQ=255)
                                             aGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGGCGCTGTTGCGCTTAACGCCGGAATACGATTACGCAg                                                     >  2:1820624/1‑100 (MQ=255)
                                             aGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAg                                                     <  2:623227/100‑1 (MQ=255)
                                             aGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAg                                                     >  2:2191228/1‑100 (MQ=255)
                                              gCCAGTGCGGCAACACGCTTTGTACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGt                                                    <  2:1842158/100‑1 (MQ=255)
                                              gCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGt                                                    <  2:1735788/100‑1 (MQ=255)
                                                   tGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                               >  2:1971103/1‑100 (MQ=255)
                                                   tGCGGCAACACGCTTTATACGGCAGGCAATCTAAACAACTACCTCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTg                                               <  1:1815215/100‑1 (MQ=255)
                                                     cGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTCACGCCGGAATACGATTACGCAGTTATTGaa                                             >  1:1959424/1‑100 (MQ=255)
                                                      ggCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAAc                                            >  1:2513335/1‑100 (MQ=255)
                                                      ggCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAAc                                            <  2:2179888/100‑1 (MQ=255)
                                                           cacGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGgc                                       >  2:2453273/1‑100 (MQ=255)
                                                           cacGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGgc                                       >  2:337706/1‑100 (MQ=255)
                                                               cTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGa                                   <  1:1425748/100‑1 (MQ=255)
                                                                 ttATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAc                                 <  2:2460253/100‑1 (MQ=255)
                                                                  tataCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAcc                                <  1:1145639/100‑1 (MQ=255)
                                                                  tataCGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAAcc                                <  1:1621402/100‑1 (MQ=255)
                                                                      cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATc                             <  2:349080/99‑1 (MQ=255)
                                                                      cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCa                            >  2:2187572/1‑100 (MQ=255)
                                                                      cGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCa                            <  1:1829147/100‑1 (MQ=255)
                                                                        gcaggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                          >  2:363111/1‑100 (MQ=255)
                                                                        gcaggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAgg                          >  1:2265917/1‑100 (MQ=255)
                                                                         caggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGATTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAggg                         <  2:1521693/100‑1 (MQ=255)
                                                                          aggcaATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGc                        <  2:457703/100‑1 (MQ=255)
                                                                             caATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaa                     >  2:2568242/1‑100 (MQ=255)
                                                                              aaTCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGaaa                    <  2:1706968/100‑1 (MQ=255)
                                                                               agcgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAAt                   >  2:2535905/5‑100 (MQ=21)
                                                                               agcgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAAt                   >  2:614412/5‑100 (MQ=21)
                                                                                tctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGACGTTATTGAATTTGGCGCGAACCATCAGGGCGaa                     >  1:179656/1‑97 (MQ=255)
                                                                                gcgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATa                  >  1:2122397/4‑100 (MQ=25)
                                                                                gcgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATa                  >  1:450575/4‑100 (MQ=25)
                                                                                 ctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                 >  1:197465/1‑100 (MQ=255)
                                                                                 ctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                 >  2:2485742/1‑100 (MQ=255)
                                                                                 ctcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGAGAAATAg                 >  2:667288/1‑100 (MQ=255)
                                                                                 cgcaacagCGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAg                 >  2:139689/3‑100 (MQ=25)
                                                                                  tcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGc                <  2:937362/100‑1 (MQ=255)
                                                                                  tcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGc                >  2:1198751/1‑100 (MQ=255)
                                                                                  tcAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGc                >  1:712261/1‑100 (MQ=255)
                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc               <  2:2027533/100‑1 (MQ=255)
                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc               <  1:2485723/100‑1 (MQ=255)
                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc               >  1:1088880/1‑100 (MQ=255)
                                                                                   caacaaCGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGc                >  1:401476/1‑99 (MQ=255)
                                                                                    aacaacGAGATCGGTGTACCGATGAGGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGcc               <  1:1966629/99‑1 (MQ=255)
                                                                                     acagcGTCATCGGTACACCGATGTCGTTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTg             >  2:1847422/1‑100 (MQ=37)
                                                                                     acaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTg             <  1:1785596/100‑1 (MQ=255)
                                                                                     acaacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTg             >  1:1306871/1‑100 (MQ=255)
                                                                                      caacGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCAGGAATACGATTACGCAGTTGTTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTgg            >  1:1410118/1‑100 (MQ=255)
                                                                                          gACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtg        <  1:1934749/100‑1 (MQ=255)
                                                                                            cATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACtgtg      >  1:994769/1‑100 (MQ=255)
                                                                                             aTCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGa     <  2:1165349/100‑1 (MQ=255)
                                                                                                ggTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTc  <  2:1545642/100‑1 (MQ=255)
                                                                                                   |                                                                                                
TCCTCCGGCAAAACCTCCGTTAAAGAGATGACGGCGGCGATTTTAAGCCAGTGCGGCAACACGCTTTATACGGCAGGCAATCTCAACAACGACATCGGTGTACCGATGACGCTGTTGCGCTTAACGCCGGAATACGATTACGCAGTTATTGAACTTGGCGCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTC  >  NC_000913/94971‑95166

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