Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 711,536 C→T intergenic (‑71/+69) fldA ← / ← ybfE flavodoxin 1/LexA‑regulated protein, CopB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913711,5360CT96.7% 180.3 / NA 62intergenic (‑71/+69)fldA/ybfEflavodoxin 1/LexA‑regulated protein, CopB family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/37);  minor base G (0/2);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGC  >  NC_000913/711396‑711669
                                                                                                                                            |                                                                                                                                      
ttttAACCATTTTTGCGATATTTTCAGTATTCCCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCAtagtag                                                                                                                                      <  2:48073/143‑1 (MQ=255)
tttGAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCAtagtag                                                                                                                                      >  1:971987/1‑143 (MQ=255)
   gattCATTTTTGCGATGGTTGCGCTATGAACGGAGCCGCTGCCGGAATGGATGCCAGTGATAGCCATGAGTGAAATAATCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAAc                                                                                                                                   <  2:133377/140‑1 (MQ=255)
   gAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAAc                                                                                                                                   <  1:329929/143‑1 (MQ=255)
    aaTCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAAGCATAGTAGAACa                                                                                                                                  >  1:622040/1‑143 (MQ=255)
      tCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGAAGAACAgg                                                                                                                                >  2:999642/1‑143 (MQ=255)
                gttatTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGTATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCg                                                                                                                      >  1:55496/3‑143 (MQ=255)
                gATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCg                                                                                                                      >  1:432784/1‑143 (MQ=255)
                 atatTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCgg                                                                                                                     >  1:1000882/1‑143 (MQ=255)
                    ttttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATAAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGGCTTGCGGAAt                                                                                                                  >  1:1035248/1‑143 (MQ=255)
                    tttgCGGTATCACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTGATCGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAAt                                                                                                                  <  2:326564/143‑1 (MQ=255)
                     tttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATc                                                                                                                 <  1:1193181/143‑1 (MQ=255)
                      ttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCa                                                                                                                <  1:247168/143‑1 (MQ=255)
                      ttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCa                                                                                                                <  1:25541/143‑1 (MQ=255)
                       tCGGTATCACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAg                                                                                                               <  1:986729/143‑1 (MQ=255)
                        cGGTATTGCCGGTGTCGCTGCCGGAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGTAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGc                                                                                                              <  2:293441/143‑1 (MQ=255)
                        cGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGc                                                                                                              <  1:999497/143‑1 (MQ=255)
                          gTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGACAATCATAGTAGAACAGGCAGTCGCGCGGAATCAGCaa                                                                                                            >  2:889948/1‑143 (MQ=255)
                              tACCTGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGa                                                                                                        <  2:514230/143‑1 (MQ=255)
                                 cGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCa                                                                                                     <  2:1034250/143‑1 (MQ=255)
                                 cGGGGTCGCTGCCGAAAAAGATGCCAGTGAGAGCCGTGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGACCAGGCAGTCTTGCGGAATCAGCCAACGAGcc                                                                                                     <  2:29388/143‑2 (MQ=255)
                                    tgtCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGa                                                                                                  <  2:1003180/143‑1 (MQ=255)
                                         gtGCGGAAAATGATGCCAGGTGTAGCCATGTGTTGAATACTCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGCCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGca                                                                                             <  2:390087/142‑1 (MQ=255)
                                         cTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGca                                                                                             <  1:898118/143‑1 (MQ=255)
                                            ccGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCacac                                                                                          >  1:888837/1‑143 (MQ=255)
                                              gAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCACACGAGCAGGACTGCACACtg                                                                                        >  1:1136702/1‑143 (MQ=255)
                                              gAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACtg                                                                                        >  1:857926/1‑143 (MQ=255)
                                               aaaaaGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATACAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACtgt                                                                                       >  1:604596/1‑143 (MQ=255)
                                                aaaaGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACtgtg                                                                                      <  1:10830/143‑1 (MQ=255)
                                                  aaGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCt                                                                                    <  1:423518/143‑1 (MQ=255)
                                                    gATGCCAGTGATAGCCATGGGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGTGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTAc                                                                                  <  2:254726/143‑1 (MQ=255)
                                                              aTAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTgg                                                                       >  1:171568/1‑143 (MQ=255)
                                                                    aTGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAAttt                                                                 <  1:414948/143‑1 (MQ=255)
                                                                    aTGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAAATGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTCGCGGAATCAGCAAACGAGCAGGCCTGCACACTATGCTACATGAAAGTGGaacatt                                                                 >  2:789341/1‑139 (MQ=255)
                                                                       gtgcaaaTAAGCGCTTGGAACCTACTGAAATGGGGGTGGAAGATTGGCCACGGATAAAGGCAATCGTAGTAGACCAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                              <  2:839783/139‑1 (MQ=255)
                                                                       ggTGAAATACGCTCTTGGAACTTATTGAAAGGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGTTGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                              <  2:493250/142‑1 (MQ=255)
                                                                       cgTGAAATATCCTCTTGAAACTTCTTTAAATGGGGGTGGACAACTGCCCACGGATAAATGCAATCATAGTAGCACAGGCAGTCTTGCGGAATCAGCCAACGAGCATGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                              <  2:145070/142‑1 (MQ=255)
                                                                         tGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACg                                                            >  2:983453/1‑143 (MQ=255)
                                                                           aaaTAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGAt                                                          <  2:432784/143‑1 (MQ=255)
                                                                           aaaTAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGAt                                                          >  2:919605/1‑143 (MQ=255)
                                                                             aTAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCATACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAAGTTCAACGATGc                                                        >  2:137401/1‑143 (MQ=255)
                                                                                accttctGATACGTATTGAAATGGGGGTGGAAAATTGCGCACGGATAAGGGCAATAATAGTAGCAGAGGGAGTTTTGTGGTGGATTCAAACAATCAGGACTGCACTCTGTACTACATGGAAGTGGAAATTTAAACGATGCCCt                                                     <  2:531935/137‑1 (MQ=255)
                                                                                  ctctTGAAACTTATTGAAATGGCGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGa                                                   <  2:1115366/143‑1 (MQ=255)
                                                                                    ctTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACt                                                 <  1:668000/143‑1 (MQ=255)
                                                                                            ttATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCGTGCGAACTAGAGGGACGTTCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAgcg                                         >  1:1194326/1‑143 (MQ=255)
                                                                                               ttGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAACCGAGCAGGACTGCACACTGTGCTACCTGAAAGTGGCACTTTACACGATGCCCTGACTACGCAGcgccg                                      >  1:244244/1‑143 (MQ=255)
                                                                                                 gAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCgccgcc                                    >  1:373823/1‑143 (MQ=255)
                                                                                                  aaaTGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTGGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCa                                   >  1:743366/1‑143 (MQ=255)
                                                                                                             gAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCAtt                        >  1:651893/1‑143 (MQ=255)
                                                                                                              aaaaTTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGCCTACGCAGCGCCGCCAGTTGCTGCATTa                       <  1:333775/143‑1 (MQ=255)
                                                                                                              aaaaGTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTa                       <  2:402193/143‑1 (MQ=255)
                                                                                                               aaaTTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTaa                      <  2:285002/143‑1 (MQ=255)
                                                                                                               aaaTTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACGGCACAAAGAACAACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTaa                      <  1:45600/143‑1 (MQ=255)
                                                                                                                agttGTCCACGGATAAAGGGCATCATAGTCGGACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGTCTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAAc                     <  2:34546/141‑1 (MQ=255)
                                                                                                                        aCGGATAAAGGCAATCATAGTAGACCAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTc             <  1:441339/143‑1 (MQ=255)
                                                                                                                          ggATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATCAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGa           <  1:477132/143‑1 (MQ=255)
                                                                                                                           gATAAAGGCAACCATAGTAGAACAGGCAGTCCTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAt          <  1:1144197/143‑1 (MQ=255)
                                                                                                                           gATAAAGGCAACCATAGGAGGACGGGGAGTCGTGCGGAAGCAGCAAGCGAGTAGGACCGCACACTGGGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAt          <  1:800529/143‑1 (MQ=255)
                                                                                                                             tAAAGGCAATCATAGTAGAACAGGCAGTCCTGCGGAATCAGCACACCAGCAGGACTGCTCCCTCCCCACGTCGGTTTTGGCTGTGCAGACCAGGCCGTGACTACGAAGAGACGTCAGTTGCTGCATTAACATCTCTTCGATCa        <  2:1167391/143‑1 (MQ=255)
                                                                                                                                aGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGAGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCCTTACCATCCCGTCgctgcgat     >  2:480995/1‑136 (MQ=255)
                                                                                                                                 ggCAATCATAGTGGAACACGCGGTCTTGCGGAATCAGCAATCGAGCAGTACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTc    <  2:569857/143‑1 (MQ=255)
                                                                                                                                  ggcaTCATAGTAGAACAGGCAGTCGTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGTCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCg   <  2:425079/140‑1 (MQ=255)
                                                                                                                                  gCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCg   >  1:850950/1‑143 (MQ=255)
                                                                                                                                   cagtCATAGGAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGc  <  1:260351/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TTTGAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGC  >  NC_000913/711396‑711669

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