Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,128,232 A→T 7.7% Q1201H (CAA→CAT)  PP_1887 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,128,2320AT 7.7% 181.1 / 5.3 65Q1201H (CAA→CAT) PP_1887hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

CCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGA  >  NC_002947/2128090‑2128375
                                                                                                                                              |                                                                                                                                                 
ccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGc                                                                                                                                           >  2:103853/1‑149 (MQ=255)
ccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGc                                                                                                                                           >  6:3308/1‑149 (MQ=255)
      aCCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATT‑‑TGTGCCGGATc                                                                                                                                     <  6:282223/149‑1 (MQ=255)
        cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGCCCCGACGGGTCGTGGCAGCGGCCCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCa                                                                                                                                    >  8:218638/1‑148 (MQ=255)
        cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCAt                                                                                                                                   >  2:131103/1‑149 (MQ=255)
         ccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTTGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCAtt                                                                                                                                  <  3:263426/149‑1 (MQ=255)
           tCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTcc                                                                                                                                >  1:2510/1‑149 (MQ=255)
               gTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCt                                                                                                                            <  2:2510/149‑1 (MQ=255)
                  aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGa                                                                                                                          <  2:77035/148‑1 (MQ=255)
                   aaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATg                                                                                                                        <  5:3308/149‑1 (MQ=255)
                      cATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATTGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATgcgc                                                                                                                     <  5:220407/149‑1 (MQ=255)
                           ccAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                                <  2:164776/149‑1 (MQ=255)
                           ccAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                                <  3:71509/149‑1 (MQ=255)
                             aGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                                >  6:79272/1‑147 (MQ=255)
                               gCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGTGCTACTGgaga                                                                                                            >  2:77474/1‑149 (MQ=255)
                               gCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGgaga                                                                                                            <  1:131103/149‑1 (MQ=255)
                                cAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGAc                                                                                                           <  1:334213/149‑1 (MQ=255)
                                  gATTTTCGGCTTCGATTGCATGGACAACATCCTCTCGGTGATCACCTACAGACCTGAAGGGTCGTGGCAGCGGACCTTGTACGAATTATAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCt                                                                                                         <  5:35843/149‑1 (MQ=255)
                                   atatTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTATGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTa                                                                                                        >  2:65317/1‑149 (MQ=255)
                                     atTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGcc                                                                                                                             >  4:65369/1‑126 (MQ=255)
                                     atTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGcc                                                                                                                             <  3:65369/126‑1 (MQ=255)
                                       tCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATT‑‑TGTGCCGGAt                                                                                                                                      >  4:186619/1‑115 (MQ=255)
                                       tCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATT‑‑TGTGCCGGAt                                                                                                                                      <  3:186619/115‑1 (MQ=255)
                                        cGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGAc                                                                                                    <  7:142110/148‑1 (MQ=255)
                                               gATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACaa                                                                                                                                                   <  8:325558/96‑1 (MQ=255)
                                               gATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACaa                                                                                                                                                   >  7:325558/1‑96 (MQ=255)
                                                          caacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                                  >  6:278562/1‑148 (MQ=255)
                                                             caTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgt                                                                              >  3:30451/1‑149 (MQ=255)
                                                             caTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgt                                                                              >  1:2474/1‑149 (MQ=255)
                                                                    tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAATAATAGCGATCCGGCACAAATCTTTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgt                                                                           <  1:103853/147‑1 (MQ=255)
                                                                      ggTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGTGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCa                                                                        <  1:77474/146‑1 (MQ=255)
                                                                      ggTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCa                                                                        <  5:278562/146‑1 (MQ=255)
                                                                           tCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                                <  2:6169/149‑1 (MQ=255)
                                                                            cACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGATGTTAAAGATCGTc                                                                >  2:142619/1‑148 (MQ=255)
                                                                             aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGAc                                                                                                           <  8:302336/104‑1 (MQ=255)
                                                                             aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGAc                                                                                                           >  7:302336/1‑104 (MQ=255)
                                                                             aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCg                                                                  <  3:136778/145‑1 (MQ=255)
                                                                             aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCg                                                                  >  4:136778/1‑145 (MQ=255)
                                                                             aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                               >  2:39133/1‑148 (MQ=255)
                                                                              ccTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGgaggag                                                                               <  4:145/131‑1 (MQ=255)
                                                                                 aCAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATgc                                                                                                                       <  1:72531/88‑1 (MQ=255)
                                                                                 aCAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATgc                                                                                                                       >  2:72531/1‑88 (MQ=255)
                                                                                  cAGACGGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCATATACg                                                         >  5:239820/1‑149 (MQ=255)
                                                                                                       gCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGtt                                               <  2:12721/138‑1 (MQ=255)
                                                                                                       gCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGtt                                               >  1:12721/1‑138 (MQ=255)
                                                                                                         ggACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTg                                  >  5:175027/1‑149 (MQ=255)
                                                                                                          gACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAC‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGc                                 <  4:214939/149‑1 (MQ=255)
                                                                                                           aCCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggagga                                                                                >  8:313857/1‑101 (MQ=255)
                                                                                                           aCCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggagga                                                                                <  7:313857/101‑1 (MQ=255)
                                                                                                                 tACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATTCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacg                           <  5:79272/148‑1 (MQ=255)
                                                                                                                  aCGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTgg                                             <  5:162191/129‑1 (MQ=255)
                                                                                                                  aCGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTgg                                             >  6:162191/1‑129 (MQ=255)
                                                                                                                          aaGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATAGGAACCGCAGTTGGCTATCGATTTGCATTACGACGGCAATgg                 >  3:91727/1‑149 (MQ=255)
                                                                                                                          aaGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATgg                 >  4:124160/1‑149 (MQ=255)
                                                                                                                          aaGAATAGCGATCCGGCACAAA‑‑TGAGACGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                               >  3:103079/1‑103 (MQ=255)
                                                                                                                          aaGAATAGCGATCCGGCACAAA‑‑TGAGACGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                               <  4:103079/103‑1 (MQ=255)
                                                                                                                          aaGAAAAGAGTTCCGGAACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATg                  <  6:239820/148‑1 (MQ=255)
                                                                                                                             aaTAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATg                  <  8:173677/145‑1 (MQ=255)
                                                                                                                             aaTAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGGTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATg                  >  7:173677/1‑145 (MQ=255)
                                                                                                                              ggaatgATCCGGCTCATT‑‑TGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAAc             >  3:250038/6‑149 (MQ=255)
                                                                                                                              aTAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAAc             <  4:30451/149‑1 (MQ=255)
                                                                                                                                 gCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCAGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATgctgcgcg           >  1:122914/1‑141 (MQ=255)
                                                                                                                                   gATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt        >  2:26631/1‑149 (MQ=255)
                                                                                                                                    aTCCGGCTCATT‑‑TGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt        >  1:93369/1‑148 (MQ=255)
                                                                                                                                       cGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGt    <  1:65317/149‑1 (MQ=255)
                                                                                                                                       cGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGt    >  2:18357/1‑149 (MQ=255)
                                                                                                                                         gCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGtga  <  1:18357/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAA‑‑TGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGA  >  NC_002947/2128090‑2128375

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: ‑