Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,941,734 Δ1 bp noncoding (409/2906 nt) rrlC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,941,7340A.99.2% 411.7 / ‑0.8 235noncoding (409/2906 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (127/106);  minor base G (1/1);  total (128/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 255 total aligned reads displayed.

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTA  >  NC_000913/3941644‑3941786
                                                                                          |                                                    
gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                    >  1:3577829/1‑93 (MQ=21)
gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                    >  2:6444706/1‑93 (MQ=25)
gATACAGGGTGACAGCACCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                    >  1:398354/1‑93 (MQ=16)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                   >  2:2616117/1‑92 (MQ=21)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                   >  1:6305333/1‑92 (MQ=18)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                   <  1:7446828/93‑2 (MQ=17)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                   >  2:1087133/1‑92 (MQ=21)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                   >  2:6565787/1‑92 (MQ=21)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                 >  2:2503341/1‑93 (MQ=16)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                 >  2:2712271/1‑93 (MQ=12)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                 <  2:5894230/93‑1 (MQ=9)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                 >  2:3972175/1‑93 (MQ=9)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc                                                <  2:6635736/93‑1 (MQ=9)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc                                                >  2:4442852/1‑93 (MQ=16)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAcc                                                >  1:1522091/1‑93 (MQ=16)
     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                              >  2:6250292/1‑93 (MQ=16)
     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                              >  2:3810500/1‑93 (MQ=16)
      gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                             >  2:5486070/1‑93 (MQ=16)
      gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                             >  2:6150748/1‑93 (MQ=16)
       ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                            <  2:1989840/93‑1 (MQ=9)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa                                        <  1:783059/93‑1 (MQ=9)
            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                       <  1:5323300/93‑1 (MQ=9)
            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                       <  1:3788105/93‑1 (MQ=9)
            cAGCCCCGTACACAAAAATGCACATATTGGGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                       >  2:1375452/1‑93 (MQ=9)
            cAGCCCCGTACACAAAAATGCACATATTATGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                       >  1:2710312/1‑93 (MQ=255)
             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGGCACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                      <  2:177313/93‑1 (MQ=9)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                     <  2:4048657/93‑1 (MQ=9)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                     >  1:165410/1‑93 (MQ=12)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                     >  1:3499381/1‑93 (MQ=16)
               ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                    <  1:6056059/93‑1 (MQ=9)
               ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                    <  2:6098872/93‑1 (MQ=9)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                   >  2:3134735/1‑93 (MQ=16)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                   >  1:101453/1‑93 (MQ=16)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                   >  1:2085970/1‑93 (MQ=16)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                   <  1:4930458/93‑1 (MQ=9)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                   >  1:540082/1‑93 (MQ=16)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGTACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                  >  2:6212234/1‑93 (MQ=255)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                  >  1:5753959/1‑93 (MQ=9)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                  >  1:2802417/1‑93 (MQ=16)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                  >  1:1611139/1‑93 (MQ=16)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                 <  2:387217/93‑1 (MQ=9)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                 >  1:4746675/1‑93 (MQ=9)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                 <  1:221380/93‑1 (MQ=9)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                <  1:4315250/93‑1 (MQ=9)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                >  1:43561/1‑93 (MQ=12)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                >  1:1764093/1‑93 (MQ=16)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                >  2:2174594/1‑93 (MQ=16)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                <  1:2458717/93‑1 (MQ=9)
                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                <  2:2932811/93‑1 (MQ=9)
                   gTACACAAAAATGCACATATTGTGAGCTAGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                <  1:1637032/93‑1 (MQ=9)
                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                               <  1:6876306/93‑1 (MQ=9)
                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                               >  1:2561877/1‑93 (MQ=16)
                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                               >  2:3636382/1‑93 (MQ=16)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                              <  2:953711/93‑1 (MQ=9)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                              >  1:1051648/1‑93 (MQ=16)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                              <  1:6360297/93‑1 (MQ=9)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                              <  1:3379370/93‑1 (MQ=9)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                              >  1:1386924/1‑93 (MQ=9)
                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                             <  1:5643837/93‑1 (MQ=9)
                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                             <  1:2963207/93‑1 (MQ=9)
                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                             <  2:2055222/93‑1 (MQ=9)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                            >  2:2407939/1‑93 (MQ=9)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                            >  2:5352742/1‑93 (MQ=16)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                            >  2:5276254/1‑93 (MQ=16)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                            >  2:462851/1‑93 (MQ=16)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAAACt                            >  2:3039247/1‑93 (MQ=255)
                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                           <  2:3405029/93‑1 (MQ=9)
                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                           <  2:1658176/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          >  1:1957216/1‑93 (MQ=16)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          >  1:1275513/1‑93 (MQ=16)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          <  2:4903465/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          >  1:5480446/1‑93 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          <  2:4600411/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          >  2:3223162/1‑93 (MQ=16)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          <  2:1458820/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                          <  1:6617238/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGACCCATCCTCCAAGGCTAAATACTcc                          >  1:6031504/1‑93 (MQ=11)
                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                         <  1:4209591/93‑1 (MQ=9)
                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                         >  2:6835326/1‑93 (MQ=16)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                       >  2:2753939/1‑93 (MQ=16)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                       >  2:6276343/1‑93 (MQ=12)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                       <  1:3510097/93‑1 (MQ=9)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                      <  1:355060/93‑1 (MQ=9)
                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                     >  2:1898320/1‑93 (MQ=16)
                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                     <  2:6953185/93‑1 (MQ=9)
                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                     >  1:733867/1‑93 (MQ=12)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  2:5873926/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  >  2:5613585/1‑93 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  2:5934469/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  1:2331439/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  1:7126987/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  2:53683/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  1:6187554/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                  <  2:4503210/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCCAAATACTCCtgactgac                  >  1:6153525/1‑93 (MQ=9)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                 <  2:3061784/93‑1 (MQ=0)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                 >  2:5124397/1‑93 (MQ=0)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                 >  1:1154235/1‑93 (MQ=0)
                                       tgtgAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATa             >  1:6961631/1‑93 (MQ=18)
                                                  tGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTa  <  1:3331082/93‑1 (MQ=17)
                                                                                          |                                                    
GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTA  >  NC_000913/3941644‑3941786

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