Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,453,684 T→A 8.5% I51F (ATT→TTT)  ppiC‑I ← peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,453,6840TA 8.5% 166.6 / 5.2 59I51F (ATT→TTT) ppiC‑Ipeptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base A (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

TCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCGGGCTTGGGCTTT  >  NC_002947/3453536‑3453831
                                                                                                                                                    |                                                                                                                                                    
tCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAAt                                                                                                                                                      <  2:245838/149‑1 (MQ=255)
   ggTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTc                                                                                                                                                   >  5:35006/1‑149 (MQ=255)
   ggTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTc                                                                                                                                                   >  5:29474/1‑149 (MQ=255)
         aaaGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAAATTCGccccg                                                                                                                                             >  1:110741/1‑147 (MQ=255)
          aaGGGCTAGCGCGTTAGTAGCGGTAGAACGTTTGCACCAGGTTGTAGGCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGACGACCATCTGCCCGGGGCGAATTTCGCCCAgg                                                                                                                                            <  7:125617/149‑1 (MQ=255)
           aGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAAATTCGCCCCggg                                                                                                                                           <  2:110741/149‑2 (MQ=255)
                 aGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTcacccc                                                                                                                                     >  2:138685/1‑147 (MQ=255)
                  gcgcGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGc                                                                                                                                                 <  2:59457/136‑1 (MQ=255)
                  gcgcGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGc                                                                                                                                                 >  1:59457/1‑136 (MQ=255)
                     cgTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCTCCAAGGTCACAACcgcg                                                                                                                                 >  3:46412/1‑149 (MQ=255)
                         aGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTg                                                                                                                             >  2:267001/1‑149 (MQ=255)
                          gTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTg                                                                                                                             <  3:50254/148‑1 (MQ=255)
                             gcgGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAAATTCGccc                                                                                                                                               <  8:320906/127‑1 (MQ=255)
                                gTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATg                                                                                                                      <  3:109865/149‑1 (MQ=255)
                                    aaCGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGc                                                                                                                   >  2:68432/1‑148 (MQ=255)
                                       gTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGc                                                                                                                                                 >  2:248482/1‑115 (MQ=255)
                                       gTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGAGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCa                                                                                                                  >  3:7446/1‑146 (MQ=255)
                                            cACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAagagag                                                                                                           <  6:23956/148‑1 (MQ=255)
                                             aCCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCg                                                                                                                                                  >  5:309333/1‑108 (MQ=255)
                                             aCCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCg                                                                                                                                                  <  6:309333/108‑1 (MQ=255)
                                             aCCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGaa                                                                                                         <  1:68432/149‑1 (MQ=255)
                                              ccAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAAc                                                                                                        <  6:347188/149‑1 (MQ=255)
                                               cAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGATTCTTGAAGATCACCTGGTCGATGGAGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACt                                                                                                       >  1:210888/1‑149 (MQ=255)
                                                aGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGCACtt                                                                                                      >  4:329011/1‑149 (MQ=255)
                                                  gtggtAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACtt‑tt                                                                                                   >  3:134234/1‑149 (MQ=255)
                                                       aGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTcgc                                                                                               >  2:118860/1‑148 (MQ=255)
                                                           gAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCCTCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCa                                                                                          >  6:5424/1‑149 (MQ=255)
                                                              tttGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGt                                                                                       <  2:188006/149‑1 (MQ=255)
                                                                   tctTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGa                                                                                    >  3:165201/1‑147 (MQ=255)
                                                                     ttGATCGGCCCGTGTACCACACGTAAAGGCTTCTTGAAGATCACCTGGTCAATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGc                                                                                >  1:284563/1‑149 (MQ=255)
                                                                       gATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCt                                                                              >  6:67240/1‑149 (MQ=255)
                                                                           ggCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCACAGGTCACCACCGCGCTTGCCGGATGGGAAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGc                                                                           <  3:313393/148‑1 (MQ=255)
                                                                            gCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCg                                                                         >  2:8444/1‑149 (MQ=255)
                                                                              ccGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGCGCAATT‑TATCGCGACCAAGTCCAGCGCCTCCCCGtt                                                                       >  6:289271/1‑149 (MQ=255)
                                                                                      cacaCGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATgc                                                               <  1:26257/149‑1 (MQ=255)
                                                                                       acacGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATgcg                                                              >  8:192228/1‑149 (MQ=255)
                                                                                           gCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCgctgc                                                          <  6:29474/149‑1 (MQ=255)
                                                                                           gCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCgctgc                                                          <  6:35006/149‑1 (MQ=255)
                                                                                             aaaGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCtt                                                        >  6:127580/1‑149 (MQ=255)
                                                                                                 gCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGc                                                    <  1:267001/149‑1 (MQ=255)
                                                                                                 gCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGc                                                    >  6:244186/1‑149 (MQ=255)
                                                                                                    tcttGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCttt                                                 >  6:111155/1‑149 (MQ=255)
                                                                                                      ttGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTc                                               >  5:330542/1‑149 (MQ=255)
                                                                                                        gAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGGCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTc                                               <  1:138685/147‑1 (MQ=255)
                                                                                                        gAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCgg                                             >  5:278403/1‑149 (MQ=255)
                                                                                                        gAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCg                                              <  6:330542/148‑1 (MQ=255)
                                                                                                        gAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCg                                              <  5:111155/148‑1 (MQ=255)
                                                                                                            aTCACCTGGTCGATGGCGCCGACCATCTGCGCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTc                                         <  4:134234/149‑1 (MQ=255)
                                                                                                              cACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACGT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGt                                       >  7:174836/1‑149 (MQ=255)
                                                                                                                       ttttgcgcCGAGCATCTGTGCGGGGGGATTATGGCCCAGGTGACACCCGCGCGTGGGGGTGTTGCTCAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCt                              <  5:289271/145‑1 (MQ=255)
                                                                                                                       gATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTc                               <  7:192228/148‑1 (MQ=255)
                                                                                                                             gcCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCa                        <  5:127580/149‑1 (MQ=255)
                                                                                                                                   cATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTTCTTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCgctgc                                                          <  1:263101/110‑1 (MQ=255)
                                                                                                                                   cATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCgctgc                                                          >  2:263101/1‑109 (MQ=255)
                                                                                                                                         cccGGGGCGAATTTCGCCCAGGGCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGGCTTGACCAGAATGTGGCg              <  5:5424/147‑1 (MQ=255)
                                                                                                                                          cctgggCGAAATTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCgggc           >  1:348656/4‑149 (MQ=255)
                                                                                                                                           ctgggCGAAATTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGATGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGCCCAGAATGTGGCgggct          >  8:188234/3‑149 (MQ=255)
                                                                                                                                              ggCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCgggcttgg       <  4:165201/149‑1 (MQ=255)
                                                                                                                                              ggCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCgggcttgg       <  2:284563/149‑1 (MQ=255)
                                                                                                                                               gCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCgggcttggg      <  5:244186/149‑1 (MQ=255)
                                                                                                                                                    tttCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCGGGCTTGGGCttt  >  6:327717/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATCGGCCCGTGTACCACACGCAAAGGCTTCTTGAAGATCACCTGGTCGATGGCGCCGACCATCTGCCCGGGGCGAATTTCGCCCAGGTCACCACCGCGCTTGCCGGATGGGCAAAGAGAGAACTT‑TTTCGCGAGCACGTCGAAAGCCTCGCCGTTGGCGATGCGCTGCTTGAGCTTTTCGGCTTCGTCAGCGGTCTTGACCAGAATGTGGCGGGCTTGGGCTTT  >  NC_002947/3453536‑3453831

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