Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,867,513 A→G D422D (GAT→GAC cysN ← sulfate adenylyltransferase subunit CysN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,867,5130AG96.9% 222.2 / 0.1 66coding (1266/1428 nt)cysNsulfate adenylyltransferase subunit CysN
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (26/37);  minor base . (1/1);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTA  >  NZ_CP009273/2867380‑2867648
                                                                                                                                     |                                                                                                                                       
cccAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGcggtcgat                                                                                                                                    <  1:615811/139‑1 (MQ=255)
        ggAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAAt                                                                                                                              >  2:558192/1‑137 (MQ=255)
          aaaGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGc                                                                                                                          >  2:213481/1‑139 (MQ=255)
             gTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCa                                                                                                                       <  2:443065/139‑1 (MQ=255)
              tGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCAc                                                                                                                      <  2:462210/139‑1 (MQ=255)
               ggCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGcccccc                                                                                                                     >  2:1035701/1‑139 (MQ=255)
               ggCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCAcc                                                                                                                     >  2:51341/1‑139 (MQ=255)
                    cgaACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCa                                                                                                                <  1:1210938/139‑1 (MQ=255)
                    cgaACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCG‑‑CGATAAAAATCAGCCCACCCGTCAcc                                                                                                              >  2:326189/1‑139 (MQ=255)
                            agCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCAc                                                                                                               <  2:1257473/132‑1 (MQ=255)
                            agCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCAc                                                                                                               >  1:1257473/1‑132 (MQ=255)
                              cATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTttg                                                                                                      <  2:35976/139‑1 (MQ=255)
                                   aaTTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGAt                                                                                                 <  1:1176417/139‑1 (MQ=255)
                                     ttCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               <  2:304843/139‑1 (MQ=255)
                                          gTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATc                                                                                          >  1:422019/1‑139 (MQ=255)
                                           ttCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCt                                                                                         <  1:527627/139‑1 (MQ=255)
                                                            cAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGt                                                                         <  1:1269400/138‑1 (MQ=255)
                                                                  gcgcAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGt                                                                  >  2:85124/1‑139 (MQ=255)
                                                                             gCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCa                                                                       >  2:376778/1‑123 (MQ=255)
                                                                             gCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCa                                                                       <  1:376778/123‑1 (MQ=255)
                                                                             gCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCAc                                                         >  2:1205643/1‑137 (MQ=255)
                                                                             gCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGCGAGATCCAc                                                         >  1:577655/1‑137 (MQ=255)
                                                                                               tGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               >  1:565078/1‑81 (MQ=255)
                                                                                               tGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               <  2:565078/81‑1 (MQ=255)
                                                                                                 cACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGc                                   <  1:596632/139‑1 (MQ=255)
                                                                                                 cACCATACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGc                                   <  2:471141/139‑1 (MQ=255)
                                                                                                      tACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTgag                                                               >  1:106071/1‑106 (MQ=255)
                                                                                                      tACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTgag                                                               <  2:106071/106‑1 (MQ=255)
                                                                                                      tACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    <  1:1106436/117‑1 (MQ=255)
                                                                                                      tACCGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    >  2:1106436/1‑117 (MQ=255)
                                                                                                        ccGGCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGtttt                            <  2:929880/139‑1 (MQ=255)
                                                                                                           gCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               >  2:91737/1‑69 (MQ=255)
                                                                                                           gCACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               <  1:91737/69‑1 (MQ=255)
                                                                                                            cACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               >  2:908156/1‑68 (MQ=255)
                                                                                                            cACCCACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               <  1:908156/68‑1 (MQ=255)
                                                                                                              cccACGGTCACATTGCTCAGGCGGTCGATAAATATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCg                                                  >  1:1124849/1‑111 (MQ=255)
                                                                                                              cccACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCAc                                                         <  1:1091993/104‑1 (MQ=255)
                                                                                                              cccACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCAc                                                         >  2:1091993/1‑104 (MQ=255)
                                                                                                              cccACGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCg                                                  <  2:1124849/111‑1 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   >  1:625210/1‑139 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   <  2:625210/139‑1 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   >  1:811089/1‑139 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   <  1:567326/139‑1 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   <  1:51341/139‑1 (MQ=255)
                                                                                                                 aCGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCAc                   <  1:1205643/139‑1 (MQ=255)
                                                                                                                  cGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACg                  <  1:85124/139‑1 (MQ=255)
                                                                                                                  cGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACg                  <  1:360891/139‑1 (MQ=255)
                                                                                                                  cGGTCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACg                  <  2:395136/139‑1 (MQ=255)
                                                                                                                     tCACATTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTg               >  2:291573/1‑139 (MQ=255)
                                                                                                                       acaTTGCTCAGGCG‑‑CGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGt            <  1:326189/138‑1 (MQ=255)
                                                                                                                        caTTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGt            <  2:773773/139‑1 (MQ=255)
                                                                                                                        caTTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGt            <  2:877302/139‑1 (MQ=255)
                                                                                                                        caTTGCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGt            <  2:422019/139‑1 (MQ=255)
                                                                                                                        caTTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGt            <  2:984913/139‑1 (MQ=255)
                                                                                                                            gCTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAg         <  1:213481/138‑1 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    <  1:1048316/94‑1 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    >  2:1122970/1‑94 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    >  2:1048316/1‑94 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGc                                                    <  1:1122970/94‑1 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACg                  >  1:967737/1‑128 (MQ=255)
                                                                                                                             cTCAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACg                  <  2:967737/128‑1 (MQ=255)
                                                                                                                               cAGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGtta     >  1:1021631/1‑139 (MQ=255)
                                                                                                                                aGGCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGttat    <  1:291573/139‑1 (MQ=255)
                                                                                                                                 ggCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               >  2:526772/1‑47 (MQ=255)
                                                                                                                                 ggCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATaa                                                                                               <  1:526772/47‑1 (MQ=255)
                                                                                                                                 ggCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGttat    >  1:264067/1‑138 (MQ=255)
                                                                                                                                  gCGGTCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGttatta  <  2:1022341/139‑1 (MQ=255)
                                                                                                                                     atcgatAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCAc                                                         <  2:24822/81‑1 (MQ=255)
                                                                                                                                     atcgatAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCAc                                                         >  1:24822/1‑81 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTA  >  NZ_CP009273/2867380‑2867648

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