Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,057,664 T→C G95G (GGA→GGG cobT ← nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,057,6640TC96.7% 325.1 / 1.2 92G95G (GGA→GGGcobTnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (49/40);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

GGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACA  >  NZ_CP009273/2057577‑2057750
                                                                                       |                                                                                      
ggggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTccc                                                                                      <  2:2848860/90‑1 (MQ=255)
 gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCa                                                                                     >  2:1681929/1‑90 (MQ=255)
 gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCa                                                                                     >  1:2043180/1‑90 (MQ=255)
 gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCa                                                                                     >  1:1960080/1‑90 (MQ=255)
  ggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCAc                                                                                    <  1:2993046/90‑1 (MQ=255)
    aTAGGCTCAGCAGTATCAATACCAACATCAATTACGTCGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACgg                                                                                  >  1:2657906/1‑90 (MQ=255)
      aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGt                                                                                >  1:2770156/1‑90 (MQ=255)
      aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGt                                                                                >  2:813612/1‑90 (MQ=255)
       ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTc                                                                               <  2:1765086/90‑1 (MQ=255)
       ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCAGGTTGTCCCACGGGTc                                                                               >  2:1688578/1‑90 (MQ=255)
         cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCat                                                                             <  2:2346402/90‑1 (MQ=255)
           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCTACGGGTCatat                                                                           >  2:2225403/1‑90 (MQ=255)
           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                           >  2:1780382/1‑90 (MQ=255)
           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                           >  1:2256939/1‑90 (MQ=255)
             gcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATAttt                                                                         >  2:2276224/1‑90 (MQ=255)
                gTATCAATACCAACATCAATTACGTGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     <  1:161141/90‑1 (MQ=255)
                 tATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     <  1:2998691/90‑1 (MQ=255)
                 tATCAATACCAACATCAATTACCTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     <  2:1357653/90‑1 (MQ=255)
                  aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGc                                                                    >  2:1658460/1‑90 (MQ=255)
                  aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGc                                                                    >  1:1142962/1‑90 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   <  2:2807465/90‑1 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   <  1:1158646/90‑1 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   <  2:2942101/90‑1 (MQ=255)
                     aaTACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTg                                                                 >  1:288217/1‑90 (MQ=255)
                     aaTACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTg                                                                 >  2:785657/1‑90 (MQ=255)
                     aaTACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGGCATATTTTCAGCCTg                                                                 >  2:1615677/1‑90 (MQ=255)
                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                               >  1:2052276/1‑90 (MQ=255)
                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                            >  1:2000033/1‑90 (MQ=255)
                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                            >  1:1203964/1‑90 (MQ=255)
                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                            >  1:1181979/1‑90 (MQ=255)
                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                            >  1:940153/1‑90 (MQ=255)
                           aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   >  2:2390200/1‑82 (MQ=255)
                           aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   <  1:2390200/82‑1 (MQ=255)
                           aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGc                                                           <  1:2822701/90‑1 (MQ=255)
                           aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGc                                                           <  1:2654247/90‑1 (MQ=255)
                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                          >  2:477196/1‑90 (MQ=255)
                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                          <  1:534433/90‑1 (MQ=255)
                                cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                      >  1:1083707/1‑90 (MQ=255)
                                cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                      >  1:2251867/1‑90 (MQ=255)
                                 aaTTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAc                                                     >  2:579151/1‑90 (MQ=255)
                                  aTTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                    <  2:665606/90‑1 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     <  2:664784/68‑1 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     >  1:664784/1‑68 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   >  2:873303/1‑70 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                   <  1:873303/70‑1 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                               >  1:1417439/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                               >  1:614547/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                               >  2:54645/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                               >  2:1138448/1‑90 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                              <  1:2655997/90‑1 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                              >  1:383019/1‑90 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                              <  1:2066853/90‑1 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                              <  1:1921427/90‑1 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                            >  1:2747417/1‑90 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                            >  1:86307/1‑90 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                            >  1:2874491/1‑90 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                            >  1:966493/1‑90 (MQ=255)
                                               tAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGaaa                                       >  1:882098/1‑90 (MQ=255)
                                                 gcgcCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     >  2:437379/1‑58 (MQ=255)
                                                 gcgcCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                     <  1:437379/58‑1 (MQ=255)
                                                  cgcCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAg                                    <  2:2747417/90‑1 (MQ=255)
                                                   gcCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGc                                   >  1:845003/1‑90 (MQ=255)
                                                    cccGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCg                                  >  1:1903271/1‑90 (MQ=255)
                                                     ccGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGa                                 <  2:298464/90‑1 (MQ=255)
                                                       gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               <  2:2770156/90‑1 (MQ=255)
                                                       gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               >  1:2390465/1‑90 (MQ=255)
                                                           gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               >  1:696422/1‑86 (MQ=255)
                                                           gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               <  2:696422/86‑1 (MQ=255)
                                                           gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                           <  2:1120435/90‑1 (MQ=255)
                                                           gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                           <  1:1208867/90‑1 (MQ=255)
                                                                gCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTtcctcc                      >  1:290674/1‑90 (MQ=255)
                                                                         cacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               >  1:2835960/1‑72 (MQ=255)
                                                                         cacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                               <  2:2835960/72‑1 (MQ=255)
                                                                         cacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCg             >  1:1783800/1‑90 (MQ=255)
                                                                          acacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGt            >  2:287940/1‑90 (MQ=255)
                                                                             cGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGAt         >  1:1839935/1‑90 (MQ=255)
                                                                              gCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc        >  2:1529368/1‑90 (MQ=255)
                                                                              gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc        <  2:1783800/90‑1 (MQ=255)
                                                                                 ggTTGTCCCCCGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:786001/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  1:1103531/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:86307/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:1960080/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:1181979/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:1203964/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  1:579151/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:288217/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:2842822/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  1:2276224/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  1:477196/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:2727051/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     <  2:2631185/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     >  1:2842822/1‑90 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc     >  2:2342874/1‑90 (MQ=255)
                                                                                   ttgtcCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGcac   >  1:1094624/1‑90 (MQ=255)
                                                                                    tgtcCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGcaca  <  2:1142962/90‑1 (MQ=255)
                                                                                       |                                                                                      
GGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACA  >  NZ_CP009273/2057577‑2057750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: