Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,480,831 T→C T111A (ACC→GCC)  kefG ← protein required for KefB activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,480,8310TC99.6% 950.5 / ‑4.8 230T111A (ACC→GCC) kefGprotein required for KefB activity
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (113/116);  minor base A (0/1);  total (113/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 232 total aligned reads displayed.

GCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGTAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGC  >  NC_000913/3480733‑3480928
                                                                                                  |                                                                                                 
gCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCaa                                                                                                 >  1:5105144/1‑101 (MQ=255)
gCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCaa                                                                                                 <  1:5104617/101‑1 (MQ=255)
  acaTGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATc                                                                                               >  2:5105147/1‑101 (MQ=255)
  acaTGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATc                                                                                               <  2:5104779/101‑1 (MQ=255)
   caTGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCa                                                                                              >  1:5105148/1‑101 (MQ=255)
    aTGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCAc                                                                                             <  1:5104567/101‑1 (MQ=255)
     tGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACg                                                                                            >  2:5105150/1‑101 (MQ=255)
        ccgccgCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTa                                                                                         >  2:5105156/1‑101 (MQ=255)
        ccgccgCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTa                                                                                         >  2:5105154/1‑101 (MQ=255)
         cgccgcCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTAc                                                                                        <  1:5104687/101‑1 (MQ=255)
         cgccgcCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTAc                                                                                        >  1:5105157/1‑101 (MQ=255)
            cgccAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGcc                                                                                     <  1:5104743/101‑1 (MQ=255)
             gccAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCa                                                                                    <  2:5104667/101‑1 (MQ=255)
             gccAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCa                                                                                    >  1:5105160/1‑101 (MQ=255)
             gccAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCa                                                                                    <  1:5104883/101‑1 (MQ=255)
               cAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGt                                                                                  <  2:5104835/101‑1 (MQ=255)
               cAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGt                                                                                  >  2:5105163/1‑101 (MQ=255)
                aGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTa                                                                                 >  1:5105164/1‑101 (MQ=255)
                 gTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGcc                                                                                     >  2:5105171/1‑96 (MQ=255)
                 gTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTAc                                                                                >  1:5105168/1‑101 (MQ=255)
                 gTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTAc                                                                                >  1:5105170/1‑101 (MQ=255)
                 gTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTAc                                                                                <  1:5104820/101‑1 (MQ=255)
                   tCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACtt                                                                              <  1:5104465/101‑1 (MQ=255)
                    cAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACttt                                                                             >  1:5105173/1‑101 (MQ=255)
                     aaaGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTc                                                                            >  2:5105175/1‑101 (MQ=255)
                     aaaGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTc                                                                            <  2:5104710/101‑1 (MQ=255)
                      aaGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTcc                                                                           <  1:5104787/101‑1 (MQ=255)
                      aaGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTcc                                                                           >  1:5105178/1‑101 (MQ=255)
                           gcgcAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCa                                                                      <  2:5104698/101‑1 (MQ=255)
                           gcgcAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCa                                                                      <  2:5104586/101‑1 (MQ=255)
                            cgcAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAg                                                                     >  2:5105180/1‑101 (MQ=255)
                             gcagcaCATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTa                                                                                         >  2:5105185/1‑80 (MQ=255)
                             gcagcaCATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGt                                                                    >  2:5105183/1‑101 (MQ=255)
                              cagcaCATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGtt                                                                   >  1:5105186/1‑101 (MQ=255)
                               agcaCATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTg                                                                  <  1:5104769/101‑1 (MQ=255)
                                  acaTCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGtt                                                               <  2:5104829/101‑1 (MQ=255)
                                    aTCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTtc                                                             <  1:5104840/101‑1 (MQ=255)
                                    aTCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTtc                                                             >  1:5105189/1‑101 (MQ=255)
                                    aTCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTtc                                                             <  1:5104427/101‑1 (MQ=255)
                                     tCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTtcc                                                            <  1:5104862/101‑1 (MQ=255)
                                        cTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTtcctcc                                                         >  1:5105193/1‑101 (MQ=255)
                                          cATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCg                                                       >  2:5105195/1‑101 (MQ=255)
                                          cATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCg                                                       <  1:5104916/101‑1 (MQ=255)
                                           aTCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCgg                                                      >  1:5105197/1‑101 (MQ=255)
                                            tCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGc                                                     >  2:5105199/1‑101 (MQ=255)
                                             cGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGcc                                                    >  1:5105201/1‑101 (MQ=255)
                                              gggTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGccc                                                   <  1:5104633/101‑1 (MQ=255)
                                              gggTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGccc                                                   >  2:5105203/1‑101 (MQ=255)
                                              gggTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGc                                                     <  2:5104867/99‑1 (MQ=255)
                                               ggTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCg                                                  <  2:5104879/101‑1 (MQ=255)
                                                 tAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCt                                                <  1:5104802/101‑1 (MQ=255)
                                                 tAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCt                                                <  1:5104859/101‑1 (MQ=255)
                                                  aGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTg                                               <  2:5104813/101‑1 (MQ=255)
                                                   gCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTgg                                              <  1:5104756/101‑1 (MQ=255)
                                                    cGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGc                                             >  1:5105206/1‑101 (MQ=255)
                                                     gATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCa                                            <  1:5104966/101‑1 (MQ=255)
                                                     gATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCa                                            <  2:5104804/101‑1 (MQ=255)
                                                      aTTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCaa                                           >  2:5105209/1‑101 (MQ=255)
                                                       ttCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCaaa                                          >  1:5105211/1‑101 (MQ=255)
                                                        tCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCaaaa                                         <  2:5105024/101‑1 (MQ=255)
                                                         cAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAAc                                        >  2:5105213/1‑101 (MQ=255)
                                                          aGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAAcc                                       >  1:5105214/1‑101 (MQ=255)
                                                           gcgcGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCa                                      <  2:5104795/101‑1 (MQ=255)
                                                           gcgcGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCa                                      <  1:5104796/101‑1 (MQ=255)
                                                            cgcgTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCAc                                     >  1:5105218/1‑101 (MQ=255)
                                                            cgcgTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCAc                                     <  1:5105031/101‑1 (MQ=255)
                                                             gcgTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACg                                    >  2:5105220/1‑101 (MQ=255)
                                                             gcgTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACg                                    <  1:5104708/101‑1 (MQ=255)
                                                               gTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGAc                                  <  2:5104738/101‑1 (MQ=255)
                                                                tCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACt                                 >  2:5105226/1‑101 (MQ=255)
                                                                tCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACt                                 >  1:5105225/1‑101 (MQ=255)
                                                                tCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACt                                 >  2:5105223/1‑101 (MQ=255)
                                                                 cATAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACtt                                <  2:5105037/101‑1 (MQ=255)
                                                                  aTAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTa                               >  1:5105229/1‑101 (MQ=255)
                                                                  aTAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTa                               >  2:5105230/1‑101 (MQ=255)
                                                                   tAACGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTAc                                                                                        <  2:5104921/43‑1 (MQ=255)
                                                                    aaCGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACtt                                <  2:5104837/98‑1 (MQ=255)
                                                                    aaCGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAAt                             <  1:5104816/101‑1 (MQ=255)
                                                                     aCGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATa                            <  2:5104903/101‑1 (MQ=255)
                                                                     aCGGTAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATa                            >  2:5105234/1‑101 (MQ=255)
                                                                         tAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCg                        >  2:5105236/1‑101 (MQ=255)
                                                                         tAAGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCg                        <  2:5104905/101‑1 (MQ=255)
                                                                          aaGCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCgg                       >  2:5105239/1‑101 (MQ=255)
                                                                            gCACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTc                     >  1:5105240/1‑101 (MQ=255)
                                                                             cACTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAAt                             <  2:5104877/92‑1 (MQ=255)
                                                                              aCTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAAc                                        <  2:5104746/80‑1 (MQ=255)
                                                                               cTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGc                                             >  1:5105245/1‑74 (MQ=255)
                                                                               cTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTaa                              <  1:5104913/89‑1 (MQ=255)
                                                                               cTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTcc                    <  2:5104929/99‑1 (MQ=255)
                                                                               cTTTCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAg                  >  2:5105244/1‑101 (MQ=255)
                                                                                tttCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGc                 <  2:5104552/101‑1 (MQ=255)
                                                                                 ttCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGcc                >  1:5105248/1‑101 (MQ=255)
                                                                                 ttCCGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGcc                >  2:5105249/1‑101 (MQ=255)
                                                                                   ccGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCAc              <  1:5104898/101‑1 (MQ=255)
                                                                                    cGGCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAAt                             <  1:5104981/85‑1 (MQ=255)
                                                                                      gCTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCActct           >  2:5105253/1‑101 (MQ=255)
                                                                                       cTCGCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCtt          >  1:5105254/1‑101 (MQ=255)
                                                                                          gCCGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCa       <  1:5105007/101‑1 (MQ=255)
                                                                                           ccGGTGGCAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAg      <  1:5104962/101‑1 (MQ=255)
                                                                                              gtggcAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAg   >  1:5105258/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACTTTCCGGCTCGCCGGTGGTAATCACGCTACGCCAGTACTTTCCCGCCAGTTGGTTTCCTCCCGGCCCGCTGGCAAAACCACGACTTAATACCCGGTCCAGCCACTCTTTCAGTAGC  >  NC_000913/3480733‑3480928

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