Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,505 G→T 6.4% E29* (GAG→TAG) ‡ ygiC → putative enzyme with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,5050GT 6.4% 174.4 / 3.3 63E29* (GAG→TAG) ‡ygiCputative enzyme with ATPase activity
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

ACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAA  >  NC_000913/3180365‑3180645
                                                                                                                                            |                                                                                                                                             
aCTATGCAGCGTAGTGCACCCGGTTCCTCTTCTCGTCCAACGTGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAg                                                                                                                                             <  1:420230/143‑1 (MQ=255)
     gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTa                                                                                                                                        <  1:1026793/143‑1 (MQ=255)
      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                       >  2:241595/1‑143 (MQ=255)
      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                       <  1:429365/143‑1 (MQ=255)
          tcaGTGCAACCGGTCCCCCTTCTGGTTCAATGGGTAGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAAGACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                  <  2:1091662/141‑6 (MQ=255)
          gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                  <  1:861441/143‑1 (MQ=255)
              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACGGgtgggc                                                                                                                              >  2:687643/1‑140 (MQ=255)
              tGCAACCGGGACTGATTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGACAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAACACGGTCTCCATTTTACCACCATGTACGGCGAGCCGTACTGGTGTGa                                                                                                                              <  1:484432/143‑1 (MQ=255)
                cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAg                                                                                                                            >  1:739111/1‑143 (MQ=255)
                  aCCGGTACCTCTTCTCGCTCATAGGCTAACTGATACCGATGGATATAGTGATTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGGTGTGAAAAt                                                                                                                          <  1:1138357/143‑1 (MQ=21)
                   ccGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCCGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTACAATAATCAAACCATGTACGGCGAGCCGTACTGGTGTGAAGATg                                                                                                                         >  1:331235/1‑143 (MQ=255)
                            cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACa                                                                                                                >  1:320885/1‑143 (MQ=255)
                              tctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAg                                                                                                              <  1:77195/143‑1 (MQ=255)
                              cctcGTCCAACGGGTGGTTGATACCGAGGGAAAGAGTCAGTATTACCGAGCGCCCGGACGGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAg                                                                                                              <  2:474592/142‑1 (MQ=255)
                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCGCGGACTGGCGGGAGAAAGCCCACGAATACGGTTTCAAGTTTCACACCATGTACACCGAGCCGGACTGGTGTGAACATGCAGACTACAAGATGa                                                                                                          >  2:55682/1‑143 (MQ=255)
                                   ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGCTGCTGCTTCCTACACGTTGAc                                                                                                         >  1:437191/1‑143 (MQ=255)
                                   ttCAATGGGTGGCTGATACCGATGGAAAGAGTAAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGCGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                         >  1:18770/1‑143 (MQ=255)
                                   ttCAAGGGGTGGCGGATACCGATGGAAAGGGGCATGATTACCGGGCGCTGCGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                         <  2:835294/143‑1 (MQ=255)
                                         gggTGGCATATATCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATTCGGTTTCAATGTTCACCCCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGc                                                                                                   <  2:242445/143‑1 (MQ=255)
                                          ggTGGCTGATACCGAGGGAAAGGGTCAGTATTACCGGGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTGCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGcc                                                                                                  <  2:859841/143‑1 (MQ=255)
                                           gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTAAGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTGGAACCACGccc                                                                                                 >  2:486468/1‑143 (MQ=255)
                                            tGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCGGAGCTGGCGATAAGATGCTTACTACAAGTTGACCCTCGCCCa                                                                                                >  2:263572/1‑143 (MQ=255)
                                             ggCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAg                                                                                               <  1:1242099/143‑1 (MQ=255)
                                             ggCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAg                                                                                               <  1:571585/143‑1 (MQ=255)
                                               cTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGACCCTTCATACAAGCTGACCCTCGCCCACGt                                                                                             >  2:953492/1‑143 (MQ=255)
                                                   tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaa                                                                                         <  1:356321/143‑1 (MQ=255)
                                                    aCCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaaa                                                                                        >  1:478889/1‑143 (MQ=255)
                                                        aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAGGAGGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                    <  2:146941/143‑1 (MQ=255)
                                                        aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACCACTATGTCGTCAATCTCCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                    <  1:1271507/143‑1 (MQ=255)
                                                          ggAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                  <  1:1279331/143‑1 (MQ=255)
                                                           gacagagTCGTTGTTACCGAGCGCCCGGACTGGCGTGCGATAGCCCACTAATACGGTTTCAATTTTCACACGATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGga                                                                                 <  2:1131480/140‑1 (MQ=255)
                                                            aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCCATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTCCAAGTTGACCCTCGCCCAGCTTGAACAGCTGgaa                                                                                >  1:223227/1‑143 (MQ=255)
                                                            aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCGGGTTGAAAAGCTGgag                                                                                >  1:1016277/1‑142 (MQ=255)
                                                            aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                <  2:974641/143‑1 (MQ=255)
                                                               aggTCAGTATTACCGAGCGCGCGGACTGGCGTGAGAATGCCTACGAATATTGTTTCAATTTTCATACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                             <  1:376234/141‑1 (MQ=255)
                                                                agtaagtaTTATCGAGCGTCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAGTGGTGGGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGaagaag                                                                            <  1:318270/143‑1 (MQ=255)
                                                                      tgttcccGAGCGCCCGGACGGGCGTGAGGAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccg                                                                      <  2:860366/138‑1 (MQ=255)
                                                                        ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                    >  1:1000226/1‑143 (MQ=255)
                                                                                   ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGACAAGCTGGAAGAAGGCCCCGCCGCACTGCACCa                                                         >  2:189614/1‑143 (MQ=255)
                                                                                       aCTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                     <  1:1032952/143‑1 (MQ=255)
                                                                                         tGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgtg                                                   >  1:738585/1‑143 (MQ=255)
                                                                                           gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGCAGATGCTTACTACAAGTTCACCCTCGCCCAGGATGAAAAGCTGGCATAAGTCACCGCCGAGCTGCACCAGCTGTGcc                                                 >  2:1157018/1‑143 (MQ=255)
                                                                                                 gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTCCAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGCAGCCAACGCCGAACTGCCCCAGATGTGCCTGACAg                                           >  2:293361/1‑143 (MQ=255)
                                                                                                    aGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGCACTGCACCAGATGTGCCAGAAAgagg                                        >  2:595453/1‑143 (MQ=255)
                                                                                                     gCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                       <  1:1025544/143‑1 (MQ=255)
                                                                                                     gCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTCGCCCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                       <  2:965903/143‑1 (MQ=255)
                                                                                                      cccACGAATACGGTTTCAATTTCCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAAGTTGCAGATCTTGAAGAAGGCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                      <  2:273262/143‑1 (MQ=255)
                                                                                                          cGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAATAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaa                                  <  1:666204/143‑1 (MQ=255)
                                                                                                           gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                                 <  1:486468/143‑1 (MQ=255)
                                                                                                            aaTACGGTTTAAAGTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCGCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                <  1:528752/143‑1 (MQ=255)
                                                                                                              tACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAAATGAGAACTACAAGTTGACCCTGGACTAGGTGGAAAAGCGGGAAGAAGTCACCATCTAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                              >  1:474923/1‑143 (MQ=255)
                                                                                                                       aaTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAACGCCAg                     >  2:943023/1‑143 (MQ=255)
                                                                                                                       aaTTTGCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACACCCGCCCAGGGTGAAAAGCTGGCAGAAGTCACCGCCGACCTGCGCCAGATGTGCCCGAAAGTGGGGGAAAAAGGGATCGCCAg                     >  2:943030/1‑143 (MQ=255)
                                                                                                                        aTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAATGATCGccccc                    >  1:1243080/1‑140 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAACGTGGTGGAAACAGTGATCGCCCGCGAtaaact            >  1:228376/1‑138 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt            >  1:271484/1‑143 (MQ=255)
                                                                                                                                 accagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATGGCCAGCGATGAGCt            >  2:55412/5‑143 (MQ=255)
                                                                                                                                  cATGTACGGCGGGCCGTACTGGTGCGACGATGCTTACCACACGCTGCCCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga          <  1:498965/143‑1 (MQ=255)
                                                                                                                                   gtGTACGGCGTGGGGTGCAAGTTTGAAGGTGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat         <  1:526818/142‑1 (MQ=255)
                                                                                                                                   cagTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga          >  1:267926/3‑143 (MQ=255)
                                                                                                                                      tACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc      >  1:1248810/1‑143 (MQ=255)
                                                                                                                                       aCGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGACGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGGTGAGCTGATGAcc     >  2:804993/1‑143 (MQ=255)
                                                                                                                                          gCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaaa  <  1:1030181/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAA  >  NC_000913/3180365‑3180645

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