Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,019,626 A→T 5.6% V143D (GTC→GAC)  ompA ← outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,019,6260AT 5.6% 213.2 / 3.1 71V143D (GTC→GAC) ompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base T (2/2);  total (33/38)
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.85e-01

ACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTT  >  NC_000913/1019478‑1019774
                                                                                                                                                    |                                                                                                                                                    
aCCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGa                                                                                                                                                      >  1:146645/1‑149 (MQ=255)
    aaaCGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGGGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCgg                                                                                                                                                  <  7:170124/149‑1 (MQ=255)
     aaCGGTAGGAAACACCCAGGCGCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGGAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGga                                                                                                                                                 <  2:270908/149‑1 (MQ=255)
            ggAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGc                                                                                                                                          >  8:58138/1‑149 (MQ=255)
            ggAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGc                                                                                                                                          >  7:41045/1‑149 (MQ=255)
                 cacCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGgtg                                                                                                                                     <  3:244876/149‑1 (MQ=255)
                    ccAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCg                                                                                                                                  >  7:31433/1‑149 (MQ=255)
                     cAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGt                                                                                                                                 >  1:338986/1‑149 (MQ=255)
                            aGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGtttt                                                                                                                           >  6:271587/1‑148 (MQ=255)
                            aGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGttt                                                                                                                            <  2:154688/147‑1 (MQ=255)
                                tGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCa                                                                                                                      >  5:142510/1‑149 (MQ=255)
                                 gCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCa                                                                                                                      >  6:93473/1‑148 (MQ=255)
                                   cGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATa                                                                                                                    >  6:231793/1‑148 (MQ=255)
                                    gTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATaaa                                                                                                                  <  2:338986/149‑1 (MQ=255)
                                        tCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGtt                                                                                                              <  1:242695/149‑1 (MQ=255)
                                           ggACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGa                                                                                                           <  8:242285/149‑1 (MQ=255)
                                               gAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTa                                                                                                       <  2:270719/149‑1 (MQ=255)
                                               gAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTa                                                                                                       <  1:131954/149‑1 (MQ=255)
                                               gAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTa                                                                                                       <  8:31433/149‑1 (MQ=255)
                                                  tGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAgt                                                                                                     <  5:231793/148‑1 (MQ=255)
                                                    ccGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGGGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTc                                                                                                  >  5:232255/1‑149 (MQ=255)
                                                     cGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCt                                                                                                 <  7:58138/149‑1 (MQ=255)
                                                          gtgtgtGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCAc                                                                                             <  8:41045/148‑1 (MQ=255)
                                                             tgtgCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCa                                                                                         >  2:322037/1‑149 (MQ=255)
                                                               tgCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATa                                                                                       <  3:191784/149‑1 (MQ=255)
                                                                  gTCACCGATGTTTTTGGTCCACTGGTATTCAAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTAAACACCGCCAGCGAACACCGGAGAAACGCAGGTGTCGTGGTTTTTACCATAAACTATGGATTTAGTGTCTGCCCGCCATACCa                                                                                    >  8:112851/1‑149 (MQ=255)
                                                                   tCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCAt                                                                                   <  6:142510/149‑1 (MQ=255)
                                                                   tCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCAt                                                                                   <  3:151937/149‑1 (MQ=255)
                                                                     aCCGATGTTGTTGGTCCACTGCTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGc                                                                                 >  6:321069/1‑149 (MQ=255)
                                                                       cGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCa                                                                               >  8:639/1‑149 (MQ=255)
                                                                       cGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCa                                                                               >  3:184564/1‑149 (MQ=255)
                                                                            ttgttgGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCa                                                                           <  1:81021/148‑1 (MQ=255)
                                                                                 ggTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCATACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGa                                                                      <  3:216666/148‑1 (MQ=255)
                                                                                  gTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAgt                                                                    >  7:2784/1‑149 (MQ=255)
                                                                                  gTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAgt                                                                    <  5:93473/149‑1 (MQ=255)
                                                                                     cACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAgtg                                                                   >  4:275464/1‑147 (MQ=255)
                                                                                     cACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTa                                                                 <  2:146645/149‑1 (MQ=255)
                                                                                     cACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAAAGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTATTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTa                                                                 <  8:248782/149‑1 (MQ=255)
                                                                                      aCTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTa                                                                 <  6:232255/148‑1 (MQ=255)
                                                                                       cTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCGCCCAGACGAGTGTAGa                                                               <  2:327478/149‑1 (MQ=255)
                                                                                       cTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGa                                                               <  7:112851/149‑1 (MQ=255)
                                                                                             ttCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAATACCGGAGAAACGCCGGTGTAGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTc                                                           >  6:975/1‑147 (MQ=255)
                                                                                               ccAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAgg                                                       >  1:14651/1‑149 (MQ=255)
                                                                                                cAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAgg                                                       <  5:271587/148‑1 (MQ=255)
                                                                                                cAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCAAGgt                                                      >  1:282291/1‑149 (MQ=255)
                                                                                                       gTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGt                                               <  1:322037/149‑1 (MQ=255)
                                                                                                       gTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGt                                               >  6:63474/1‑149 (MQ=255)
                                                                                                        tAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgtcgtc                                                  >  1:304271/1‑145 (MQ=255)
                                                                                                        tAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgtcgtc                                                  <  2:304271/145‑1 (MQ=255)
                                                                                                            gATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACAATAAACGTTGTATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCATGTCGTCAGTGATTg                                          >  8:189369/1‑149 (MQ=255)
                                                                                                                      gTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAg                                <  8:307766/149‑1 (MQ=255)
                                                                                                                         aTCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACg                                                                       >  7:339301/1‑107 (MQ=255)
                                                                                                                         aTCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACg                                                                       <  8:339301/107‑1 (MQ=255)
                                                                                                                              gTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGGTTGGGTAACCCAGTTTAGCgg                        <  8:2784/149‑1 (MQ=255)
                                                                                                                               tACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATa                                                                                       >  5:251170/1‑85 (MQ=255)
                                                                                                                               tACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATa                                                                                       <  6:251170/85‑1 (MQ=255)
                                                                                                                               tACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGt                       >  2:151529/1‑149 (MQ=255)
                                                                                                                                aCTCAACACCGCCAGCGAAGTCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGt                       >  1:155046/1‑148 (MQ=255)
                                                                                                                                aCTCAACACCGCCAGCGAAGTCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGt                       >  1:150073/1‑148 (MQ=255)
                                                                                                                                 cTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGc                          <  5:46736/144‑1 (MQ=255)
                                                                                                                                 cTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGc                          >  6:46736/1‑144 (MQ=255)
                                                                                                                                    aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGa                                                                      <  1:264003/97‑1 (MQ=255)
                                                                                                                                    aaCACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGa                                                                      >  2:264003/1‑97 (MQ=255)
                                                                                                                                     acacCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCAt                                                                                        >  4:311593/1‑78 (MQ=255)
                                                                                                                                     acacCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCAt                                                                                        <  3:311593/78‑1 (MQ=255)
                                                                                                                                       acCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGaa               <  1:317488/149‑1 (MQ=255)
                                                                                                                                        ccGCCAGCGAAGTCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAAc              <  2:155046/149‑1 (MQ=255)
                                                                                                                                        ccGCCAGCGAAGTCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAAc              <  2:150073/149‑1 (MQ=255)
                                                                                                                                        ccGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAAc              <  7:639/149‑1 (MQ=255)
                                                                                                                                             aGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGAATTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGAGTAGATGTCCAGGACGTCAGTGATAGTGAAACCCAGATTAGCGGTCAGTTGAACGCCCt         >  4:287198/1‑149 (MQ=255)
                                                                                                                                                gAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACTCCATACCATGCCACCCAGACGAGTGAAGATCTCCATATACTCAGTAAAATTAAAAAACAGTTTAGCGGTCAGTTGAAAGCCCt         >  6:339672/1‑146 (MQ=255)
                                                                                                                                                  aGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCt    >  1:333874/1‑149 (MQ=255)
                                                                                                                                                    aCCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCttt  >  5:362923/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGAAACGGTAGGAAACACCCAGGCTCAGCATGCCGTTGTCCGGACGAGTGCCGATGGTGTGTGCGTCACCGATGTTGTTGGTCCACTGGTATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTT  >  NC_000913/1019478‑1019774

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