Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,486,985 A→C E353D (GAA→GAC)  tsgA → MFS transporter TsgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,486,9850AC83.8% 202.8 / 29.6 68E353D (GAA→GAC) tsgAMFS transporter TsgA
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base C (35/22);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGACCTCTCATTAATTGCTCA  >  NZ_CP009273/3486850‑3487115
                                                                                                                                       |                                                                                                                                  
cgcgATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCa                                                                                                                                  >  2:716865/1‑138 (MQ=255)
    aTCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCg                                                                                                                             >  1:862835/1‑139 (MQ=255)
     tCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCgg                                                                                                                            >  1:875563/1‑139 (MQ=255)
         taCCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGc                                                                                                                              <  1:26639/133‑1 (MQ=255)
               catcatCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcg                                                                                                                  >  2:289316/1‑139 (MQ=255)
               catcatCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcg                                                                                                                  >  1:784100/1‑139 (MQ=255)
               catcatCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcg                                                                                                                  >  1:56823/1‑139 (MQ=255)
               catcatCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcg                                                                                                                  >  1:421599/1‑139 (MQ=255)
               catcatCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcg                                                                                                                  >  1:375900/1‑139 (MQ=255)
                 tcatcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAggcgg                                                                                                                 >  2:544111/1‑138 (MQ=255)
                 tcatcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCg                                                                                                                             >  1:590133/1‑126 (MQ=255)
                 tcatcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCg                                                                                                                             <  2:590133/126‑1 (MQ=255)
                   atcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCAc                                                                                                              >  2:1179573/1‑139 (MQ=255)
                   atcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCAc                                                                                                              >  2:440496/1‑139 (MQ=255)
                   atcaCCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCa                                                                                                               <  1:1033559/138‑1 (MQ=255)
                        ccTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTcc                                                                                                                         >  1:918217/1‑123 (MQ=255)
                        ccTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTcc                                                                                                                         <  2:918217/123‑1 (MQ=255)
                                 acaGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAAc                                                                                                <  1:238339/139‑1 (MQ=255)
                                      agACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCt                                                                                           >  1:24322/1‑139 (MQ=255)
                                        aCCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAgg                                                                                                                    <  1:197950/112‑1 (MQ=255)
                                        aCCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAgg                                                                                                                    >  2:197950/1‑112 (MQ=255)
                                              gTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTg                                                                                   <  1:289316/139‑1 (MQ=255)
                                                    tCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCt                                                                             >  2:314244/1‑139 (MQ=255)
                                                    tCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCt                                                                             >  1:1176576/1‑139 (MQ=255)
                                                      gCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCttt                                                                           <  2:737233/139‑1 (MQ=255)
                                                       ccAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTtg                                                                          >  1:166430/1‑139 (MQ=255)
                                                          aaaCTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGa                                                                       <  1:1179573/139‑1 (MQ=255)
                                                           aaCTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGAt                                                                      <  1:716865/139‑1 (MQ=255)
                                                                gTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 <  1:544111/139‑1 (MQ=255)
                                                                gTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 <  2:875563/139‑1 (MQ=255)
                                                                 ttAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 <  1:249478/138‑1 (MQ=255)
                                                                 ttAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 >  2:249478/1‑138 (MQ=255)
                                                                   aaCTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTc                                                               <  2:862835/138‑1 (MQ=255)
                                                                        tGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGt                                                           >  1:1234807/1‑137 (MQ=255)
                                                                                    cGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGt                                                 <  2:1007494/135‑1 (MQ=255)
                                                                                    cGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGt                                                 >  1:1007494/1‑135 (MQ=255)
                                                                                    cGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGc                                             >  1:1138742/1‑139 (MQ=255)
                                                                                        aCCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtc                                         <  2:394181/139‑1 (MQ=255)
                                                                                              ggTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcaccgt                                    <  2:150578/138‑1 (MQ=255)
                                                                                              ggTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcaccgt                                    <  2:24322/138‑1 (MQ=255)
                                                                                                     tGTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATaa                            >  1:490418/1‑139 (MQ=255)
                                                                                                     tGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCAt                              <  2:1234807/137‑1 (MQ=255)
                                                                                                      gTTGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATaa                            >  2:328059/1‑138 (MQ=255)
                                                                                                       ttGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcacc                                      <  2:1034425/127‑1 (MQ=255)
                                                                                                       ttGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcacc                                      <  2:998148/127‑1 (MQ=255)
                                                                                                       ttGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcacc                                      >  1:1034425/1‑127 (MQ=255)
                                                                                                       ttGACCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtcacc                                      >  1:998148/1‑127 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGc                                                                                     <  1:132768/77‑1 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGc                                                                                     >  2:132768/1‑77 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 >  2:952685/1‑97 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCt                                                                 <  1:952685/97‑1 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCAt                              >  1:1061972/1‑132 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCAt                              <  2:1061972/132‑1 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  2:536916/1‑139 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  2:1205197/1‑139 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  2:677093/1‑139 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  2:1212222/1‑139 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  1:807322/1‑139 (MQ=255)
                                                                                                          aCCTTTGTGGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACAccc                       >  2:1008531/1‑139 (MQ=255)
                                                                                                               tgtgGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTTTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCAACGTCAGCATAACACCCTGAcc                  >  2:510814/1‑139 (MQ=255)
                                                                                                               tgtgGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATa                             >  2:704130/1‑128 (MQ=255)
                                                                                                               tgtgGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATa                             <  1:704130/128‑1 (MQ=255)
                                                                                                               tgtgGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGAcc                  >  2:898641/1‑139 (MQ=255)
                                                                                                               tgtgGTTACCGGCCCGATCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGAc                   <  1:257921/138‑1 (MQ=255)
                                                                                                                      aCCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGACCTCTCAt            >  1:1145745/1‑138 (MQ=255)
                                                                                                                               aTCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtc                                         <  1:358170/100‑1 (MQ=255)
                                                                                                                               aTCGTTGACCATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGccgtc                                         >  2:358170/1‑100 (MQ=255)
                                                                                                                               aTCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGACCTCTCATTAATTGCTCa  >  2:282462/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CGCGATCTATACCACCATCATCACCCTGGGTTCACAGCAGACCAAAGTACCGTCGCCAAAACTGGTTAACTTTGTCCTGACCTGCGGAACCATCGGTACTATGTTGACCTTTGTGGTTACCGGCCCGATCGTTGAACATAGCGGTCCGCAGGCGGCACTGCTTACAGCAAACGGTCTGTACGCTGTCGTCTTTGTGATGTGCTTCCTGTTAGGTTTCGTCAGCCGTCACCGTCAGCATAACACCCTGACCTCTCATTAATTGCTCA  >  NZ_CP009273/3486850‑3487115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑